Thoracic aortic aneurysm (TAA) is a localized dilatation of the aorta resulting from compromised wall composition, structure, and function, which can lead to life-threatening dissection or rupture. Several genetic mutations and predisposing factors that contribute to TAA have been studied in mouse models to characterize specific changes in aortic microstructure and material properties that result from a wide range of mechanobiological insults. By contrast, assessments of TAA progression in vivo are largely limited to measurements of aneurysm size and growth rate. It has been shown, however, that aortic geometry alone is not sufficient to predict the patient-specific progression of TAA but computational modeling of the evolving biomechanics of the aorta could predict future geometry and properties from initiating insults. In this work, we present an integrated framework to train a deep operator network (DeepONet)-based surrogate model to identify contributing factors for TAA by using synthetic finite element-based datasets of aortic growth and remodeling (G&R) resulting from prescribed mechanobiological insults. For training data, we investigate multiple types of TAA risk factors and spatial distributions within a computationally efficient constrained mixture model to generate axial-azimuthal maps of aortic dilatation and distensibility. The trained network is then capable of predicting the initial distribution and extent of the insult from a given set of dilatation and distensibility information, which in turn can be used to determine subsequent aortic geometry and mechanical properties. Two DeepONet frameworks are proposed, one trained on sparse information and one on full-field grayscale images, to gain insight into a preferred neural operator-based approach. Performance of the surrogate models is evaluated through multiple simulations carried out on insult distributions varying from fusiform (analytically defined) to complex (randomly generated). We show that this integrated continuous learning modeling approach can predict the patient-specific mechanobiological insult profile associated with any given dilatation and distensibility map with a high accuracy, particularly when based on full-field images. Our findings demonstrate the feasibility of