Evaluating the Rheological Controls on Topography Development During Craton Stabilization: Objective Approaches to Comparing Geodynamic Models

Abstract

Surface topography is an important yet largely neglected aspect of the early evolution of cratons. The lateral accretion of cratonic nuclei inevitably forms orogenic belts that subsequently provide a sediment source for large, resource‐rich intracratonic basins, but to date, geodynamic models have focused exclusively on lithospheric root processes. Here we use two‐dimensional thermal‐mechanical models to study the topography and lithospheric deformation during 50 Myr of compression of a cratonic nucleus, to simulate the lateral accretion phase of craton growth in the Neoarchean. Although the cratonic nucleus thickens slightly during the compression phase, most of the deformation occurs in the regions adjacent to the nucleus that have weaker lithosphere. Here, crustal thickness triples developing high topography in excess of 10 km without active erosion. Models with different initial rheological parameters will have different final topography and lithosphere geometry, but in general it is difficult to shorten and deform the depleted cratonic nucleus, unless there are significantly weak heterogeneities in the mantle lithosphere. We apply two quantitative analysis techniques to objectively evaluate a multitude of model outputs. Cross‐correlation clustering (CCC) measures the degree of similarity between topography profiles and categorizes models based on the general topographic character. Six different topography families are possible in the context of our models and crustal strength is the most important parameter affecting the shape. From principal component analysis (PCA) we identify four dominant lithosphere geometries. When used together, these two methods provide distinct yet complementary information about the surface and subsurface deformation features in our models.