Europa and Ganymede, Galilean satellites of Jupiter, exhibit geologic activity in their outer H₂O ice shells that might convey material from water oceans within the satellites to their surfaces. Imagery from the Voyager and Galileo spacecraft reveal surfaces rich with tectonic deformation, including dilational bands on Europa and groove lanes on Ganymede. These features are generally attributed to the extension of a brittle ice lithosphere overlaying a possibly convecting ice asthenosphere. To explore band formation and interaction with interior oceans, we employ fully visco-elasto-plastic 2-D models of faulting and convection with complex, realistic pure ice rheologies, and track the change in dominant deformation mechanisms at each model element. In these models, material entering from below is tracked and considered to be “fossilized ocean,” ocean material that has frozen into the ice shell and evolves through geologic time. We vary fault localization mechanisms and fault annealing rates, ice shell thickness, melting-temperature ice viscosity, and the presence of pre-existing weaknesses. We allow the spontaneous brittle failure of the lithosphere through cohesion loss to localize strain onto fault surfaces, and explore the role of fault annealing timescales on terrain development. We also consider fault localization through tidal heating, frictional fault heating, and viscous dissipation. Mechanisms which act to weaken the ice shell and thin the lithosphere (e.g. vigorous convection, thinner shells, pre-existing weaknesses) tend to plastically yield to form smooth bands at high strains, and are more likely to expose fossil ocean material. In contrast, lithosphere strengthened by rapid fault annealing or increased viscosity, for example, exhibits large-scale tectonic rifting at low strains superimposed over pre-existing terrains, and inhibiting the delivery of fossil ocean material to the surface. Thus, our results identify a spectrum of extensional terrain formation mechanisms as linked to lithospheric strength, rather than any specific mechanism being unique to each type of band. This presentation focuses specifically on how fault localization assumptions affect the types of extensional terrains produced.