Abstract:We introduce a simplified ecohydrological model to quantitatively assess the resiliency and stability of ecosystems over long periods. The model couples a hydrological soil moisture balance with a set of spatiotemporal systems dynamics and agent-based algorithms to represent the interactions among several plant populations in a gridded area under different water, soil and temperature constraints. We characterize the plant populations by allometric rules (i.e., power laws for generational and reproductive times; linear approximations for water and temperature gains, losses and optimal values; and a set of intra and interspecific interaction rules based on high, optimal and low competition responses among the populations), that represent different plant phenotypes. We define the disturbances by a clearance of populations in an area within the model's domain and calculate the resiliency and stability with simple indices to determine the ability of the ecosystem to recover from a disturbance. The indices evaluated on each population and over the structure of the entire ecosystem show how the populations respond differently to disturbances, following patterns similar to those expected in nature. The model can represent the spatial and temporal succession of the ecosystem after being disturbed, suggesting how the differences in the phenotypical characteristics of plant populations can be advantageous or disadvantageous for the ecosystem recovery.Resumen: Se presenta un modelo ecohidrológico simplificado para evaluar cuantitativamente la resiliencia y estabilidad de los ecosistemas durante largos períodos. El modelo acopla el balance hidrológico de humedad del suelo con un conjunto de dinámicas espaciotemporales de sistemas y algoritmos basados en agentes para representar las interacciones entre varias poblaciones de plantas en un área de grilla bajo diferentes restricciones de agua, suelo y temperatura. Las poblaciones de plantas se caracterizan por reglas alométricas (por ejemplo, leyes de potencia para tiempos generacionales y reproductivos; aproximaciones lineales para ganancias de agua y temperatura; pérdidas y valores óptimos; y un conjunto de reglas de interacción intra e interespecíficas basadas en respuestas de competencia altas, óptimas y bajas entre las poblaciones), que representan diferentes fenotipos de plantas. Los trastornos se determinan por una remoción de las poblaciones en un área dentro del dominio del modelo. La resistencia y la estabilidad se calculan con índices simples para determinar la capacidad del ecosistema para recuperarse de una perturbación. Los índices evaluados en cada población y sobre la estructura de todo el ecosistema muestran cómo las poblaciones responden de manera diferente a los disturbios, siguiendo patrones similares a los esperados en la naturaleza. El modelo puede representar la sucesión espacial y temporal del ecosistema después de ser perturbado, lo que sugiere que las diferencias en las características fenotípicas de las poblaciones de plantas pueden ser ventaj...