The ubiquity of tectonic features formed in extension, and the apparent absence of ones formed in contraction, has led to the hypothesis that Ganymede has undergone global expansion in its past. Determining the magnitude of such expansion is challenging however, and extrapolation of locally or regionally inferred strains to global scales often results in strain estimates that exceed those based on global constraints. Here we use numerical simulations of groove terrain formation to develop a strain history for Ganymede that is generally consistent at local, regional, and global scales. These simulations reproduce groove-like amplitudes, wavelengths, and average slopes at modest regional extensions (10–15%). The modest strains are more consistent with global constraints on Ganymede’s expansion. Yet locally, we also find that surface strains can be much larger (30–60%) in the same simulations, consistent with observations of highly-extended impact craters. Thus our simulations satisfy both the smallest-scale and largest-scale inferences of strain on Ganymede. The growth rate of the topography is consistent with (or exceeds) predictions of analytical models, and results from the use of a non-associated plastic rheology that naturally permits localization of brittle failure (plastic strain) into linear fault-like shear zones. These fault-like zones are organized into periodically-spaced graben-like structures with stepped, steeply-dipping faults. As in previous work, groove amplitudes and wavelengths depend on both the imposed heat flux and surface temperature, but because our brittle strength increases with depth, we find (for the parameters explored) that the growth rate of topography is initially faster for lower heat flows. We observe a transition to narrow rifting for higher heat flows and larger strains, which is a potential pathway for breakaway margin or band formation.