Abstract

GMD

Geoscientific Model Development

GMD

Geosci. Model Dev.

1991-9603

Copernicus Publications

Göttingen, Germany

10.5194/gmd-19-6829-2026

copan:LPJmL: a new hybrid modelling framework for dynamic land use and agricultural management

copan:LPJmL: a hybrid land use and management framework

Breier

Jannes

jannes.breier@pik-potsdam.de Schwarz

Luana

Prawitz

Hannah

von Bloh

Werner

Müller

Christoph

https://orcid.org/0000-0002-9491-3550

Wirth

Stephen Björn

https://orcid.org/0000-0003-3090-3318

Bechthold

Max

Gerten

Dieter

https://orcid.org/0000-0002-6214-6991

Donges

Jonathan F.

https://orcid.org/0000-0001-5233-7703

1Potsdam Institute for Climate Impact Research (PIK), Member of the Leibniz Association, P.O. Box 60 12 03, 14412 Potsdam, Germany 2Humboldt Universität zu Berlin, Department of Geography, Unter den Linden 6, 10099 Berlin, Germany 3Integrative Earth System Science, Max Planck Institute of Geoanthropology, Kahlaische Str. 10, 07745 Jena, Germany 4University of Potsdam, Institute of Environmental Science and Geography, Karl-Liebknecht-Str. 24–25, 14476 Potsdam, Germany 5University of Potsdam, Institute of Physics and Astronomy, Karl-Liebknecht-Str. 24–25, 14476 Potsdam, Germany 6Stockholm Resilience Centre, Stockholm University, Kräftriket 2B, 106 91 Stockholm, Sweden

Jannes Breier (jannes.breier@pik-potsdam.de)

28July2026

19 14 68296855 11September2025 28October2025 21May2026 3June2026

2026

This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/

This article is available from https://gmd.copernicus.org/articles/19/6829/2026/gmd-19-6829-2026.html

The full text article is available as a PDF file from https://gmd.copernicus.org/articles/19/6829/2026/gmd-19-6829-2026.pdf

Abstract

Dynamic Global Vegetation Models (DGVMs) are established in environmental and agricultural sciences for many purposes, e.g., modelling plant growth and productivity, water and carbon cycles, and biosphere-climate interactions. Nevertheless, DGVMs are still rather limited in terms of simulating mutual interactions between biospheric and human processes. While DGVMs such as the Lund Potsdam Jena managed Land (LPJmL) model have been successfully connected to Integrated Assessment Models (IAMs), the model couplings often remain loose and static over the simulation period. The copan:LPJmL modelling framework is an extension of the copan:CORE framework for integrated and dynamic human-Earth system modelling, and addresses this issue by integrating LPJmL via a new interface, consisting of an LPJmL coupling library and a Python library pycoupler, which together enable LPJmL inputs and outputs to be coupled in copan:LPJmL during the simulation period. It uses the copan:CORE entities and integrates the coupled data into the World (simulation space as a whole) and the (grid) Cell entity, allowing other entities such as Individuals, e.g., for agent-based modelling (ABM), to access them. Besides ABM, this framework allows for a broad range of modelling approaches to be represented with copan:LPJmL, of which we introduce three examples: (1) The model of Integrated Social-Ecological Resilient Land Systems (InSEEDS), which uses a classical ABM approach to model management decisions by farmers, (2) an adaption of an established crop calendar model (Crop Calendar), and (3) a novel Large Language Model (LLM)-driven ABM approach (LLM Fertilization).

HORIZON EUROPE Framework Programme

101081661

Carl-Zeiss-Stiftung

n/a

European Research Council

ERC-2016-ADG-743080

Bundesministerium für Forschung, Technologie und Raumfahrt

01LS2001A

1Introduction

The Anthropocene is a new epoch in the evolution of the Earth system characterized by increasingly strong and entangled coevolutionary dynamics of biogeophysical Earth system processes with human societies . The recognition of biogeophysical and social processes as intertwined has emerged from historical developments in the field of Earth System Science. While James Lovelock already published the Gaia hypothesis in the early 1970s , the broader scientific uptake of the term “Earth System Science” took place in the 1980s and 1990s . Central publications in this decade, such as the Bretherton Diagram and the Brundtland report , acknowledge that human societies are tightly connected to Earth system dynamics. Followed by the realisation that (some) human societies now act as the major driving force of change on our planet , while at the same time being shaped and impacted by the ecological conditions they are embedded in , diverse novel conceptualisations of this intertwinedness, like the “technosphere” have emerged.

