Abstract. Increased availability and use of 3D-rendered geological models have provided society with predictive capabilities, supporting natural resource assessments, hazard awareness, and infrastructure development. The Geological Survey of Canada, along with other such institutions, has been trying to standardize and operationalize this modelling practice. Knowing what is in the subsurface, however, is not an easy exercise, especially when it is difficult or impossible to sample at greater depths. Existing approaches for creating 3D geological models involve developing surface components that represent spatial geological features, horizons, faults, and folds, and then assembling them into a framework model as context for downstream property modelling applications (e.g. geophysical inversions, thermo-mechanical simulations, and fracture density models). The current challenge is to develop geologically reasonable starting framework models from regions with sparser data when we have more complicated geology. This study explores the problem of geological data sparsity and presents a new approach that may be useful to open up the logjam in modelling the more challenging terrains using an agent-based approach. Semi-autonomous software entities called spatial agents can be programmed to perform spatial and property interrogation functions, estimations and construction operations for simple graphical objects, that may be usable in building 3D geological surfaces. These surfaces form the building blocks from which full geological and topological models are built and may be useful in sparse-data environments, where ancillary or a priori information is available. Critical in developing natural domain models is the use of gradient information. Increasing the density of spatial gradient information (fabric dips, fold plunges, and local or regional trends) from geologic feature orientations (planar and linear) is the key to more accurate geologic modelling and is core to the functions of spatial agents presented herein. This study, for the first time, examines the potential use of spatial agents to increase gradient constraints in the context of the Loop project (https://loop3d.github.io/, last access: 1 October 2021) in which new complementary methods are being developed for modelling complex geology for regional applications. The spatial agent codes presented may act to densify and supplement gradient as well as on-contact control points used in LoopStructural (https://www.github.com/Loop3d/LoopStructural, last access: 1 October 2021) and Map2Loop (https://doi.org/10.5281/zenodo.4288476, de Rose et al., 2020). Spatial agents are used to represent common geological data constraints, such as interface locations and gradient geometry, and simple but topologically consistent triangulated meshes. Spatial agents can potentially be used to develop surfaces that conform to reasonable geological patterns of interest, provided that they are embedded with behaviours that are reflective of the knowledge of their geological environment. Initially, this would involve detecting simple geological constraints: locations, trajectories, and trends of geological interfaces. Local and global eigenvectors enable spatial continuity estimates, which can reflect geological trends, with rotational bias, using a quaternion implementation. Spatial interpolation of structural geology orientation data with spatial agents employs a range of simple nearest-neighbour to inverse-distance-weighted (IDW) and quaternion-based spherical linear rotation interpolation (SLERP) schemes. This simulation environment implemented in NetLogo 3D is potentially useful for complex-geology–sparse-data environments where extension, projection, and propagation functions are needed to create more realistic geological forms.