ALBERT

Description: Unravelling the long-term evolution of the subglacial landscape of Antarctica is vital for understanding past ice sheet dynamics and stability, particularly in marine-based sectors of the ice sheet. Here we model the evolution of the bedrock topography beneath the Recovery catchment, a sector of the East Antarctic Ice Sheet characterized by fast-flowing ice streams that occupy overdeepened subglacial troughs. We use 3-D flexural models to quantify the effect of erosional unloading and mechanical unloading associated with motion on border faults in driving isostatic bedrock uplift of the Shackleton Range and Theron Mountains, which are flanked by the Recovery, Slessor, and Bailey ice streams. Inverse spectral (free-air admittance) and forward modeling of topography and gravity anomaly data allow us to constrain the effective elastic thickness of the lithosphere (Te) in the Shackleton Range region to ~20 km. Our models indicate that glacial erosion, and the associated isostatic rebound, has driven 40–50% of total peak uplift in the Shackleton Range and Theron Mountains. A further 40–50% can be attributed to motion on normal fault systems of inferred Jurassic and Cretaceous age. Our results indicate that the flexural effects of glacial erosion play a key role in mountain uplift along the East Antarctic margin, augmenting previous findings in the Transantarctic Mountains. The results suggest that at 34 Ma, the mountains were lower and the bounding valley floors were close to sea level, which implies that the early ice sheet in this region may have been relatively stable.

Description: SUMMARY Recent modelling studies have shown that laboratory-derived rheology is too strong to reproduce observations of flexure at the Hawaiian Islands, while the same rheology appears consistent with outer rise—trench flexure at circum-Pacific subduction zones. Collectively, these results indicate that the rheology of an oceanic plate boundary is stronger than that of its interior, which, if correct, presents a challenge to understanding the formation of trenches and subduction initiation. To understand this dilemma, we first investigate laboratory-derived rheology using fully dynamic viscoelastic loading models and find that it is too strong to reproduce the observationally inferred elastic thickness, Te, at most plate interior settings. The Te can, however, be explained if the yield stress of low-temperature plasticity is significantly reduced, for example, by reducing the activation energy from 320 kJ mol−1, as in Mei et al., to 190 kJ mol−1 as was required by previous studies of the Hawaiian Islands, implying that the lithosphere beneath Hawaii is not anomalous. Second, we test the accuracy of the modelling methods used to constrain the rheology of subducting lithosphere, including the yield stress envelope (YSE) method, and the broken elastic plate model (BEPM). We show the YSE method accurately reproduces the model Te to within ∼10 per cent error with only modest sensitivity to the assumed strain rate and curvature. Finally, we show that the response of a continuous plate is significantly enhanced when a free edge is introduced at or near an edge load, as in the BEPM, and is sensitive to the degree of viscous coupling at the free edge. Since subducting lithosphere is continuous and generally mechanically coupled to a sinking slab, the BEPM may falsely introduce a weakness and hence overestimate Te at a trench because of trade-off. This could explain the results of recent modelling studies that suggest the rheology of subducting oceanic plate is stronger than that of its interior. However, further studies using more advanced thermal and mechanical models will be required in the future in order to quantify this.