Cortical neurodevelopment is sensitive to disruption following preterm birth, with lasting impact on cognitive outcomes. The creation of generative models of neurodevelopment could aid clinicians in identifying at-risk subjects but is complicated by the degree of subject variability in cortical folding, and significant heterogeneity in the effect of preterm birth. In this work, we propose a graph convolutional generative adversarial network (GAN) and a training scheme to simulate neonatal cortical surface developmental trajectories. The proposed model is used to smoothly modify two cortical phenotypes: post-menstrual age at scan (PMA) and gestational age at birth (GA) on data from the developing Human Connectome Project (dHCP). The synthetic images were validated with an independently trained regression network, and compared against follow up scans, indicating that the model can realistically age individuals whilst preserving subject-specific cortical morphology. Deviation between simulated healthy scans, and preterm follow up scans generated a metric of individual atypicality, which improved prediction of 18-month cognitive outcome over GA alone.