The consideration of the role of humans in the Earth System has thus progressed towards a coevolutionary, bidirectional approach, now evident in different strands of simulation modelling that can be summarised under the term World-Earth (System) modelling (WEM) or integrated human-Earth system modelling, with the entirety of human societies being referred to as the World . WEM can be understood as a specific strand of social–ecological systems (SES) modelling that explicitly focuses on representing bidirectional feedbacks between human societies and the Earth system at regional to global scales . There is growing research calling for and adopting a social-ecological, integrated perspective of humans embedded in the Earth System (e.g., . Different modelling communities approach this integration in distinct ways, and focus on different aspects of social-ecological dynamics. One central challenge of these efforts is moving beyond simple proof-of-concept models towards more complex, integrated models .

Land systems are a key example of the inherently coevolutionary nature of social-ecological interactions in World-Earth Systems . Agricultural production is the single largest driver of transgressions of multiple planetary boundaries: land-system change, freshwater use, biogeochemical flows and biosphere integrity , also contributing to the transgression of other planetary boundaries, such as climate change. In socio-ecological terms, global food demand drives responses in land use and production, including cropland and pasture expansion, and intensification through management practices such as irrigation and fertilization, and the reorganization of supply chains through trade. On the production side, farmers' management decisions, such as the use of crop rotations, cover cropping, reduced tillage, and integrated pest management can significantly influence the adverse environmental impacts that manifest themselves in the transgression of planetary boundaries . Improved knowledge about the dynamic and coevolutionary development of agricultural systems at spatial scales from local to global is therefore imperative to foster our understanding of fundamental Anthropocene dynamics.

Land use systems alone have been a research subject in Earth system science for years, leading to their incorporation into DGVMs (such as LPJmL), crop models (e.g., DSSAT) as well as Earth System Models (ESMs, (e.g. MPI-ESM)) . Specific (such as tillage) or bundled management practices (like conservation agriculture) have been represented in detail and studied in local and global applications of such models . Additionally, more holistic applications have been carried out, e.g., to show the potential of maintaining multiple planetary boundaries while feeding 10 billion people . However, the most important factor in implementing such measures or achieving such goals has been neglected: Humans are often merely represented as biophysical “disturbance” factors. This way, their decision-making capacities and social-psychological and socio-cultural complexities are therefore underrepresented, with many assumptions and dynamics, for example on land management and land use, remaining rather static . IAMs make it possible to simulate and optimize certain land use dimensions under a set of “human boundary conditions”, such as Shared Socioeconomic Pathways (SSP) scenarios . Through the application of macroeconomic and energy-economic optimization approaches, many IAMs reduce the complexity of human behaviour to predefined scenarios or boundary conditions, thereby excluding key dimensions of human decision-making uncertainty . While the SSP1 and SSP2 scenarios presented in the Sixth IPCC Assessment Report are considered economically feasible within IAMs, their assumptions regarding social-ecological feasibility – such as behavioural adaptation and governance dynamics – remain subject to significant uncertainty . Furthermore, these scenarios lack the fully coupled feedbacks of societies with the Earth system, i.e. the coevolution inherent in the underlying dynamics that is currently not representable endogenously with the majority of existing modelling approaches .

To address this gap, introduced the copan:CORE modelling framework, which supports the development of World-Earth models (WEMs). WEMs are characterised by the explicit and bidirectional coupling of social and Earth system processes, enabling the coevolution of human and natural systems to be represented within a single modelling framework. They aim to go beyond traditional IAMs by incorporating human agency, social heterogeneity, and feedbacks between human decisions and biophysical dynamics . A modelling framework is a software or conceptual infrastructure that provides the structure, components, and tools needed to build and run models, without being a single model itself. Copan:CORE provides the components and structures required to build WEMs. It consists of entities, such as a simulation cell or individuals, that interact with each other via various processes . The latter are categorized by three overlapping process taxa, representing biogeophysical and biogeochemical (ENV, e.g. biophysical conditions, crop growth), socio-metabolic (MET, e.g. crop harvest, fertilization) and socio-cultural processes (CUL, e.g., governance, social learning, social norms dynamics, or individual cognitive-behavioural processes like attitude formation) . We here have advanced the copan:CORE framework by integrating LPJmL as the ENV taxon , through which we can represent diverse natural as well as managed land systems, such as forests, grassland, and cropland and link them bidirectionally to social processes that can be of economic nature but also go beyond that, for example including behavioural change, collective decision-making or political processes. Based on this, MET and CUL-based components of different modelling approaches can be realised by applying copan:LPJmL to represent various dimensions of social-ecological systems. By building on the diverse entities and process/model templates of copan:CORE, complex and detailed systems can be integrated to map underlying dynamics such as human, social, and societal processes using ABMs, dynamical system, or rule-based approaches.

In this paper, we introduce the copan:LPJmL modelling framework as a flexible and extensible platform for the development of WEMs, including a process-based DGVM Earth in a land system context. To demonstrate its applicability and versatility, we present three example applications:

InSEEDS, a global agent-based framework for simulating the adoption of regenerative agricultural practices by farmer agents ;

Crop Calendar, an integrated (in runtime calculation) update of the rule-based model published by ; and

LLM fertilization, representing the application of crop fertilizer by LLM farmer agents. These examples illustrate how copan:LPJmL enables the integration of diverse socio-cultural and socio-metabolic processes with detailed Earth system dynamics.

2Framework description

copan:LPJmL is a newly developed, enhanced WEM framework that integrates the DGVM LPJmL and the existing open modelling framework copan:CORE, thereby creating a modelling environment in which biophysical Earth-system and social world processes can be represented and bidirectionally linked to enable fully integrated and global gridded social-ecological modelling. It is designed to bring together the details and complexity of a process-based land biosphere model with the flexibility of an open modelling framework, thereby enabling the build-up of new types of integrated models. It allows selection from various modelling approaches such as agent-based, rule-based, dynamical system, or data-driven (statistical). copan:CORE forms the core framework in which the modular structure and building blocks, such as basic entities, are defined that can interact with each other. In copan:LPJmL, LPJmL is currently integrated via annual coupling to represent processes of the ENV taxon, with future development aiming to implement more frequent coupling intervals, ranging from monthly to daily. This is facilitated by LPJmL itself, which has been extended with a coupler library that is available in LPJmL from version 5.6 onwards, and applied in copan:LPJmL via the Python interface pycoupler .

Figure 1

Basic scheme of copan:LPJmL featuring the global grid, the process taxonomy as well as main entities with the three application examples: InSEEDS, Crop Calendar, and LLM Fertilization.

2.1The CORE framework

copan:CORE, introduced by , has established a novel systematic approach for building WEMs in a Python-based modelling environment with its corresponding Python library pycopancore. All copan:CORE model elements are represented as modular, object-oriented entities. These include agents (such as households, firms, or governments), social institutions, components of the Earth system or even more abstract entities such as a rulebook (basis for Crop Calendar). Additional entities can be flexibly introduced using the framework. In copan:CORE, each of these entities is involved in model processes that can be categorized by applying the aforementioned taxonomy system described by (ENV, MET and CUL). Another essential feature of copan:CORE is its support for heterogeneous agent populations and scalable network structures that define the relationships and interactions among agents, referred to as Individuals, and between agents and their linked entities. Besides Individuals these available entities in copan:CORE are the simulation space, the World, the underlying elemental spatial units, the Cells, the Social Systems such as countries or cities . Recently, added an additional Group entity as a social structure in which Individuals can organize themselves. In contrast to social systems, groups are not tied to a specific location, are loosely connected, and share only certain (temporary) characteristics, such as social norms. For copan:LPJmL, this variety is especially valuable, as the underlying hierarchies and networks provide a basic, real-life reflecting structure that aligns well with the simulation space (world) of LPJmL, in which Cells also reflect the elementary spatial unit of simulation. This way, the framework can represent the LPJmL organisational structure within the copan:CORE World. The library, pycopancore, also included model components, models and studies outside of the CORE definition . Within the process of building new software around copan:CORE, such as copan:LPJmL, this structure was reorganised to distinguish the CORE from other software and components related to copan:CORE, leaving only the components used by as an exemplary modelling approach .

2.2LPJmL as a model component

The DGVM LPJmL simulates the carbon, nitrogen, water and energy cycles of the terrestrial biosphere in coupling with the growth and productivity of natural and agricultural ecosystems, forced by climate, crop distribution, and other globally gridded input data . It can represent various agricultural management practices (see Sect. ) that requires additional information on the spatial and temporal distribution of management systems . The smallest entity representing these processes is a grid cell, by default with a spatial extent of 0.5×0.5°, such that 67 420 cells represent the land surface globally. This spatial resolution reflects the design of LPJmL as a global-scale model that aims to capture large-scale biophysical dynamics while maintaining computational tractability. At this scale, individual grid cells represent aggregated socio-ecological systems rather than single land holdings or individual actors. In the current example applications, one representative agent is associated with each grid cell, which can be interpreted as representing the aggregated behaviour of multiple real-world actors within that region. However, the copan:LPJmL framework does not restrict the number of agents per grid cell. Applied models can introduce multiple agents per cell, represent heterogeneous actor types, or implement social interaction networks within and across cells. In this way, the spatial resolution of the biophysical LPJmL model can be combined with finer-grained representations of social processes depending on the research question and data availability. Plants are modelled according to the concept of Plant Functional Types (PFT) . The same principle is applied to Crop Functional Types (CFT) used for modelling agriculture on prescribed shares of each cell, called stands . copan:LPJmL introduces LPJmL as a new model component that covers the ENV taxon. It represents the LPJmL model as a World entity WLPJmL (LPJmL simulation space, i.e. global), that integrates the coupled LPJmL model 1WLPJmL=(LPJmL,{Ci∣i=1,2,…,n}) and its corresponding Cell entities CLPJmL as shown in Fig. following the concept of . CLPJmL includes the input and output of LPJmL for the cell as well as the grid information (longitude, latitude) and a network of its neighbouring cells N(Ci): 2Ci=(input,output,grid,N(Ci)), where N(Ci) is defined as 3N(Ci)={Cj∣d(Ci,Cj)≤r,j≠i and Cj∈WLPJmL} and is calculated based on the grid information of each cell, with d(Ci,Cj) denoting the Euclidean distance between cell centroids, and r=2⋅0.5° corresponding to the maximum distance between adjacent cells in the 0.5° grid (including diagonal neighbours). In the copan:LPJmL implementation, the World entity stores global LPJmL input, output, and grid information, while the associated Cell entities provide access to the corresponding spatial subsets of these data. This structure allows processes to operate either at the global (World) level or at the level of individual spatial units (Cell), enabling both aggregated and local interactions between biophysical and social model components. The specific technical implementation of these data structures and their efficient handling is documented in the copan:LPJmL software repository .

Figure 2

Stylised flow chart of the integrated copan:LPJmL WEM framework. The orange highlighting indicates the parts and attributes that can be modified, extended or added. This applies to each entity type represented within its structural hierarchy and modelling interfaces of the framework described in detail in Sect. . “A WEM model” can be implemented by using and/or extending the World and corresponding Cell entities and by creating, for example, Individual A and B with distinct features and processes (a). The entities are hierarchically linked, while their interactions and processes are defined within the respective coupled model implementation. The three dots indicate that entities and attributes can be further extended. By creating A_WEM_model_config.yaml, those can be set up and added to the LPJmL configuration (LPJmLConfig) in which further simulation parameters and settings are specified (b). The coupled simulation is started based on the generated <sim_name>_config.json file (c). The LPJmLCoupler handles the coupling between LPJmL (Earth) and “A WEM model” by receiving and sending coupled data on an annual basis. The circular arrows represent the annual exchange of coupled LPJmL input and output data. “A WEM model” can compute intra-annual processes based on coupled LPJmL output to generate new coupled LPJmL input for closed feedback loops throughout the simulation period.

2.3Bidirectional coupling with LPJmL

To enable a bidirectional, annual data coupling, LPJmL has been extended such that input and output data can be exchanged via TCP/IP socket connections during runtime. This option is now available for all inputs and outputs in the LPJmL configuration. During preprocessing, this configuration incorporates additional model settings such as the specification of coupled inputs and outputs. An excerpt of such a configuration with coupled inputs and outputs is shown as in Listing , also containing additional information such as the name of the coupled model, the start year of coupling, as well as the hostname and the port number of the coupled model. The coupling itself is provided by the LPJmL internal coupler library, which was implemented to encapsulate any socket communication in LPJmL. The coupling follows a prescribed protocol, as described in Listing of the Appendix.

Conversely, to enable simple and practical coupling with LPJmL, we have developed pycoupler, a Python interface that provides tools around the coupling as well as handling of LPJmL in Python environments (Fig. ). It serves as a configuration interface for coupled as well as stand-alone LPJmL simulations as well as functionality to send and receive xarray-based data objects from LPJmL. The coupling follows the prescribed coupling protocol, for which pycoupler contains the module LPJmLCoupler. The protocol is provided in detail in Listing of the Appendix. With LPJmLCoupler, LPJmL is coupled on the Python side, and both sending and receiving methods are provided. It creates the underlying data structures (input, output, grid, country, etc.) and in the case of input data, initialises it with data read in at the start of the simulation. While all static data (e.g., grid, country) are only sent at the start of the simulation, time transient LPJmL output data is received from LPJmL on an annual basis and updated during simulation time via pycoupler (Sect. ). Likewise, the input data can be changed within copan:LPJmL and sent back to LPJmL vice versa. pycoupler also allows for configuring (coupled) LPJmL simulations by providing an LPJmLConfig module, which simplifies reading, handling, modifying and writing of LPJmL configurations, the files of which are subsequently used to run the corresponding LPJmL simulations following . Sending and receiving data is possible during the runtime of LPJmL on an annual basis. After a defined coupling year, LPJmL waits until the required input data is sent via LPJmLCoupler to continue the simulation until the end of the year to send back the output data. For the following years, the procedure is repeated. For a coupled simulation, the modules are usually applied together and in sequence in an individual run script. This way pycoupler supports a wide range of simulation and coupling setups. The pycoupler package is available on GitHub with more detailed function documentation and examples . Further technical details of the coupling implementation are documented in the software repository. Beyond that, pycoupler serves further utility functions around LPJmL, such as getting the neighbour cells of a cell or subsetting the grid for country-specific LPJmL simulations.

2.4Development of WEMs with copan:LPJmL

The copan:LPJmL framework with its corresponding pycopanlpjml library unites all presented software components pycopancore (Sect. ), LPJmL (Sect. ) including the coupler library and the Python interface pycoupler (Sect. ). For model developers, copan:LPJmL provides the hierarchical structure of entities illustrated in Fig. . The framework defines the global World entity and its associated Cell entities as the spatial unit structure and interface to LPJmL. Additional entities such as Individuals (e.g. farmer agents in an ABM) can be defined and linked to these spatial units. Model logic is implemented through initialisation and update methods that define how entities interact and evolve. Concrete application examples for this are presented in Sect. . All model parameters and settings can be specified via the corresponding model configuration files, which are combined with the LPJmL configuration to form a complete simulation setup, shown in Listing . Each configuration defines a specific simulation experiment that can be executed using pycoupler, enabling the coupled execution of the specific model that is implemented in the copan:LPJmL framework and LPJmL (Listing ). With this structure, copan:LPJmL provides a modelling framework that enables researchers from diverse disciplinary fields to engage in modelling human-environmental interactions in the land-use sector. It allows users to specify custom decision rules, agent interactions, or feedback mechanisms, without requiring modification of the LPJmL model. This lowers the entry barrier and fosters interoperability, supporting a wide range of applications and interdisciplinary collaboration.

Within copan:LPJmL, LPJmL, the DGVM, defines the spatial scope and spatial resolution in this setup. It constitutes the ENV taxon of the WEM, and is complemented by the taxa MET and optionally CUL (Fig. ). While MET processes are typically represented in many integrated modelling approches (i.e., IAMs), the CUL taxon is, as indicated in Sect. , often neglected . We here provide a platform to include and explicitly represent these social processes of different actors and their interactions that influence environmental dynamics, for example, through land-use and agricultural management decisions that affect the inputs to the ENV component and its feedback to the CUL taxon. It facilitates the representation of “the social” from different perspectives. The CUL taxon can cover individual decision-making, social learning, but also economic dynamics, policy-making, and much more. A broader discussion of the types of social and decision-making processes that can be represented within this framework is provided in Sect. . Within the copan:LPJmL framework, MET and CUL processes can be represented with the full flexibility of copan:CORE. Their potential interactions with ENV are bound to the data interface of the terrestrial Earth, the inputs and outputs of LPJmL (Sect. ). These are, on the one hand, the exposed inputs, represented in the World and Cell entity as an input attribute. In general, all available time-variable inputs in LPJmL can be exposed via pycoupler and thus coupled, an overview is given in Table .

Table 1

Overview and categorization of key inputs to LPJmL (ENV).

Category Key Variables Climate Temperature (daily mean, min, max), precipitation, radiation, atmospheric CO₂ concentration, wind speed Land Use Land-use composition (fraction of cultivated CFT), irrigation systems (for each CFT) Management Tillage systems, crop residue management, fertilizer & manure nitrogen, sowing dates & harvest dates Society Population density, country affiliation, water use of household, industry & livestock Miscellaneous Soil texture type, soil acidity, elevation, livestock density, lakes & reservoirs, river drainage direction, neighbour irrigation network, water demand for households, livestock and industry, atmospheric nitrogen deposition

A corresponding change to one of these inputs, which is applied in LPJmL in the following year, results in a corresponding biophysical response that is reflected in the outputs. Vice versa, the range of available LPJmL outputs can be utilized in the MET and CUL taxon via the output attribute in World and Cell to inform and shape the respective processes. Here, too, all outputs are available in copan:LPJmL, the most important of which are listed in Table .

Table 2

Overview and categorization of key LPJmL outputs (ENV).

Category Key Variables Vegetation Aboveground biomass, total vegetation carbon, net primary productivity (per PFT available), gross primary productivity, respiration, Soil & Litter Soil carbon & nitrogen content (per soil layer available) , soil temperature (per soil layer available), maximum thaw depth, N₂O emissions, nitrogen leaching, ammonium volatilization Water Soil evaporation, transpiration, canopy interception, local runoff, river discharge, potential evapotranspiration, soil moisture, irrigation water use per crop Fire & Disturbance Burned area fraction, fire emissions, fire carbon released, deforestation emissions Land use & Agriculture For each simulated crop: yield, growing area, sowing date, harvest date, fertilizer amount, fertilizer & application

In principle, inputs can be extended, for example, by converting further parameters or settings into inputs. Any state variable or flux can also be written as an output, providing even more opportunities for model coupling. By design, copan:LPJmL allows the integration of a broad range of modelling paradigms of different domains that represent socio-ecological dynamics and feedbacks at varying levels of complexity and abstraction. Table shows an overview of approaches that potentially could be represented in CUL and MET using copan:LPJmL. These include ABMs for simulating heterogeneous actors and interactions, rule-based systems for capturing institutional or behavioural heuristics, optimisation-based approaches for identifying efficient or goal-oriented management strategies, and surrogate or machine learning models for data-driven decision-making. System Dynamics models can be incorporated by formulating their underlying equation systems within the Dynamical Systems paradigm, ensuring compatibility with its feedback-oriented architecture. The flexibility enables the exploration of diverse real-world processes ranging from farmer decision-making to governance interventions.

Table 3

Modelling paradigms and application domains supported by copan:LPJmL, with representative examples.

Paradigm Application Domain(s) Representative Examples Agent-Based Models (ABMs) Socio-ecological systems such as land-use change and adaptation, or household and (farmer) decision-making; socio-hydrological systems, such as mitigation of extremes and water management; higher-level organisation, e.g., policy scenarios InSEEDS: heterogeneous farmer agents with LPJmL feedback ; CRAFTY: competition-based land allocation ; Household-level farming responses ; Social-hydrological studies Rule-Based Models (RBMs) Systems that follow rule-based logic, such as agricultural management and land use; higher-level organisations, e.g., acting schemes of governmental institutions Crop Calendar ; Land use allocation ; Crop rotation and fallow scheduling Optimization Models Economic land and resource optimization; yield and input efficiency; food security and environmental trade-offs; policy-guided scenario analysis Examples: MAgPIE, GLOBIOM, or IMAGE . E.g. for water and nutrient optimization or land-based mitigation Dynamical System Models System-level feedbacks and tipping dynamics; coupled socio-environmental processes; bioeconomic models; evolutionary game theoretic models Tipping point and regime shift models , Socio-epidemic models, complex contagion Human-climate models Data-Driven Models Emulators; Pattern discovery; Scenario validation and calibration; remote-sensing based approaches Statistical crop yield emulators ; Surrogate models ; Remote-sensing LLM-Enhanced Models Decision emulation; Stakeholder reasoning; Adaptive behavior LLM-based farmer management decisions LLM agents as institutional policy-makers

In addition, the framework supports a wide range of application areas that interact with the land system, in particular its ecological conditions, and resources. This can well result in more complex (social) structures in which, for example, only a first layer of agents interacts with the land system, while other actors are only connected to them, as in a supply and demand model for agricultural products involving farmers, food producers, and consumers. In general, potential domains are land-use change, agricultural management, climate adaptation, food-water-energy dynamics, policy evaluation, or socio-ecological transitions (Table ). Through the modular interface and configurable data exchange with LPJmL, copan:LPJmL can serve as a backbone for integrated assessments across spatial scales and decision contexts.

3Application examples

In order to illustrate a broad range of modelling approches that can be covered using copan:LPJmL, we present in the following three different representative applications examples of varying complexity, which are based on the different model paradigms and application areas shown in Table . While those examples give a taste of the possibilities of copan:LPJmL, they are not exhaustive, and many more applications of the framework are possible.

3.1InSEEDS: a new agent-based World-Earth model 3.1.1Modelling approach

In the field of SES science, different modelling approaches have emerged to capture the intertwined dynamics between human and biophysical spheres . Arguably, the most prominent approach to social-ecological modelling, ABM, is ideal to be used for modelling with copan:LPJmL . ABMs simulate interactions between agents and their environments over time. These agents can represent individuals, households, organizations, or other entities. The simulation of these micro-level interactions gives rise to different macro-level outcomes, like spatial adoption patterns of a certain land use . Several features of ABMs make them a particularly useful methodological choice for the investigation of SES. SES are often understood as complex adaptive systems and therefore are inherently characterized by dynamical adaptation to changing behaviors and environments . Furthermore, ABMs facilitate the study of macro-level phenomena emerging through micro-level dynamics. Lastly, ABMs are able to capture agent heterogeneity in human and biophysical spheres .

3.1.2Model design

InSEEDS, described in detail in and implemented in , is a WEM created using the copan:LPJmL framework that uses an ABM component to capture farmer agent decision-making (Fig. ). InSEEDS was originally designed to investigate the SES dynamics at play in transitions from conventional farming to regenerative farming practices such as conservation tillage. The CUL taxon comprises farmer management decision-making processes that are based purely on social interactions (i.e., evaluation of social norms) as well as information obtained from social-ecological processes (i.e., observing the environment). The social network of farmers forms the basis for these social dynamics and is initialised based on the LPJmL model grid (Sect. ). In this realisation, one representative farmer agent is assigned to each cell. This means that each farmer has a maximum of eight direct neighbours who form their neighbourhood, which is currently the only explicitly represented social network (Eq. ). However, this neighbourhood structure can be extended or complemented by additional non-local networks at regional or global scales. Such networks follow the implementation of acquaintance networks in copan:CORE and can, for example, represent groups of agents distributed globally or within specific regions. An example of an implementation of a social-ecological feedback for a farmer in InSEEDS is given in Fig. . It shows a simple attitude formation process of farmer agents as an evaluation of their farming performance. In this realisation, farming performance is measured using the average crop yield and soil carbon content at different points in time. It is assumed that farmers have perfect knowledge of both variables in order to enable a preferably deterministic behaviour for analysis purposes only. copan:LPJmL allows for a generic processing of any LPJmL output, like agricultural soil carbon or CFT-specific crop yield demonstrates the flexibility copan:LPJmL provides in combination with the multiple features of xarray in subsetting and aggregating the underlying data. Thereby, it enables direct feedback functions to be set up to simulate important social-ecological aspects, such as farmers' attitudes towards their land.

Figure 3

Detailed scheme of InSEEDS highlighted in Fig. representing the feedback mechanisms among the taxons ENV, MET and CUL. Individuals who are either of type pioneer or of type traditionalist observe crop yields and soil carbon and make decisions on conservation/conventional tillage to then be re-evaluated again.

Figure 4

Social-ecological processes and feedbacks in copan:LPJmL. (a) shows an example as a subset of InSEEDS – a farmer agent estimates their attitude toward their land based on two ENV variables (LPJmL output, see Sect. ). The agent compares the topsoil carbon content and average crop yield between simulation steps. (b) shows a similar process in Crop Calendar in which the calendar is informed by processed climate data to estimate the ideal sowing dates for the next year. Social and ecological processes are highlighted in red and blue.

Such processes are part of the MET taxon, constituting the cross-section of socio-cultural and biophysical processes. Figure a shows a typical example of an observation of ecological variables, as well as the initiation of potential management decisions: Information on average crop yield and topsoil carbon content is calculated by the LPJmL ENV component and provided via the output attribute in the cell to farmer agents within a MET process. Vice versa, farmers' decisions are forwarded to LPJmL as input for the simulation via the MET component. Following the agent-based logic, the current main actors in the InSEEDS model are individual Farmer agents. The Farmer agent class itself is a child class of the Individual agent class in the copan:CORE modelling framework, inheriting the logic of Individuals described in . As an additional property of the Farmer, we introduce two agent-functional types (AFTs) , a traditionalist and a pioneer farmer, who differ in their respective weighting of different parameters in the decision-making function (Fig. . The decision-making process is based on a formalisation of the Theory of Planned Behaviour originally described by . In principle, the decision-making process can be replaced by processes from other theories as well which is discussed in Sect. .

3.1.3Example simulation results

InSEEDS can simulate social-ecological model dynamics on a wide spectrum of scales up to global scale. Analysis possibilities include distributed and accumulated social and ecological outcomes of variables such as attitude, social norm, or soil carbon, crop yield, or even adoption patterns of certain management practices. InSEEDS is the first model realisation using copan:LPJmL that simulates the aforementioned coevolutionary social-ecological dynamics through closed feedback loops. Figure illustrates the underlying coevolution via three variables: The top row (a) shows the spatio-temporal dynamics of management practice adoption globally. The middle and lower row (b, c) depict the biophysical response of these behaviour changes in topsoil carbon content and average crop yield, compared to a business-as-usual simulation. In some areas, such as Kazakhstan, the spreading of conservation tillage and the underlying coevolution is particularly evident. Here, the adoption started in the southern regions and spread north-eastwards with moderate increases in soil carbon and significant increases in average crop yield. This synergistic effect results in a certain irreversibility in the modelled system. In our simulation, we find that in many regions, the adoption of conservation tillage has a positive effect on these variables, even though the results vary strongly at the local level. To better understand the underlying dynamics of specific cases, a more regional perspective is needed, which can be found together with the decision equations, parameterisation, parameter sensitivity, and detailed simulation results on coevolutionary model dynamics in .

The development of InSEEDS is currently in an early phase. Future projects such as mapping non-local networks, social systems with multiple layers of complexity, and more social-ecological feedback processes can build on this approach. copan:CORE as described by provides entities, networks and components to achieve this. At this point, at the latest, the connection to the IAMs and the existing paradigmatic problems as described in (Sect. ), such as missing closed social-ecological feedback loops can be drawn. By explicitly representing such closed social–ecological feedbacks, InSEEDS and the underlying copan:LPJmL framework provide a platform to explore the coevolution of social decision-making and environmental dynamics, offering new opportunities to investigate adaptive land-use transitions and management responses under changing environmental conditions.

Figure 5

InSEEDS simulation of (a) global spreading of conservation tillage represented as years since last management switch (backwards from simulation end year 2100). The adoption and spatial spreading of conservation tillage is depicted in orange, and conventional farming adoption in blue. If the saturation of both colours decreases over time, this indicates that there will be no further change in strategy. Vice versa, the more saturated the colour, the more recent the change. (b) and (c) show the corresponding changes in topsoil carbon content and average crop yield in 2100 compared to a simulation without management changes (business-as-usual) until 2100.

3.2Integrated rule-based approaches: Crop calendar 3.2.1Modelling approach

In many land-use and agricultural systems, management decisions are not made by autonomous agents, but instead follow fixed logics, institutional guidelines, or context-dependent thresholds. These RBMs operate via explicit if–then conditions, fuzzy logic or temporal schedules, allowing for transparent and interpretable decision structures . Their structured logic makes them particularly suitable for encoding expert knowledge, empirical heuristics, or scenario-specific governance interventions, especially when interactions between agents are minimal or absent. As such, RBMs offer an efficient and reproducible way to represent adaptive but non-agentive processes across socio-environmental domains. Nevertheless, the boundaries are fluid, and rule-based and agent-based systems often overlap. While global, top-down approaches with Boolean logic can be clearly assigned to RBMs, there are bottom-up, autonomous approaches, such as cellular automata, that can be located between these paradigms . Within copan:LPJmL, the RBM paradigm is particularly useful when system feedbacks should arise directly from dynamic ENV quantities – such as climate or resource availability – rather than from emergent behaviour. This makes them well-suited to simulate adaptive but non-agentic responses to environmental change.

Figure 6

Detailed scheme of the crop calendar implementation highlighted in Fig. and representing the feedback mechanisms between the taxons ENV and MET. Climate data is forwarded by LPJmL to be used by the Crop Calendar to determine sowing and harvest dates for the upcoming year. Under changing climatic conditions, choosing the right sowing and harvesting times is crucial for maximizing crop yields, indicated by the grey and orange maize plants.

3.2.2Model design

The adjustment of growing seasons in response to changes in climatic conditions is such an example and a central element in agricultural adaptation strategies. While changes in sowing dates may already be implemented by the farmers based on their experience , cultivar choices are subject to availability and breeding . proposed a modelling approach for simulating changes in growing seasons based on changes in climatic conditions only. This approach was used to create climate-scenario-specific time series of sowing dates and cultivar parameters as inputs for simulations with LPJmL . Implementing the algorithms of using copan:LPJmL allows for a flexible application of adaptive growing seasons during runtime without requiring to previously compute growing seasons and corresponding crop parameters for each climate scenario . In this setup, which is exemplary for rule-based management decisions, LPJmL passes only climate information as output to the growing season rules, which represent the MET taxon (Fig. ). A CUL taxon is not involved in this setup, as there is no interaction between individual rule-based decision-making per crop and grid cell. To compute the required multi-year averages of monthly temperature and precipitation, the data is stored over a 10-year rolling window in the world's and cell's output attribute, updated annually via pycoupler, and continuously averaged during the simulation. Together with a vectorised global xarray-based implementation of the rule set at World level this approach is comparatively concise and computationally efficient compared to its reference . Figure b shows the corresponding processes of the Crop Calendar rule-based model logic.

3.2.3Example simulation results

To verify the suitability of copan:LPJmL as a framework for Crop Calendar, it was implemented in its entirety and applied with one new climate input data set (SSP460, climate model: IPSL-CM6A-LR) from ISIMIP3b to reproduce the original approach from . Similar to this study, we conducted a comparison of recent and future (2080–2099) sowing and harvest dates for two important crop types, maize and temperate cereals (Fig. ). For maize, the implementation dynamically adapts sowing dates in response to climatic changes, enabling earlier sowing in large areas of the temperate latitudes (Fig. a), while in the tropical regions, patterns are more heterogeneous, with both earlier and later sowing, depending on local conditions. Maturity dates (Fig. b) show greater spatial variability, often, but not always, following the pattern of earlier sowing. In some cases (e.g., India and China), maturity is reached later despite earlier sowing.

For temperate cereals, two varieties – a spring and a winter variety – are distinguished, which are sown in different seasons. The model not only adapts sowing dates for each variety but also allows for variety switching based on climate thresholds. While major variety switches remain rare due to modest warming, the model captures spatially nuanced shifts in sowing timing, with earlier sowing in parts of Canada, Europe, Russia, China, and India, and later sowing in eastern Europe and the USA. Maturity dates shift accordingly, with earlier maturity in the USA and southern Europe, and later dates, especially in higher latitudes.

These changes emerge from the dynamic rule-based responses to the climate input as indicated in Fig. a and b, showing the change in average annual temperature and annual precipitation between the recent and future time steps. For example, increases in spring temperatures extend the viable growing season and enable earlier sowing in temperature-sensitive regions. The calculation of the maturity date involves a multi-step approach ( for more details): (i) climate-sensitive harvest rule classification and (ii) harvest date and reason determination based on thresholds such as wet season or the warmest days. Figure c–f shows that harvest rules for maize shift the harvest reason in response to climate change, for instance, a transition in harvest reason from “mid temperature/precipitation” to “high temperature/precipitation” in South America, and from “mid temperature/mixed” to “high temperature/mixed” in parts of the USA and China (Fig. c and d). Although this is only a comparison between two time steps, the underlying runtime algorithm provides insight into each simulation year.

The results demonstrate that the copan:LPJmL-based crop calendar implementation successfully reproduces the climate-responsive adaptation of sowing and maturity dates of . However, unlike the reference implementation of , which requires additional preprocessing steps of the corresponding data products, the copan:LPJmL-based implementation enables a direct, runtime calculation of sowing and harvest dates under changing conditions. This makes the system more suitable for ensemble climate simulations or further SES applications in models such as InSEEDS, where crop production needs to coevolve endogenously to changing biophysical conditions.

Finally, this approach can also be applied to other existing rule-based model systems, such as livestock densities as a function of past grassland performance or whether to plant cover crops in the off-season .

Figure 7

Recreated figure after (Fig. 1) showing the difference in simulated sowing (a) and maturing (b) between no adaptation and timely adaptation for a climate period 2080–2099 in an SSP460 scenario using ISIMIP3b data .

3.3Enhancing classical modelling with LLMs: LLM fertilization 3.3.1Modelling approach

In recent years, there has been an increasing number of modelling approaches using LLMs to emulate agent behaviour in various use cases, spanning from modelling mobility choices to agents' behaviour in online forums . This strong synergy with ABM stems from the fact that LLMs are inherently trained to model human language, reasoning, and decision-making patterns . As such, they are well-suited to serve as proxies for heterogeneous agents, whether individuals, households, or institutions – by generating context-sensitive decisions, goals, or narratives based on inputs from their environment. This capacity makes LLMs particularly compatible with the core idea of ABMs: Simulating the interactions and adaptive behaviour of autonomous entities in a shared environment. However, the integration of LLMs is not limited to ABMs. Their ability to translate between qualitative knowledge and formal rules makes them suitable for enhancing RBMs, for example, by extracting management logic or institutional rules from text sources. Yet, it is in agent-based environments where the conversational and decision-oriented nature of LLMs most directly reflects the modelled processes, making ABMs the current frontier for LLM integration. In the field of land use and agricultural management system, this approach has been taken up by to simulate institutional agency of land use dynamics . Using the copan:LPJmL framework, such an approach can also be used to represent the underlying coevolutionary social-ecological dynamics.

Figure 8

Detailed scheme of the LLM-based fertilization approach highlighted in Fig. representing the feedback mechanisms between the taxons ENV and MET. The LLM gets a prompt with the context of the application, information about location and the cultivated crops, as well as ENV output on crop-specific fertilizer applications and corresponding crop yields. This information helps the LLM to decide whether to increase, maintain, or decrease fertilizer levels for the next simulation year.

3.3.2Model design

As proof of concept for being able to enhance WEMs with LLM agents, we have developed a simple model for nitrogen fertilizer application based on copan:LPJmL , whose approach is illustrated in Fig. . The basis is similar to that of InSEEDS (Sect. ), in which one farming agent is initialized per grid cell to make decisions about local farming practices based on LPJmL (ENV) observations. However, except for the interface between observations and fertilizer application located in the MET taxon, the farmer or decision-making process has been completely outsourced to the LLM. This way CUL processes are not explicitly represented, neither through modelled interactions between farmers, nor through predefined decision rules based on social norms or beliefs. For simplicity, we have therefore decided to omit the CUL taxon in Fig. , even though CUL processes might be part of the LLM reasoning. At the beginning of each simulation year, the farming agent is given a prompt like in Fig. (full prompt available at https://github.com/jnnsbrr/landmanager, last access: 20 May 2026) including ENV observations of the farmer and the request to make decisions based on this knowledge. In the fertilization example described here, the initialized LLM-based farming agents know their geographical position, the crops (CFTs) they grow, the share of land they cultivate, the amount of nitrogen fertilizer applied in each of the last 10 years, and the resulting crop yield for the corresponding years. Based on this knowledge, the LLM-farmers decide on the amount of nitrogen fertilizer they want to apply to each of their crops in the next year, to “increase the crop yield by increasing the application of nitrogen fertilizer as long as it is reasonable” (full prompt in Fig. S2). Moreover, the LLM-farmers are asked to provide their reasoning for the decision taken.

3.3.3Example simulation results

Results of this model for Togo are depicted in Fig. . The map (a) shows that the Nitrogen fertilization level varies strongly from cell to cell and thus proves that the implemented LLM-farmers adopt the Nitrogen fertilization level differently depending on their local needs. After the coupling in 2024, the Nitrogen fertilization levels (b) of all three displayed crops are strongly increased by the LLM-farmer with respective higher harvests (c), with the LLM-farmer reasoning “Increased nitrogen for rainfed rice, maize, and tropical cereals due to low historical application and potential yield gains” (full prompt in Fig. S2). After that, the LLM-agents are only making smaller adjustments to the N fertilization levels, reacting to declines in the harvest. While the maize harvest seems to be stable over the simulated time, rice and tropical cereals show more fluctuating harvests. Thus, the LLM Agent in this cell holds the Nitrogen fertilization level stable for maize, while it adjusts its application to the other two crops. This shows that the LLM-farmer can adapt to changes in the conditions without taking unreasonable decisions.

While this only constitutes a demonstration case, it illustrates the potential of integrating LLM-based agents with copan:LPJmL to introduce more responsive and context-sensitive management decisions. This may offer a more flexible alternative to static assumptions, such as fixed fertilizer levels – in traditional model configurations.

Figure 9

Results of the LLM-based farmer model for Togo in 2050. On the left (a), accumulated nitrogen fertilizer levels between 2024 and 2050 compared to an offline LPJmL run are displayed. On the right, the Nitrogen fertilization (b) and harvest of three crops for one cell ((c), marked as red in (a)) are displayed over time.

4Discussion

A central aim of copan:LPJmL is to provide a platform for connecting modelling paradigms, research domains, and communities (Sect. ). In particular, we seek to enable the integration of process-rich, dynamic Earth system model components with models of social dynamics, supporting use cases ranging from SES science to IAM applications. Previous studies have highlighted key limitations in many agent-based and optimization models, which often lack either the integration of detailed Earth system dynamics or the representation of fundamental social systems and their underlying processes . Conversely, ESMs and DGVMs frequently rely on static or oversimplified representations of human behaviour and decision-making (Sect. ), despite their otherwise detailed process-based structures. The copan:LPJmL framework addresses these gaps by providing a flexible and extensible WEM framework that supports the realisation of hybrid approaches, unifying natural and social systems science in a synergistic way (Sect. ). It allows, for instance, the coupling of agent-based decision-making, rule-based logic, or complex dynamical systems with the biophysical processes of LPJmL. The integrated architecture enables the exploration of social-ecological feedbacks and coevolutionary mechanisms across multiple spatial and temporal scales. The current implementations operate based on the LPJmL grid at 0.5° spatial resolution, where grid cells represent aggregated socio-ecological systems rather than individual land holdings. Depending on the research question, this setup can support both aggregated and more heterogeneous representations of social, economic, and institutional processes within and across cells. In particular, copan:LPJmL is well-suited for research questions that require explicit representation of feedbacks between environmental change and human decision-making, such as agricultural adaptation, land-use transitions, policy interventions, or the dispersal rates of management practices. Beyond the specific implementations presented here, the framework enables the representation of a wide range of social and decision-making processes. At the individual level, behavioural decision models can be implemented using agent-based approaches that draw on different theoretical perspectives, ranging from utility-based or economically rational decision-making to bounded rationality and psychologically grounded frameworks such as the Theory of Planned Behaviour, Value–Belief–Norm theory, or social identity approaches . At higher organisational levels, entities such as groups or social systems can represent collective actors including governments, institutions, markets, or supply chains, enabling the exploration of governance dynamics, policy interventions, and economic interactions. The examples presented in this paper illustrate three different paradigms of decision representation within this framework: behavioural agent-based decision-making by InSEEDS (Sect. ; ), rule-based management logic using Crop Calendar (Sect. ), and LLM-based decision emulation (Sect. ). Beyond the examples presented here, the framework supports a wide range of additional modelling approaches and combinations thereof. For example, combinations are possible in which InSEEDS farmer agents (example 1) uses (parts of) the Crop Calendar (example 2). The possibilities within the MET and CUL taxon are manifold and include individual decision-making, social learning, institutional rule systems, market interactions, and policy interventions, depending on the modelling paradigm implemented . The copan:LPJmL framework builds upon copan:CORE and its modular and open structure that has now been enhanced by the LPJmL ENV integration – all without additionally coupling LPJmL within the model code itself. This lowers the entry barrier and fosters interoperability for modellers from diverse fields enabling researchers to incorporate customised decision rules, agent interactions, additional entities such as governance agents or cooperations, and feedback mechanisms, without requiring modification of the LPJmL model. This is also the main advantage compared to many pre-existing coupler libraries e.g.,. Both the LPJmL coupler library and the Python interface pycoupler (Sect. ) were purpose-built to couple LPJmL with the copan:CORE framework, which is the basis of the integrated copan:LPJmL framework and supports the implementation of numerous projects and ideas, some of which are illustrated in Sect. . This approach eliminates the need to integrate libraries that may require advanced technical know-how – such as BMI – for model coupling. At the same time, copan:LPJmL follows the FAIR principles for research software in providing a findable, accessible, interoperable, and reusable modelling software , with extensive documentation and tutorials available at https://copanlpjml.pik-potsdam.de (last access: 20 May 2026). As with all integrated modelling frameworks, the use of copan:LPJmL involves modelling trade-offs that reflect both computational constraints and the complexity of coupled human–Earth systems. Currently, the coupling between ENV and MET/CUL is only possible on an annual basis, constrained by the LPJmL coupler library and pycoupler (Sect. ). As with all dynamically coupled modelling frameworks, the development of any of the coupled components requires testing and vetting of the coupled system and eventually co-development. Simulated decision-making in MET/CUL does not only depend on plausible decision-making mechanisms, but also on plausible and quantitatively accurate simulated responses in ENV and vice versa. The interpretation of simulation results therefore requires considering both, the behavioural assumptions of the social model and the biophysical responses of the Earth system component. This complexity often results in the continued use of well-tested legacy model versions in coupled systems e.g., . Required testing could be facilitated by a model validation tool chain to allow the harmonization and integration of data sets from various sources such as FAOSTAT for comparing these reference data with simulated data from the coupled system and stand-alone components (e.g. LPJmL) and at a later stage for calibrating the parametrizations of model components. As for now, the functionality of copan:LPJmL and each copan:LPJmL-based model presented here (Sect. ) is backed up by unit and integrity tests to verify their functional and internal validity. With copan:LPJmL, we intend to provide a modelling framework to address research questions around the complex dynamics of the Anthropocene and its coevolution of human and natural systems. While we have demonstrated several types of models based on this framework (Sect. ), we expect that more features, revisions and extensions will be necessary for future models built on this framework. The open-source basis of all model components should facilitate necessary changes . Such extensions can build on the wide range of coupling interfaces and model variables provided by LPJmL and copan:CORE, including the LPJmL input and output categories summarised in Tables and , as well as the generic entity types and process taxa of copan:CORE. This enables the implementation of additional entities, processes, and interaction structures within coupled social–ecological systems. At the same time, future developments remain constrained by the process representations and temporal coupling structure of LPJmL, computational limitations, and the need to balance model complexity with tractability and interpretability.

5Summary and Outlook

This paper introduces copan:LPJmL, a new modelling framework designed to build World-Earth models with a process-rich Dynamic Global Vegetation Model (LPJmL) and a flexible and modular core providing the structures and functionalities to represent various kinds of socio-cultural and socio-metabolic structures and processes. Models built on this framework can represent social and biophysical dynamics in a consistent and coevolutionary manner. With minimal overhead and a lightweight Python interface, copan:LPJmL enables the coupling of LPJmL to diverse types of decision logics – from top-down rule-based models to agent-based bottom-up dynamics – without the need for modifying the LPJmL model itself. By integrating LPJmL as the single component of the ENV taxon into the copan:CORE framework and hierarchy, the framework supports easy access to LPJmL outputs and enables dynamic adjustments to inputs via flexible coupling mechanisms. The three examples presented and discussed in this paper illustrate the breadth of modelling approaches that copan:LPJmL can accommodate: from the top-down, vectorized Crop Calendar, to the bottom-up, agent-based InSEEDS, and towards novel, experimental LLM-based setups. These diverse implementations underscore the framework’s adaptability to a wide range of research questions related to land-use dynamics, adaptation, and the coevolution of human and natural systems. In particular, copan:LPJmL enables the exploration of social–ecological feedbacks in which decision-making processes interact with biophysical dynamics represented by LPJmL. Future work could build on this setup to address gaps in agricultural management data for future projections, support research on planetary boundaries, or explore endogenous policy interventions, including LLM-driven decision processes. Copan:LPJmL provides a flexible foundation for these and many related applications. By adhering to FAIR principles and providing extensive documentation, copan:LPJmL invites collaboration across research domains and communities. It lowers technical barriers for incorporating social dynamics into process-based Earth system modelling and creates a space for new perspectives on topics such as food security, land-use resilience, and regenerative transformations. Ultimately, the framework aims to support the growing need for integrated tools that enable better understanding – and shaping – of coupled human–Earth dynamics in the Anthropocene.

Appendix AConfiguration of LPJmLListing A1

Excerpt of the LPJmL base configuration file lpjml_config.cjson with coupled model name (coupled_model), host (coupled_host), and port (coupled_port), start year of coupling (start_coupling) as well as input and outputs to be coupled.

Listing A2

Configuration of a coupled simulation with LPJmL via the LPJmLConfig module of pycoupler. A base configuration is read in as a LPJmLCoupler object and changed for a coupled simulation, including a regridding to simulate the Netherlands only.

Listing A3

Execution of a coupled simulation with LPJmL via the LPJmLCoupler module and run function of pycoupler that is based on the coupler extension (Sect. ). inputs and outputs are objects of class LPJmLData and can be accessed and edited following .

Appendix BAdditional figuresFigure B1

Additional crop calendar variables: (a) Global temperature and (b) precipitation change from 2025 to 2090 in the SSP460 scenario using ISIMIP3b . Harvest rules (c, d, e, f) following for both CFTs, temperate cereals and maize for 2025 and 2090. Harvest reason (g, h, i, j) also follows with a similar order.

Code and data availability

The copan:LPJmL framework is composed of four software components, each maintained in its own GitHub repository, and the version applied in this paper is archived on Zenodo. LPJmL is available at https://github.com/pik-lpjml/LPJmL (last access: 20 May 2026) and 10.5281/zenodo.17036653 , pycoupler at https://github.com/pik-lpjml/pycoupler (last access: 20 May 2026) and 10.5281/zenodo.17192633 , pycopancore at https://github.com/pik-copan/pycopancore (last access: 20 May 2026) and 10.5281/zenodo.14938316 , and pycopanlpjml at https://github.com/pik-copan/pycopanlpjml (last access: 20 May 2026) . All components are licensed under the GNU General Public License v3.0, except pycopancore, which is distributed under the BSD 2-Clause License. The models InSEEDS as well as Crop Calendar and LLM Fertilization (both part of the landmanager library) have been developed under the same license and are available at https://github.com/pik-copan/inseeds (last access: 20 May 2026) and 10.5281/zenodo.17037107 and https://github.com/jnnsbrr/landmanager (last access: 20 May 2026) and 10.5281/zenodo.17037004 . Comprehensive documentation of copan:LPJmL, including installation instructions, tutorials, a complete API overview, and usage examples, is available at https://copanlpjml.pik-potsdam.de (last access: 20 May 2026) .

The historical climate data set (GSWP-W5E5) that has been used for the InSEEDS and LLM fertilizer simulations as well as the future scenario data set, used for the Crop Calendar simulations (SSP460, IPSL-CM6A-LR) are both available on the ISIMIP homepage https://data.isimip.org (last access: 20 May 2026) , 10.48364/ISIMIP.581124.5 . All further data is either linked directly to the model that is archived together with the model code, model outputs, and scripts that have been used to produce the results presented in this paper on Zenodo 10.5281/zenodo.17054847, .

Author contributions

JB conceived and designed the study. JFD, DG, and CM contributed to the conceptual design of the framework and supervised the work. JB, WvB, and HP led the development of the copan:LPJmL software as well as the models and application examples. SBW, MB, and LS contributed to model and software development. JB, SBW, and HP performed model simulations, analysed data, and visualised results. JB and LS led the writing of the manuscript. CM, SBW, HP, DG, WvB, and JFD contributed to writing and manuscript revision. All authors discussed the simulation results and reviewed and edited the final manuscript.

Competing interests

At least one of the (co-)authors is a member of the editorial board of Geoscientific Model Development. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.

Disclaimer

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.

Acknowledgements

The authors thank the European Regional Development Fund (ERDF), the German Federal Ministry of Education and Research (BMBF), and the State of Brandenburg for providing access to the high-performance computing infrastructure at the Potsdam Institute for Climate Impact Research. The authors used LLM tools, including ChatGPT, to improve the readability of the manuscript. All content was subsequently reviewed, edited, and approved by the authors, who take full responsibility for the published work.

Financial support

J.B., L.S., H.P., D.G., and J.F.D. were supported by the Generation Foundation, the Global Challenges Foundation, and Partners for a New Economy through the Earth4All project. They further acknowledge funding from the European Union's Horizon Europe programme under grant agreement No. 101081661 (WorldTrans). J.F.D. and J.B. were supported by the Carl-Zeiss-Stiftung through the project CZS Research Groups for Earth System Modelling. J.F.D. acknowledges additional funding from the European Research Council (ERC Advanced Grant ERA – Earth Resilience in the Anthropocene, grant no. ERC-2016-ADG-743080) and from the German Federal Ministry of Education and Research (BMBF) through the project CHANGES (grant no. 01LS2001A).The article processing charges for this open-access publication were covered by the Potsdam Institute for Climate Impact Research (PIK).

Review statement

This paper was edited by Roslyn Henry and reviewed by Calum Brown and one anonymous referee.

References Adriaenssens et al.(2004)Adriaenssens, Baets, Goethals, and Pauw

Adriaenssens, V., Baets, B. D., Goethals, P. L. M., and Pauw, N. D.: Fuzzy rule-based models for decision support in ecosystem management, Sci. Total Environ., 319, 1–12, 10.1016/S0048-9697(03)00433-9, 2004.

Ajzen(1985)

Ajzen, I.: From Intentions to Actions: A Theory of Planned Behavior, in: Action Control: From Cognition to Behavior, edited by: Kuhl, J. and Beckmann, J., Springer, Berlin, Heidelberg, 11–39, ISBN 978-3-642-69746-3, 10.1007/978-3-642-69746-3_2, 1985.

Anderies et al.(2023)Anderies, Barfuss, Donges, Fetzer, Heitzig, and Rockström

Anderies, J. M., Barfuss, W., Donges, J. F., Fetzer, I., Heitzig, J., and Rockström, J.: A modeling framework for World-Earth system resilience: exploring social inequality and Earth system tipping points, Environ. Res. Lett., 18, 095001, 10.1088/1748-9326/ace91d, 2023.

Arneth et al.(2014)Arneth, Brown, and Rounsevell

Arneth, A., Brown, C., and Rounsevell, M. D. A.: Global models of human decision-making for land-based mitigation and adaptation assessment, Nat. Clim. Change, 4, 550–557, 10.1038/nclimate2250, 2014.

Arnold et al.(2018)Arnold, Bieger, White, Srinivasan, Dunbar, and Allen

Arnold, J. G., Bieger, K., White, M. J., Srinivasan, R., Dunbar, J. A., and Allen, P. M.: Use of Decision Tables to Simulate Management in SWAT+, Water, 10, 713, 10.3390/w10060713, 2018.

Asefi-Najafabady et al.(2021)Asefi-Najafabady, , , and Morgan

Asefi-Najafabady, S., Laura, V.-O., and Morgan, J.: The failure of Integrated Assessment Models as a response to “climate emergency” and ecological breakdown: the Emperor has no clothes, Globalizations, 18, 1178–1188, 10.1080/14747731.2020.1853958, 2021.

Barker et al.(2022)Barker, Chue Hong, Katz, Lamprecht, Martinez-Ortiz, Psomopoulos, Harrow, Castro, Gruenpeter, Martinez, and Honeyman

Barker, M., Chue Hong, N. P., Katz, D. S., Lamprecht, A.-L., Martinez-Ortiz, C., Psomopoulos, F., Harrow, J., Castro, L. J., Gruenpeter, M., Martinez, P. A., and Honeyman, T.: Introducing the FAIR Principles for research software, Sci. Data, 9, 622, 10.1038/s41597-022-01710-x, 2022.

Bauch et al.(2016)Bauch, Sigdel, Pharaon, and Anand

Bauch, C. T., Sigdel, R., Pharaon, J., and Anand, M.: Early warning signals of regime shifts in coupled human–environment systems, P. Natl. Acad. Sci. USA, 113, 14560–14567, 10.1073/pnas.1604978113, 2016.

Bauer et al.(2020)Bauer, Klein, Humpenöder, Kriegler, Luderer, Popp, and Strefler

Bauer, N., Klein, D., Humpenöder, F., Kriegler, E., Luderer, G., Popp, A., and Strefler, J.: Bio-energy and CO2 emission reductions: an integrated land-use and energy sector perspective, Climatic Change, 163, 1675–1693, 10.1007/s10584-020-02895-z, 2020.

Bechthold et al.(2024)Bechthold, Heitzig, Donges, Barfuss, Marc, Zimmerer, Kolb, Kolster, Müller-Hansen, and Breitbach

Bechthold, M., Heitzig, J., Donges, J., Barfuss, W., Marc, W., Zimmerer, K. B., Kolb, J., Kolster, T., Müller-Hansen, F., and Breitbach, P.: pycopancore: Reference implementation of the copan:CORE World-Earth modelling framework, Zenodo [code], 10.5281/zenodo.12654089, 2024.

Bechthold et al.(2025)Bechthold, Barfuss, Butz, Breier, Constantino, Heitzig, Schwarz, Vardag, and Donges

Bechthold, M., Barfuss, W., Butz, A., Breier, J., Constantino, S. M., Heitzig, J., Schwarz, L., Vardag, S. N., and Donges, J. F.: Social norms and groups structure safe operating spaces in renewable resource use in a social–ecological multi-layer network model, Earth Syst. Dynam., 16, 1365–1390, 10.5194/esd-16-1365-2025, 2025.

Beck and Krueger(2016)

Beck, M. and Krueger, T.: The epistemic, ethical, and political dimensions of uncertainty in integrated assessment modeling, WIREs Climate Change, 7, 627–645, 10.1002/wcc.415, _eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/wcc.415, 2016.

Beckage et al.(2020)Beckage, Lacasse, Winter, Gross, Fefferman, Hoffman, Metcalf, Franck, Carr, Zia, and Kinzig

Beckage, B., Lacasse, K., Winter, J. M., Gross, L. J., Fefferman, N., Hoffman, F. M., Metcalf, S. S., Franck, T., Carr, E., Zia, A., and Kinzig, A.: The Earth has humans, so why don’t our climate models?, Climatic Change, 163, 181–188, 10.1007/s10584-020-02897-x, 2020.

Beckage et al.(2022)Beckage, Moore, and Lacasse

Beckage, B., Moore, F. C., and Lacasse, K.: Incorporating human behaviour into Earth system modelling, Nature Human Behaviour, 6, 1493–1502, 10.1038/s41562-022-01478-5, 2022.

Beier et al.(2025)Beier, Dietrich, Heinke, Abrahao, Jeetze, Bodirsky, Crawford, Humpenöder, Merfort, Weindl, Herrero, Mason-D’Croz, Rockström, Sundiang, Wierik, Norberg, Klein, Müller, Lotze-Campen, and Popp

Beier, F. D., Dietrich, J. P., Heinke, J., Abrahao, G., Jeetze, P. v., Bodirsky, B. L., Crawford, M., Humpenöder, F., Merfort, L., Weindl, I., Herrero, M., Mason-D’Croz, D., Rockström, J., Sundiang, M., Wierik, S. t., Norberg, A., Klein, D., Müller, C., Lotze-Campen, H., and Popp, A.: Planetary boundaries under a land-based climate change mitigation scenario with a food demand transformation: a modelling study, Lancet Planetary Health, 10.1016/S2542-5196(25)00087-7, 2025.

Blanco-Gutiérrez et al.(2013)Blanco-Gutiérrez, Varela-Ortega, and Purkey

Blanco-Gutiérrez, I., Varela-Ortega, C., and Purkey, D. R.: Integrated assessment of policy interventions for promoting sustainable irrigation in semi-arid environments: A hydro-economic modeling approach, J. Environ. Manage., 128, 144–160, 10.1016/j.jenvman.2013.04.037, 2013.

Bodirsky et al.(2012)Bodirsky, Popp, Weindl, Dietrich, Rolinski, Scheiffele, Schmitz, and Lotze-Campen

Bodirsky, B. L., Popp, A., Weindl, I., Dietrich, J. P., Rolinski, S., Scheiffele, L., Schmitz, C., and Lotze-Campen, H.: N2O emissions from the global agricultural nitrogen cycle – current state and future scenarios, Biogeosciences, 9, 4169–4197, 10.5194/bg-9-4169-2012, 2012.

Bondeau et al.(2007)Bondeau, Smith, Zaehle, Schaphoff, Lucht, Cramer, Gerten, Lotze-Campen, Müller, Reichstein, and Smith

Bondeau, A., Smith, P. C., Zaehle, S., Schaphoff, S., Lucht, W., Cramer, W., Gerten, D., Lotze-Campen, H., Müller, C., Reichstein, M., and Smith, B.: Modelling the role of agriculture for the 20th century global terrestrial carbon balance, Glob. Change Biol., 13, 679–706, 10.1111/j.1365-2486.2006.01305.x, 2007.

Braakhekke et al.(2019)Braakhekke, Doelman, Baas, Müller, Schaphoff, Stehfest, and van Vuuren

Braakhekke, M. C., Doelman, J. C., Baas, P., Müller, C., Schaphoff, S., Stehfest, E., and van Vuuren, D. P.: Modeling forest plantations for carbon uptake with the LPJmL dynamic global vegetation model, Earth Syst. Dynam., 10, 617–630, 10.5194/esd-10-617-2019, 2019.

Breier(2025)

Breier, J.: copan:LPJmL homepage, https://copanlpjml.pik-potsdam.de (last access: 20 May 2026), 2025.

Breier and Prawitz(2025)

Breier, J. and Prawitz, H.: Data for GMD submission “copan:LPJmL: A new hybrid modelling framework for dynamic land use and agricultural management, Zenodo [data set], 10.5281/zenodo.17054847, 2025.

Breier and von Bloh(2025)

Breier, J. and von Bloh, W.: pycoupler: dynamic model coupling of LPJmL, Zenodo [code], 10.5281/zenodo.17192633, 2025.

Breier et al.(2024)Breier, Ostberg, Wirth, Hötten, Minoli, Stenzel, and Müller

Breier, J., Ostberg, S., Wirth, S. B., Hötten, D., Minoli, S., Stenzel, F., and Müller, C.: lpjmlkit: A toolkit for operating LPJmL and model-specific data processing, J. Open Source Softw., 9, 5447, 10.21105/joss.05447, 2024.

Breier et al.(2025a)Breier, Prawitz, and Wirth

Breier, J., Prawitz, H., and Wirth, S. B.: Example models of land management based on copan:LPJmL, Zenodo [code], 10.5281/zenodo.17037004, 2025a.

Breier et al.(2025b)Breier, Schwarz, Hotz, Bechthold, Heitzig, Donges, Barfuss, Marc, Zimmerer, Kolb, Kolster, Müller-Hansen, and Breitbach

Breier, J., Schwarz, L., Hotz, R., Bechthold, M., Heitzig, J., Donges, J., Barfuss, W., Marc, W., Zimmerer, K. B., Kolb, J., Kolster, T., Müller-Hansen, F., and Breitbach, P.: pycopancore: Reference implementation of the copan:CORE World-Earth modelling framework, Zenodo [code], 10.5281/zenodo.14938316, 2025b.

Breier et al.(2025c)Breier, Schwarz, Prawitz, Bechthold, and Hotz

Breier, J., Schwarz, L., Prawitz, H., Bechthold, M., and Hotz, R.: Model of integrated social-ecological resilient land systems (InSEEDS), Zenodo [code], 10.5281/zenodo.17037107, 2025c.

Breier et al.(2025d)

Breier, J., Schwarz, L., Prawitz, H., and Bechthold, M.: copan:LPJmL, an advanced World-Earth modeling framework extending copan:CORE, integrating LPJmL as the Earth system interface for comprehensive social-ecological simulations, Zenodo [code], 10.5281/zenodo.17036693, 2025d.

Bury et al.(2019)Bury, Bauch, and Anand

Bury, T. M., Bauch, C. T., and Anand, M.: Charting pathways to climate change mitigation in a coupled socio-climate model, PLOS Comput. Biol., 15, e1007000, 10.1371/journal.pcbi.1007000, 2019.

Calvin and Bond-Lamberty(2018)

Calvin, K. and Bond-Lamberty, B.: Integrated human-earth system modeling-state of the science and future directions, Environ. Res. Lett., 13, 063006, 10.1088/1748-9326/aac642, 2018.

Campbell et al.(2017)Campbell, Beare, Bennett, Hall-Spencer, Ingram, Jaramillo, Ortiz, Ramankutty, Sayer, and Shindell

Campbell, B. M., Beare, D. J., Bennett, E. M., Hall-Spencer, J. M., Ingram, J. S. I., Jaramillo, F., Ortiz, R., Ramankutty, N., Sayer, J. A., and Shindell, D.: Agriculture production as a major driver of the Earth system exceeding planetary boundaries, Ecology and Society, 22, https://www.jstor.org/stable/26798991 (last access: 20 May 2026), 2017.

Chen and Huang(2024)

Chen, D. and Huang, Y.: Integrating Reinforcement Learning and Large Language Models for Crop Production Process Management Optimization and Control through A New Knowledge-Based Deep Learning Paradigm, arXiv [preprint], 10.48550/arXiv.2410.09680, 2024.

Crutzen(2002)

Crutzen, P. J.: The “anthropocene”, Journal de Physique IV (Proceedings), 12, 1–5, 10.1051/jp4:20020447, 2002.

Dantas de Paula et al.(2020)Dantas de Paula, Gómez Giménez, Niamir, Thurner, and Hickler

Dantas de Paula, M., Gómez Giménez, M., Niamir, A., Thurner, M., and Hickler, T.: Combining European Earth Observation products with Dynamic Global Vegetation Models for estimating Essential Biodiversity Variables, Int. J. Digit. Earth, 13, 262–277, 10.1080/17538947.2019.1597187, 2020.

Dietrich et al.(2019)Dietrich, Bodirsky, Humpenöder, Weindl, Stevanović, Karstens, Kreidenweis, Wang, Mishra, Klein, Ambrósio, Araujo, Yalew, Baumstark, Wirth, Giannousakis, Beier, Chen, Lotze-Campen, and Popp

Dietrich, J. P., Bodirsky, B. L., Humpenöder, F., Weindl, I., Stevanović, M., Karstens, K., Kreidenweis, U., Wang, X., Mishra, A., Klein, D., Ambrósio, G., Araujo, E., Yalew, A. W., Baumstark, L., Wirth, S., Giannousakis, A., Beier, F., Chen, D. M.-C., Lotze-Campen, H., and Popp, A.: MAgPIE 4 – a modular open-source framework for modeling global land systems, Geosci. Model Dev., 12, 1299–1317, 10.5194/gmd-12-1299-2019, 2019.

Doelman et al.(2020)Doelman, Stehfest, van Vuuren, Tabeau, Hof, Braakhekke, Gernaat, van den Berg, van Zeist, Daioglou, van Meijl, and Lucas

Doelman, J. C., Stehfest, E., van Vuuren, D. P., Tabeau, A., Hof, A. F., Braakhekke, M. C., Gernaat, D. E. H. J., van den Berg, M., van Zeist, W.-J., Daioglou, V., van Meijl, H., and Lucas, P. L.: Afforestation for climate change mitigation: Potentials, risks and trade-offs, Glob. Change Biol., 26, 1576–1591, 10.1111/gcb.14887, 2020.

Donges et al.(2017)Donges, Winkelmann, Lucht, Cornell, Dyke, Rockström, Heitzig, and Schellnhuber

Donges, J. F., Winkelmann, R., Lucht, W., Cornell, S. E., Dyke, J. G., Rockström, J., Heitzig, J., and Schellnhuber, H. J.: Closing the loop: Reconnecting human dynamics to Earth System science, Anthropocene Rev., 4, 151–157, 10.1177/2053019617725537, 2017.

Donges et al.(2020)Donges, Heitzig, Barfuss, Wiedermann, Kassel, Kittel, Kolb, Kolster, Müller-Hansen, Otto, Zimmerer, and Lucht

Donges, J. F., Heitzig, J., Barfuss, W., Wiedermann, M., Kassel, J. A., Kittel, T., Kolb, J. J., Kolster, T., Müller-Hansen, F., Otto, I. M., Zimmerer, K. B., and Lucht, W.: Earth system modeling with endogenous and dynamic human societies: the copan:CORE open World–Earth modeling framework, Earth Syst. Dynam., 11, 395–413, 10.5194/esd-11-395-2020, 2020.

Donges et al.(2021)Donges, Lucht, Cornell, Heitzig, Barfuss, Lade, and Schlüter

Donges, J. F., Lucht, W., Cornell, S. E., Heitzig, J., Barfuss, W., Lade, S. J., and Schlüter, M.: Taxonomies for structuring models for World–Earth systems analysis of the Anthropocene: subsystems, their interactions and social–ecological feedback loops, Earth Syst. Dynam., 12, 1115–1137, 10.5194/esd-12-1115-2021, 2021.

FAO(2025)

FAO: FAOSTAT, https://www.fao.org/faostat/en/#data (last access: 20 May 2026), 2025.

Farahbakhsh et al.(2022)Farahbakhsh, Bauch, and Anand

Farahbakhsh, I., Bauch, C. T., and Anand, M.: Modelling coupled human–environment complexity for the future of the biosphere: strengths, gaps and promising directions, Philos. T. Roy. Soc. B, 377, 20210382, 10.1098/rstb.2021.0382, 2022.

Filatova et al.(2013)Filatova, Verburg, Parker, and Stannard

Filatova, T., Verburg, P. H., Parker, D. C., and Stannard, C. A.: Spatial agent-based models for socio-ecological systems: Challenges and prospects, Environ. Modell. Softw., 45, 1–7, 10.1016/j.envsoft.2013.03.017, 2013.

Foley et al.(2005)Foley, DeFries, Asner, Barford, Bonan, Carpenter, Chapin, Coe, Daily, Gibbs, Helkowski, Holloway, Howard, Kucharik, Monfreda, Patz, Prentice, Ramankutty, and Snyder

Foley, J. A., DeFries, R., Asner, G. P., Barford, C., Bonan, G., Carpenter, S. R., Chapin, F. S., Coe, M. T., Daily, G. C., Gibbs, H. K., Helkowski, J. H., Holloway, T., Howard, E. A., Kucharik, C. J., Monfreda, C., Patz, J. A., Prentice, I. C., Ramankutty, N., and Snyder, P. K.: Global Consequences of Land Use, Science, 309, 570–574, 10.1126/science.1111772, 2005.

Gao et al.(2024)Gao, Lan, Li, Yuan, Ding, Zhou, Xu, and Li

Gao, C., Lan, X., Li, N., Yuan, Y., Ding, J., Zhou, Z., Xu, F., and Li, Y.: Large language models empowered agent-based modeling and simulation: a survey and perspectives, Humanities and Social Sciences Communications, 11, 1259, 10.1057/s41599-024-03611-3, 2024.

Gerten et al.(2018)Gerten, Schönfeld, and Schauberger

Gerten, D., Schönfeld, M., and Schauberger, B.: On deeper human dimensions in Earth system analysis and modelling, Earth Syst. Dynam., 9, 849–863, 10.5194/esd-9-849-2018, 2018.

Gerten et al.(2020)Gerten, Heck, Jägermeyr, Bodirsky, Fetzer, Jalava, Kummu, Lucht, Rockström, Schaphoff, and Schellnhuber

Gerten, D., Heck, V., Jägermeyr, J., Bodirsky, B. L., Fetzer, I., Jalava, M., Kummu, M., Lucht, W., Rockström, J., Schaphoff, S., and Schellnhuber, H. J.: Feeding ten billion people is possible within four terrestrial planetary boundaries, Nat. Sustain., 3, 200–208, 10.1038/s41893-019-0465-1, 2020.

Hanke et al.(2016)Hanke, Redler, Holfeld, and Yastremsky

Hanke, M., Redler, R., Holfeld, T., and Yastremsky, M.: YAC 1.2.0: new aspects for coupling software in Earth system modelling, Geosci. Model Dev., 9, 2755–2769, 10.5194/gmd-9-2755-2016, 2016.

Heinke et al.(2023)Heinke, Rolinski, and Müller

Heinke, J., Rolinski, S., and Müller, C.: Modelling the role of livestock grazing in C and N cycling in grasslands with LPJmL5.0-grazing, Geosci. Model Dev., 16, 2455–2475, 10.5194/gmd-16-2455-2023, 2023.

Herzfeld et al.(2021)Herzfeld, Heinke, Rolinski, and Müller

Herzfeld, T., Heinke, J., Rolinski, S., and Müller, C.: Soil organic carbon dynamics from agricultural management practices under climate change, Earth Syst. Dynam., 12, 1037–1055, 10.5194/esd-12-1037-2021, 2021.

Holdgate(1987)

Holdgate, M. W.: Our Common Future: The Report of the World Commission on Environment and Development. Oxford University Press, Oxford & New York, xv + 347 + 35 pp., Oxford Paperback, Environmental Conservation, 14, 282–282, 10.1017/S0376892900016702, 1987.

Horsevad et al.(2022)Horsevad, Mateo, Kooij, Barrat, and Bouffanais

Horsevad, N., Mateo, D., Kooij, R. E., Barrat, A., and Bouffanais, R.: Transition from simple to complex contagion in collective decision-making, Nat. Commun., 13, 1442, 10.1038/s41467-022-28958-6, 2022.

Hoyer and Hamman(2017)

Hoyer, S. and Hamman, J.: xarray: N-D labeled Arrays and Datasets in Python, J. Open Res. Softw., 5, 10.5334/jors.148, 2017.

Huber et al.(2022)Huber, Xiong, Keller, and Finger

Huber, R., Xiong, H., Keller, K., and Finger, R.: Bridging behavioural factors and standard bio-economic modelling in an agent-based modelling framework, J. Agr. Econom., 73, 35–63, 10.1111/1477-9552.12447, 2022.

Hutton et al.(2020)Hutton, Piper, and Tucker

Hutton, E. W. h., Piper, M. D., and Tucker, G. E.: The Basic Model Interface 2.0: A standard interface for coupling numerical models in the geosciences, J. Open Source Softw., 5, 2317, 10.21105/joss.02317, 2020.

IPCC(2023)

IPCC: Climate Change 2022: Mitigation of Climate Change. Contribution of Working Group III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Shukla, P. R., Skea, J., Slade, R., Al Khourdajie, A., van Diemen, R., McCollum, D., Pathak, M., Some, S., Vyas, P., Fradera, R., Belkacemi, M., Hasija, A., Lisboa, G., Luz, S., and Malley, J., Cambridge University Press, Cambridge, UK and New York, NY, USA, 10.1017/9781009157926, 2022.

Jones et al.(2003)Jones, Hoogenboom, Porter, Boote, Batchelor, Hunt, Wilkens, Singh, Gijsman, and Ritchie

Jones, J. W., Hoogenboom, G., Porter, C. H., Boote, K. J., Batchelor, W. D., Hunt, L. A., Wilkens, P. W., Singh, U., Gijsman, A. J., and Ritchie, J. T.: The DSSAT cropping system model, Eur. J. Agron., 18, 235–265, 10.1016/S1161-0301(02)00107-7, 2003.

Jägermeyr et al.(2015)Jägermeyr, Gerten, Heinke, Schaphoff, Kummu, and Lucht

Jägermeyr, J., Gerten, D., Heinke, J., Schaphoff, S., Kummu, M., and Lucht, W.: Water savings potentials of irrigation systems: global simulation of processes and linkages, Hydrol. Earth Syst. Sci., 19, 3073–3091, 10.5194/hess-19-3073-2015, 2015.

Jägermeyr et al.(2016)Jägermeyr, Gerten, Schaphoff, Heinke, Lucht, and Rockström

Jägermeyr, J., Gerten, D., Schaphoff, S., Heinke, J., Lucht, W., and Rockström, J.: Integrated crop water management might sustainably halve the global food gap, Environ. Res. Lett., 11, 025002, 10.1088/1748-9326/11/2/025002, 2016.

Koasidis et al.(2023)Koasidis, Nikas, and Doukas

Koasidis, K., Nikas, A., and Doukas, H.: Why integrated assessment models alone are insufficient to navigate us through the polycrisis, One Earth, 6, 205–209, 10.1016/j.oneear.2023.02.009, 2023.

Kou-Giesbrecht et al.(2024)Kou-Giesbrecht, Arora, Seiler, and Wang

Kou-Giesbrecht, S., Arora, V. K., Seiler, C., and Wang, L.: The impacts of modelling prescribed vs. dynamic land cover in a high-CO₂ future scenario – greening of the Arctic and Amazonian dieback, Biogeosciences, 21, 3339–3371, 10.5194/bg-21-3339-2024, 2024.

Krawczyk and Braun(2025)

Krawczyk, F. and Braun, A. C.: Models like heroes? Making Integrated Assessment Models (IAMs) ready for deep decarbonization and a socio-economic transformation, Energy Research & Social Science, 121, 103959, 10.1016/j.erss.2025.103959, 2025.

Kreibich et al.(2025)Kreibich, Sivapalan, AghaKouchak, Addor, Aksoy, Arheimer, Arnbjerg-Nielsen, Vail-Castro, Cudennec, Madruga de Brito, Di Baldassarre, Finger, Fowler, Knoben, Krueger, Liu, Macdonald, McMillan, Mendiondo, Montanari, Muller, Pande, Tian, Viglione, Wei, Castellarin, Loucks, Oki, Polo, Savenije, Van Loon, Agarwal, Alvarez-Garreton, Andreu, Barendrecht, Brunner, Cavalcante, Cavus, Ceola, Chaffe, Chen, Coxon, Dandan, Davary, Dembélé, Dewals, Frolova (Bibikova), Gain, Gelfan, Ghoreishi, Grabs, Guan, Hannah, Helmschrot, Höllermann, Hounkpè, Koebele, Konar, Kratzert, Lindersson, Llasat, Matanó, Mazzoleni, Mejia, Mendoza, Merz, Mukherjee, Nasiri Saleh, Nlend, Nonki, Orieschnig, Papagiannaki, Penny, Petrucci, Pimentel, Pool, Ridolfi, Rusca, Sairam, Sankaran Namboothiri, Sarmento Buarque, Savelli, Schoppa, Schröter, Scolobig, Shafiei, Sikorska-Senoner, Smigaj, Teutschbein, Thaler, Todorovic, Tootoonchi, Tootoonchi, Toth, van Nooijen, Vanelli, Vásquez, Walker, Wens, Yu, Zarei, Zhou, and Blöschl

Kreibich, H., Sivapalan, M., AghaKouchak, A., Addor, N., Aksoy, H., Arheimer, B., Arnbjerg-Nielsen, K., Vail-Castro, C., Cudennec, C., Madruga de Brito, M., Di Baldassarre, G., Finger, D. C., Fowler, K., Knoben, W., Krueger, T., Liu, J., Macdonald, E., McMillan, H., Mendiondo, E. M., Montanari, A., Muller, M. F., Pande, S., Tian, F., Viglione, A., Wei, Y., Castellarin, A., Loucks, D. P., Oki, T., Polo, M. J., Savenije, H., Van Loon, A. F., Agarwal, A., Alvarez-Garreton, C., Andreu, A., Barendrecht, M. H., Brunner, M., Cavalcante, L., Cavus, Y., Ceola, S., Chaffe, P., Chen, X., Coxon, G., Dandan, Z., Davary, K., Dembélé, M., Dewals, B., Frolova (Bibikova), T., Gain, A. K., Gelfan, A., Ghoreishi, M., Grabs, T., Guan, X., Hannah, D. M., Helmschrot, J., Höllermann, B., Hounkpè, J., Koebele, E., Konar, M., Kratzert, F., Lindersson, S., Llasat, M. C., Matanó, A., Mazzoleni, M., Mejia, A., Mendoza, P., Merz, B., Mukherjee, J., Nasiri Saleh, F., Nlend, B., Nonki, R. M., Orieschnig, C., Papagiannaki, K., Penny, G., Petrucci, O., Pimentel, R., Pool, S., Ridolfi, E., Rusca, M., Sairam, N., Sankaran Namboothiri, A., Sarmento Buarque, A. C., Savelli, E., Schoppa, L., Schröter, K., Scolobig, A., Shafiei, M., Sikorska-Senoner, A. E., Smigaj, M., Teutschbein, C., Thaler, T., Todorovic, A., Tootoonchi, F., Tootoonchi, R., Toth, E., van Nooijen, R., Vanelli, F. M., Vásquez, N., Walker, D. W., Wens, M., Yu, D. J., Zarei, H., Zhou, C., and Blöschl, G.: Panta Rhei: a decade of progress in research on change in hydrology and society, Hydrolog. Sci. J., 70, 1210–1236, 10.1080/02626667.2025.2469762, 2025.

Krey et al.(2020)Krey, Havlik, Kishimoto, Fricko, Zilliacus, Gidden, Strubegger, Kartasasmita, Ermolieva, Forsell, Guo, Gusti, Huppmann, Johnson, Kikstra, Kindermann, Kolp, Lovat, McCollum, Min, Pachauri, Parkinson, Rao, Rogelj, Ünlü, Valin, Wagner, Zakeri, Obersteiner, and Riahi

Krey, V., Havlik, P., Kishimoto, P., Fricko, O., Zilliacus, J., Gidden, M., Strubegger, M., Kartasasmita, G., Ermolieva, T., Forsell, N., Guo, F., Gusti, M., Huppmann, D., Johnson, N., Kikstra, J., Kindermann, G., Kolp, P., Lovat, F., McCollum, D., Min, J., Pachauri, S., Parkinson, S., Rao, S., Rogelj, J., Ünlü, G., Valin, H., Wagner, P., Zakeri, B., Obersteiner, M., and Riahi, K.: MESSAGEix-GLOBIOM Documentation – 2020 release, Laxenburg, Austria, 10.22022/IACC/03-2021.17115, 2020.

Lange et al.(2023)Lange, Mengel, Treu, and Büchner

Lange, S., Mengel, M., Treu, S., and Büchner, M.: ISIMIP3a atmospheric climate input data (v1.2), ISIMIP Repository [data set], 10.48364/ISIMIP.982724.2, 2023.

Lange et al.(2024)Lange, Quesada-Chacón, and Büchner

Lange, S., Quesada-Chacón, D., and Büchner, M.: Secondary ISIMIP3b bias-adjusted atmospheric climate input data (v1.5), ISIMIP Repository [data set], 10.48364/ISIMIP.581124.5, 2024.

Lenton(2016)

Lenton, T.: Earth System Science: A Very Short Introduction, Oxford University Press, ISBN 978-0-19-871887-1, google-Books-ID: ijs6CwAAQBAJ, 2016.

Li et al.(2021)Li, Xu, Chen, Gao, Cai, Shen, He, Zhuang, and Chen

Li, R., Xu, M., Chen, Z., Gao, B., Cai, J., Shen, F., He, X., Zhuang, Y., and Chen, D.: Phenology-based classification of crop species and rotation types using fused MODIS and Landsat data: The comparison of a random-forest-based model and a decision-rule-based model, Soil Till. Res., 206, 104838, 10.1016/j.still.2020.104838, 2021.

Li et al.(2016)Li, Yu, Sohl, Clinton, Li, Zhu, Liu, and Gong

Li, X., Yu, L., Sohl, T., Clinton, N., Li, W., Zhu, Z., Liu, X., and Gong, P.: A cellular automata downscaling based 1 km global land use datasets (2010–2100), Sci. Bull., 61, 1651–1661, 10.1007/s11434-016-1148-1, 2016.

Liang et al.(2021)Liang, Guan, Clarke, Liu, Wang, and Yao

Liang, X., Guan, Q., Clarke, K. C., Liu, S., Wang, B., and Yao, Y.: Understanding the drivers of sustainable land expansion using a patch-generating land use simulation (PLUS) model: A case study in Wuhan, China, Computers, Environ. Urban Syst., 85, 101569, 10.1016/j.compenvurbsys.2020.101569, 2021.

Liu et al.(2023)Liu, Ye, Müller, Jägermeyr, Franke, Stephens, and Chen

Liu, W., Ye, T., Müller, C., Jägermeyr, J., Franke, J. A., Stephens, H., and Chen, S.: The statistical emulators of GGCMI phase 2: responses of year-to-year variation of crop yield to CO₂, temperature, water, and nitrogen perturbations, Geosci. Model Dev., 16, 7203–7221, 10.5194/gmd-16-7203-2023, 2023.

Lovelock and Margulis(1974)

Lovelock, J. E. and Margulis, L.: Atmospheric homeostasis by and for the biosphere: the gaia hypothesis, Tellus, 26, 2–10, 10.3402/tellusa.v26i1-2.9731, 1974.

Lutz et al.(2019)Lutz, Herzfeld, Heinke, Rolinski, Schaphoff, von Bloh, Stoorvogel, and Müller

Lutz, F., Herzfeld, T., Heinke, J., Rolinski, S., Schaphoff, S., von Bloh, W., Stoorvogel, J. J., and Müller, C.: Simulating the effect of tillage practices with the global ecosystem model LPJmL (version 5.0-tillage), Geosci. Model Dev., 12, 2419–2440, 10.5194/gmd-12-2419-2019, 2019.

Mathias et al.(2020)Mathias, Debeljak, Deffuant, Diemer, Dierickx, Donges, Gladkykh, Heitzig, Holtz, Obergassel, Pellaud, Sánchez, Trajanov, and Videira

Mathias, J.-D., Debeljak, M., Deffuant, G., Diemer, A., Dierickx, F., Donges, J. F., Gladkykh, G., Heitzig, J., Holtz, G., Obergassel, W., Pellaud, F., Sánchez, A., Trajanov, A., and Videira, N.: Grounding Social Foundations for Integrated Assessment Models of Climate Change, Earth's Future, 8, e2020EF001573, 10.1029/2020EF001573, 2020.

Meyfroidt et al.(2022)Meyfroidt, de Bremond, Ryan, Archer, Aspinall, Chhabra, Camara, Corbera, DeFries, Díaz, Dong, Ellis, Erb, Fisher, Garrett, Golubiewski, Grau, Grove, Haberl, Heinimann, Hostert, Jobbágy, Kerr, Kuemmerle, Lambin, Lavorel, Lele, Mertz, Messerli, Metternicht, Munroe, Nagendra, Nielsen, Ojima, Parker, Pascual, Porter, Ramankutty, Reenberg, Roy Chowdhury, Seto, Seufert, Shibata, Thomson, Turner, Urabe, Veldkamp, Verburg, Zeleke, and zu Ermgassen

Meyfroidt, P., de Bremond, A., Ryan, C. M., Archer, E., Aspinall, R., Chhabra, A., Camara, G., Corbera, E., DeFries, R., Díaz, S., Dong, J., Ellis, E. C., Erb, K.-H., Fisher, J. A., Garrett, R. D., Golubiewski, N. E., Grau, H. R., Grove, J. M., Haberl, H., Heinimann, A., Hostert, P., Jobbágy, E. G., Kerr, S., Kuemmerle, T., Lambin, E. F., Lavorel, S., Lele, S., Mertz, O., Messerli, P., Metternicht, G., Munroe, D. K., Nagendra, H., Nielsen, J. O., Ojima, D. S., Parker, D. C., Pascual, U., Porter, J. R., Ramankutty, N., Reenberg, A., Roy Chowdhury, R., Seto, K. C., Seufert, V., Shibata, H., Thomson, A., Turner, B. L., Urabe, J., Veldkamp, T., Verburg, P. H., Zeleke, G., and zu Ermgassen, E. K. H. J.: Ten facts about land systems for sustainability, P. Natl. Acad. Sci. USA, 119, e2109217118, 10.1073/pnas.2109217118, 2022.

Minoli et al.(2019)Minoli, Egli, Rolinski, and Müller

Minoli, S., Egli, D. B., Rolinski, S., and Müller, C.: Modelling cropping periods of grain crops at the global scale, Global Planet. Change, 174, 35–46, 10.1016/j.gloplacha.2018.12.013, 2019.

Minoli et al.(2022a)Minoli, Jägermeyr, Asseng, Urfels, and Müller

Minoli, S., Jägermeyr, J., Asseng, S., Urfels, A., and Müller, C.: Global crop yields can be lifted by timely adaptation of growing periods to climate change, Nat. Communi., 13, 7079, 10.1038/s41467-022-34411-5, 2022a.

Minoli et al.(2022b)Minoli, Jägermeyr, Asseng, Urfels, and Müller

Minoli, S., Jägermeyr, J., Asseng, S., Urfels, A., and Müller, C.: Workflow for Global crop yields can be lifted by timely adaptation of growing periods to climate change, Zenodo, 10.5281/zenodo.7038163, 2022b.

Moore et al.(2014)Moore, Holzworth, Herrmann, Brown, de Voil, Snow, Zurcher, and Huth

Moore, A. D., Holzworth, D. P., Herrmann, N. I., Brown, H. E., de Voil, P. G., Snow, V. O., Zurcher, E. J., and Huth, N. I.: Modelling the manager: Representing rule-based management in farming systems simulation models, Environ. Modell. Softw., 62, 399–410, 10.1016/j.envsoft.2014.09.001, 2014.

Moore et al.(2022)Moore, Lacasse, Mach, Shin, Gross, and Beckage

Moore, F. C., Lacasse, K., Mach, K. J., Shin, Y. A., Gross, L. J., and Beckage, B.: Determinants of emissions pathways in the coupled climate–social system, Nature, 603, 103–111, 10.1038/s41586-022-04423-8, 2022.

Murray-Rust et al.(2014)Murray-Rust, Brown, van Vliet, Alam, Robinson, Verburg, and Rounsevell

Murray-Rust, D., Brown, C., van Vliet, J., Alam, S. J., Robinson, D. T., Verburg, P. H., and Rounsevell, M.: Combining agent functional types, capitals and services to model land use dynamics, Environ. Modell. Softw., 59, 187–201, 10.1016/j.envsoft.2014.05.019, 2014.

Müller et al.(2016)Müller, Stehfest, Van Minnen, Strengers, Von Bloh, Beusen, Schaphoff, Kram, and Lucht

Müller, C., Stehfest, E., Van Minnen, J. G., Strengers, B., Von Bloh, W., Beusen, A. H. W., Schaphoff, S., Kram, T., and Lucht, W.: Drivers and patterns of land biosphere carbon balance reversal, Environ. Res. Lett., 11, 044002, 10.1088/1748-9326/11/4/044002, 2016.

Müller et al.(2024)Müller, Lange, Fischer, König, Kelbling, Leal Rojas, and Thober

Müller, S., Lange, M., Fischer, T., König, S., Kelbling, M., Leal Rojas, J. J., and Thober, S.: FINAM is not a model (v1.0): a new Python-based model coupling framework, Geosci. Model Dev., 18, 4483–4498, 10.5194/gmd-18-4483-2025, 2025.

Natel et al.(2025)Natel, Belda, Anthoni, Haß, Rabin, and Arneth

Natel, C., Belda, D. M., Anthoni, P., Haß, N., Rabin, S., and Arneth, A.: Emulating grid-based forest carbon dynamics using machine learning: an LPJ-GUESS v4.1.1 application, Geosci. Model Dev., 18, 4317–4333, 10.5194/gmd-18-4317-2025, 2025.

National Research Council(1986)

National Research Council: Earth System Science: Overview – A Program for Global Change, National Academy Press, 10.17226/19210, 1986.

Ngwira et al.(2014)Ngwira, Aune, and Thierfelder

Ngwira, A. R., Aune, J. B., and Thierfelder, C.: DSSAT modelling of conservation agriculture maize response to climate change in Malawi, Soil Till. Res., 143, 85–94, 10.1016/j.still.2014.05.003, 2014.

Porwollik et al.(2021)Porwollik, Rolinski, Heinke, von Bloh, Schaphoff, and Müller

Porwollik, V., Rolinski, S., Heinke, J., von Bloh, W., Schaphoff, S., and Müller, C.: The role of cover crops for cropland soil carbon, nitrogen leaching, and agricultural yields – A global simulation study with LPJmL (V. 5.0-tillage-cc), Biogeosciences Discuss. [preprint], 10.5194/bg-2021-215, 2021.

Porwollik et al.(2022)Porwollik, Rolinski, Heinke, von Bloh, Schaphoff, and Müller

Porwollik, V., Rolinski, S., Heinke, J., von Bloh, W., Schaphoff, S., and Müller, C.: The role of cover crops for cropland soil carbon, nitrogen leaching, and agricultural yields – a global simulation study with LPJmL (V. 5.0-tillage-cc), Biogeosciences, 19, 957–977, 10.5194/bg-19-957-2022, 2022.

Preiser et al.(2018)Preiser, Biggs, De Vos, and Folke

Preiser, R., Biggs, R., De Vos, A., and Folke, C.: Social-ecological systems as complex adaptive systems: organizing principles for advancing research methods and approaches, Ecol. Soc., 23, 10.5751/ES-10558-230446, 2018.

Reick et al.(2013)Reick, Raddatz, Brovkin, and Gayler

Reick, C. H., Raddatz, T., Brovkin, V., and Gayler, V.: Representation of natural and anthropogenic land cover change in MPI-ESM, J. Adv. Model. Earth Sy., 5, 459–482, 10.1002/jame.20022, 2013.

Rockström et al.(2009)Rockström, Steffen, Noone, Persson, Chapin, Lambin, Lenton, Scheffer, Folke, Schellnhuber, Nykvist, de Wit, Hughes, van der Leeuw, Rodhe, Sörlin, Snyder, Costanza, Svedin, Falkenmark, Karlberg, Corell, Fabry, Hansen, Walker, Liverman, Richardson, Crutzen, and Foley

Rockström, J., Steffen, W., Noone, K., Persson, A., Chapin, F. S., Lambin, E. F., Lenton, T. M., Scheffer, M., Folke, C., Schellnhuber, H. J., Nykvist, B., de Wit, C. A., Hughes, T., van der Leeuw, S., Rodhe, H., Sörlin, S., Snyder, P. K., Costanza, R., Svedin, U., Falkenmark, M., Karlberg, L., Corell, R. W., Fabry, V. J., Hansen, J., Walker, B., Liverman, D., Richardson, K., Crutzen, P., and Foley, J. A.: A safe operating space for humanity, Nature, 461, 472–475, 10.1038/461472a, 2009.

Rolinski et al.(2018)Rolinski, Müller, Heinke, Weindl, Biewald, Bodirsky, Bondeau, Boons-Prins, Bouwman, Leffelaar, te Roller, Schaphoff, and Thonicke

Rolinski, S., Müller, C., Heinke, J., Weindl, I., Biewald, A., Bodirsky, B. L., Bondeau, A., Boons-Prins, E. R., Bouwman, A. F., Leffelaar, P. A., te Roller, J. A., Schaphoff, S., and Thonicke, K.: Modeling vegetation and carbon dynamics of managed grasslands at the global scale with LPJmL 3.6, Geosci. Model Dev., 11, 429–451, 10.5194/gmd-11-429-2018, 2018.

Rosol et al.(2017)Rosol, Nelson, and Renn

Rosol, C., Nelson, S., and Renn, J.: Introduction: In the machine room of the Anthropocene, Anthropocene Rev., 4, 2–8, 10.1177/2053019617701165, 2017.

Rounsevell et al.(2012)Rounsevell, Robinson, and Murray-Rust

Rounsevell, M. D. A., Robinson, D. T., and Murray-Rust, D.: From actors to agents in socio-ecological systems models, Philos. T. Roy. Soc. B, 367, 259–269, 10.1098/rstb.2011.0187, 2012.

Rounsevell et al.(2014)Rounsevell, Arneth, Alexander, Brown, de Noblet-Ducoudré, Ellis, Finnigan, Galvin, Grigg, Harman, Lennox, Magliocca, Parker, O'Neill, Verburg, and Young

Rounsevell, M. D. A., Arneth, A., Alexander, P., Brown, D. G., de Noblet-Ducoudré, N., Ellis, E., Finnigan, J., Galvin, K., Grigg, N., Harman, I., Lennox, J., Magliocca, N., Parker, D., O'Neill, B. C., Verburg, P. H., and Young, O.: Towards decision-based global land use models for improved understanding of the Earth system, Earth Syst. Dynam., 5, 117–137, 10.5194/esd-5-117-2014, 2014.

Sakschewski et al.(2016)Sakschewski, von Bloh, Boit, Poorter, Peña-Claros, Heinke, Joshi, and Thonicke

Sakschewski, B., von Bloh, W., Boit, A., Poorter, L., Peña-Claros, M., Heinke, J., Joshi, J., and Thonicke, K.: Resilience of Amazon forests emerges from plant trait diversity, Nat. Clim. Change, 6, 1032–1036, 10.1038/nclimate3109, 2016.

Schaphoff et al.(2018)Schaphoff, von Bloh, Rammig, Thonicke, Biemans, Forkel, Gerten, Heinke, Jägermeyr, Knauer, Langerwisch, Lucht, Müller, Rolinski, and Waha

Schaphoff, S., von Bloh, W., Rammig, A., Thonicke, K., Biemans, H., Forkel, M., Gerten, D., Heinke, J., Jägermeyr, J., Knauer, J., Langerwisch, F., Lucht, W., Müller, C., Rolinski, S., and Waha, K.: LPJmL4 – a dynamic global vegetation model with managed land – Part 1: Model description, Geosci. Model Dev., 11, 1343–1375, 10.5194/gmd-11-1343-2018, 2018.

Schaphoff et al.(2025)Schaphoff, von Bloh, Rammig, Thonicke, Biemans, Forkel, Gerten, Heinke, Jägermeyr, Langerwisch, Lucht, Rolinski, Waha, Ostberg, Wirth, Fader, Drüke, Jans, Lutz, Herzfeld, Minoli, Porwollik, Stehfest, de Waal, Beringer (Erbrecht), Rost (Jachner), Gumpenberger, Heyder, Werner, Braun, Breier, Stenzel, Mathesius, Hemmen, Billing, Oberhagemann, Sakschewski, and Müller

Schaphoff, S., von Bloh, W., Rammig, A., Thonicke, K., Biemans, H., Forkel, M., Gerten, D., Heinke, J., Jägermeyr, J., Langerwisch, F., Lucht, W., Rolinski, S., Waha, K., Ostberg, S., Wirth, S. B., Fader, M., Drüke, M., Jans, Y., Lutz, F., Herzfeld, T., Minoli, S., Porwollik, V., Stehfest, E., de Waal, L., Beringer (Erbrecht), T., Rost (Jachner), S., Gumpenberger, M., Heyder, U., Werner, C., Braun, J., Breier, J., Stenzel, F., Mathesius, S., Hemmen, M., Billing, M., Oberhagemann, L., Sakschewski, B., and Müller, C.: LPJmL for InSEEDS and further copan:LPJmL models, Zenodo [code], 10.5281/zenodo.17036653, 2025.

Schellnhuber(1999)

Schellnhuber, H. J.: “Earth system” analysis and the second Copernican revolution, Nature, 402, C19–C23, 10.1038/35011515, 1999.

Schill et al.(2019)Schill, Anderies, Lindahl, Folke, Polasky, Cárdenas, Crépin, Janssen, Norberg, and Schlüter

Schill, C., Anderies, J. M., Lindahl, T., Folke, C., Polasky, S., Cárdenas, J. C., Crépin, A.-S., Janssen, M. A., Norberg, J., and Schlüter, M.: A more dynamic understanding of human behaviour for the Anthropocene, Nat. Sustain., 2, 1075–1082, 10.1038/s41893-019-0419-7, 2019.

Schleussner et al.(2024)Schleussner, Ganti, Lejeune, Zhu, Pfleiderer, Prütz, Ciais, Frölicher, Fuss, Gasser, Gidden, Kropf, Lacroix, Lamboll, Martyr, Maussion, McCaughey, Meinshausen, Mengel, Nicholls, Quilcaille, Sanderson, Seneviratne, Sillmann, Smith, Steinert, Theokritoff, Warren, Price, and Rogelj

Schleussner, C.-F., Ganti, G., Lejeune, Q., Zhu, B., Pfleiderer, P., Prütz, R., Ciais, P., Frölicher, T. L., Fuss, S., Gasser, T., Gidden, M. J., Kropf, C. M., Lacroix, F., Lamboll, R., Martyr, R., Maussion, F., McCaughey, J. W., Meinshausen, M., Mengel, M., Nicholls, Z., Quilcaille, Y., Sanderson, B., Seneviratne, S. I., Sillmann, J., Smith, C. J., Steinert, N. J., Theokritoff, E., Warren, R., Price, J., and Rogelj, J.: Overconfidence in climate overshoot, Nature, 634, 366–373, 10.1038/s41586-024-08020-9, 2024.

Schlüter et al.(2012)Schlüter, Mcallister, Arlinghaus, Bunnefeld, Eisenack, Hölker, Milner-Gulland, Müller, Nicholson, Quaas, and Stöven

Schlüter, M., Mcallister, R. R. J., Arlinghaus, R., Bunnefeld, N., Eisenack, K., Hölker, F., Milner-Gulland, E., Müller, B., Nicholson, E., Quaas, M., and Stöven, M.: New Horizons for Managing the Environment: A Review of Coupled Social-Ecological Systems Modeling, Nat. Resour. Model., 25, 219–272, 10.1111/j.1939-7445.2011.00108.x, 2012.

Schlüter et al.(2017)Schlüter, Baeza, Dressler, Frank, Groeneveld, Jager, Janssen, McAllister, Müller, Orach, Schwarz, and Wijermans

Schlüter, M., Baeza, A., Dressler, G., Frank, K., Groeneveld, J., Jager, W., Janssen, M. A., McAllister, R. R. J., Müller, B., Orach, K., Schwarz, N., and Wijermans, N.: A framework for mapping and comparing behavioural theories in models of social-ecological systems, Ecol. Econ., 131, 21–35, 10.1016/j.ecolecon.2016.08.008, 2017.

Schlüter et al.(2021)Schlüter, Lindkvist, Wijermans, and Polhill

Schlüter, M., Lindkvist, E., Wijermans, N., and Polhill, G.: Agent-based modelling, in: The Routledge Handbook of Research Methods for Social-Ecological Systems, Routledge, ISBN 978-1-003-02133-9, 2021.

Schrieks et al.(2021)Schrieks, Botzen, Wens, Haer, and Aerts

Schrieks, T., Botzen, W. J. W., Wens, M., Haer, T., and Aerts, J. C. J. H.: Integrating Behavioral Theories in Agent-Based Models for Agricultural Drought Risk Assessments, Front. Water, 3, 10.3389/frwa.2021.686329, 2021.

Schulze et al.(2017)Schulze, Müller, Groeneveld, and Grimm

Schulze, J., Müller, B., Groeneveld, J., and Grimm, V.: Agent-Based Modelling of Social-Ecological Systems: Achievements, Challenges, and a Way Forward, J. Artif. Soc. S., 20, 10.18564/jasss.3423, 2017.

Schwarz et al.(2025)Schwarz, Breier, Prawitz, Bechthold, von Bloh, Constantino, Gerten, Heitzig, Hotz, John, Müller, Rockström, and Donges

Schwarz, L., Breier, J., Prawitz, H., Bechthold, M., von Bloh, W., Constantino, S. M., Gerten, D., Heitzig, J., Hotz, R., John, L., Müller, C., Rockström, J., and Donges, J. F.: From Farm to Planet: The InSEEDS World-Earth Model for Simulating Transitions to Regenerative Agriculture, EGUsphere [preprint], 10.5194/egusphere-2025-4079, 2025.

Smith et al.(1993)Smith, Shugart, Woodward, and Burton

Smith, T. M., Shugart, H. H., Woodward, F. I., and Burton, P. J.: Plant Functional Types, in: Vegetation Dynamics & Global Change, edited by Solomon, A. M. and Shugart, H. H., Springer US, Boston, MA, 272–292, ISBN 978-1-4615-2816-6, 10.1007/978-1-4615-2816-6_14, 1993.

Springmann et al.(2018)Springmann, Clark, Mason-D’Croz, Wiebe, Bodirsky, Lassaletta, de Vries, Vermeulen, Herrero, Carlson, Jonell, Troell, DeClerck, Gordon, Zurayk, Scarborough, Rayner, Loken, Fanzo, Godfray, Tilman, Rockström, and Willett

Springmann, M., Clark, M., Mason-D’Croz, D., Wiebe, K., Bodirsky, B. L., Lassaletta, L., de Vries, W., Vermeulen, S. J., Herrero, M., Carlson, K. M., Jonell, M., Troell, M., DeClerck, F., Gordon, L. J., Zurayk, R., Scarborough, P., Rayner, M., Loken, B., Fanzo, J., Godfray, H. C. J., Tilman, D., Rockström, J., and Willett, W.: Options for keeping the food system within environmental limits, Nature, 562, 519–525, 10.1038/s41586-018-0594-0, 2018.

Steffen et al.(2011)Steffen, Persson, Deutsch, Zalasiewicz, Williams, Richardson, Crumley, Crutzen, Folke, Gordon, Molina, Ramanathan, Rockström, Scheffer, Schellnhuber, and Svedin

Steffen, W., Persson, Å., Deutsch, L., Zalasiewicz, J., Williams, M., Richardson, K., Crumley, C., Crutzen, P., Folke, C., Gordon, L., Molina, M., Ramanathan, V., Rockström, J., Scheffer, M., Schellnhuber, H. J., and Svedin, U.: The Anthropocene: From Global Change to Planetary Stewardship, AMBIO, 40, 739–761, 10.1007/s13280-011-0185-x, 2011.

Steffen et al.(2020)Steffen, Richardson, Rockström, Schellnhuber, Dube, Dutreuil, Lenton, and Lubchenco

Steffen, W., Richardson, K., Rockström, J., Schellnhuber, H. J., Dube, O. P., Dutreuil, S., Lenton, T. M., and Lubchenco, J.: The emergence and evolution of Earth System Science, Nat. Rev. Earth Environ., 1, 54–63, 10.1038/s43017-019-0005-6, 2020.

Stehfest(2014)

Stehfest, E. (Ed.): Integrated assessment of global environmental change with IMAGE 3.0: model description and policy applications, PBL Netherlands Environmental Assessment Agency, The Hague, ISBN 978-94-91506-71-0, 2014.

Stern et al.(1999)Stern, Dietz, Abel, Guagnano, and Kalof

Stern, P. C., Dietz, T., Abel, T., Guagnano, G. A., and Kalof, L.: A Value-Belief-Norm Theory of Support for Social Movements: The Case of Environmentalism, Human Ecol. Rev., 6, 81–97, https://www.jstor.org/stable/24707060 (last access: 20 May 2026), 1999.

Szalai et al.(2014)Szalai, Kiss, Kövér, and Toepfer

Szalai, M., Kiss, J., Kövér, S., and Toepfer, S.: Simulating crop rotation strategies with a spatiotemporal lattice model to improve legislation for the management of the maize pest Diabrotica virgifera virgifera, Agr. Syst., 124, 39–50, 10.1016/j.agsy.2013.10.009, 2014.

Tajfel and Turner(2004)

Tajfel, H. and Turner, J. C.: The Social Identity Theory of Intergroup Behavior, in: Political Psychology, Psychology Press, 18 pp., 10.4324/9780203505984, 2004.

Verburg and Overmars(2009)

Verburg, P. H. and Overmars, K. P.: Combining top-down and bottom-up dynamics in land use modeling: exploring the future of abandoned farmlands in Europe with the Dyna-CLUE model, Landscape Ecol., 24, 1167–1181, 10.1007/s10980-009-9355-7, 2009.

von Bloh et al.(2018)von Bloh, Schaphoff, Müller, Rolinski, Waha, and Zaehle

von Bloh, W., Schaphoff, S., Müller, C., Rolinski, S., Waha, K., and Zaehle, S.: Implementing the nitrogen cycle into the dynamic global vegetation, hydrology, and crop growth model LPJmL (version 5.0), Geosci. Model Dev., 11, 2789–2812, 10.5194/gmd-11-2789-2018, 2018.

Waha et al.(2012)Waha, van Bussel, Müller, and Bondeau

Waha, K., van Bussel, L. G. J., Müller, C., and Bondeau, A.: Climate-driven simulation of global crop sowing dates, Global Ecol. Biogeogr., 21, 247–259, 10.1111/j.1466-8238.2011.00678.x, 2012.

Wens et al.(2020)Wens, Veldkamp, Mwangi, Johnson, Lasage, Haer, and Aerts

Wens, M., Veldkamp, T. I. E., Mwangi, M., Johnson, J. M., Lasage, R., Haer, T., and Aerts, J. C. J. H.: Simulating Small-Scale Agricultural Adaptation Decisions in Response to Drought Risk: An Empirical Agent-Based Model for Semi-Arid Kenya, Front. Water, 2, 10.3389/frwa.2020.00015, 2020.

Wirth et al.(2024a)Wirth, Braun, Heinke, Ostberg, Rolinski, Schaphoff, Stenzel, von Bloh, Taube, and Müller

Wirth, S. B., Braun, J., Heinke, J., Ostberg, S., Rolinski, S., Schaphoff, S., Stenzel, F., von Bloh, W., Taube, F., and Müller, C.: Biological nitrogen fixation of natural and agricultural vegetation simulated with LPJmL 5.7.9, Geosci. Model Dev., 17, 7889–7914, 10.5194/gmd-17-7889-2024, 2024a.

Wirth et al.(2024b)Wirth, Poyda, Taube, Tietjen, Müller, Thonicke, Linstädter, Behn, Schaphoff, von Bloh, and Rolinski

Wirth, S. B., Poyda, A., Taube, F., Tietjen, B., Müller, C., Thonicke, K., Linstädter, A., Behn, K., Schaphoff, S., von Bloh, W., and Rolinski, S.: Connecting competitor, stress-tolerator and ruderal (CSR) theory and Lund Potsdam Jena managed Land 5 (LPJmL 5) to assess the role of environmental conditions, management and functional diversity for grassland ecosystem functions, Biogeosciences, 21, 381–410, 10.5194/bg-21-381-2024, 2024b.

Ye et al.(2024)Ye, Liu, Swamy, Gao, Moallemi, Rydzak, and Eker

Ye, Q., Liu, Q., Swamy, D., Gao, L., Moallemi, E. A., Rydzak, F., and Eker, S.: FeliX 2.0: An integrated model of climate, economy, environment, and society interactions, Environ. Modell. Softw., 179, 106121, 10.1016/j.envsoft.2024.106121, 2024.

Zabel et al.(2021)Zabel, Müller, Elliott, Minoli, Jägermeyr, Schneider, Franke, Moyer, Dury, Francois, Folberth, Liu, Pugh, Olin, Rabin, Mauser, Hank, Ruane, and Asseng

Zabel, F., Müller, C., Elliott, J., Minoli, S., Jägermeyr, J., Schneider, J. M., Franke, J. A., Moyer, E., Dury, M., Francois, L., Folberth, C., Liu, W., Pugh, T. A., Olin, S., Rabin, S. S., Mauser, W., Hank, T., Ruane, A. C., and Asseng, S.: Large potential for crop production adaptation depends on available future varieties, Glob. Change Biol., 27, 3870–3882, 10.1111/gcb.15649, 2021.

Zeng et al.(2025)Zeng, Brown, Raymond, Byari, Hotz, and Rounsevell

Zeng, Y., Brown, C., Raymond, J., Byari, M., Hotz, R., and Rounsevell, M.: Exploring the opportunities and challenges of using large language models to represent institutional agency in land system modelling, Earth Syst. Dynam., 16, 423–449, 10.5194/esd-16-423-2025, 2025.