Learned Human Aortic Morphodynamics: A Dynamical-Systems Model of Post-EVAR Sac Remodeling

Biological tissue reorganizes in response to changes in its mechanical environment, yet quantitative models of that reorganization are almost always built where observations are dense. At organ-scale in living humans the observational data is sparse and noisy: a few irregularly spaced images per individual, acquired with variable instrumentation, with loss to follow-up and patient attrition. We ask whether governing equations for organ-scale remodeling can be recovered from this regime, and pose what those equations may reveal about the underlying biology. We study the human abdominal aorta following endovascular aneurysm repair (EVAR), in which a stent graft excludes the aneurysm sac from arterial flow and thereby changes the mechanical loading of the sac wall. Each computed tomography (CT) scan is reduced to a two-dimensional state: surface area A (size) and fluctuation in integrated Gaussian curvature [Formula] (shape), both normalized to a non-pathological aortic population. Across 100 patients and 220 scans pooled from two international centers, we use Z-SINDy, a sparse identification method with statistical-mechanical uncertainty quantification, to infer affine systems of ordinary differential equations (ODEs) governing [Formula] for anatomy that remodels successfully (regressing sacs) and anatomy that does not (stable sacs). Both outcome classes approach stable fixed points, but at different locations in the morphological state space and on different timescales. Regressing anatomy approaches [Formula] = (2.0,1.0) with characteristic timescales of 1.3 and 4.6 years; stable anatomy approaches its fixed point of (3.8, 2.5) with an order of magnitude slower timescale (10.9 and 16.6 years). Embedding the learned equations within Bayesian classifiers demonstrates that trajectory-based information resolves patient outcome classification earlier than positional evidence alone. Toward model utilization in future studies, these advantages are also observed when empirically observed loss to follow-up is carried through the priors. Recovering the equations that govern remodeling, rather than tracking its instantaneous state, therefore provides a foundation for clinical models built from the sparse records that routine imaging surveillance actually produces.