Viscous overthrusting versus folding: 2-D quantitative modeling and its application to the Helvetic and Jura fold and thrust belts

摘要

This paper presents two-dimensional (2-D) numerical simulations of the shortening of a stiff viscous layer, with a pre-existing weak zone, that is embedded in a weaker viscous matrix. The model bottom remains straight during shortening and represents a detachment surface. Four deformation styles are observed that depend on the ratio of the layer and matrix viscosities, μ_L/μ_m; the thickness ratio of the bottom matrix and the overlying stiff layer, H_m/H_i, and the power-law stress exponent of the matrix, n_M. The numerical results are used to quantify the conditions for which each deformation style occurs: (1) pure shear-dominated deformation occurs for μ_l/μ_m < ～50 and n_m = 1 (i.e., linear viscous); (2) overthrusting-dominated deformation occurs for a power-law viscous matrix (n_m = 5), μ_L/μ_M > ～50 and H_m/H_l < ～ 0.5; (3) folding-dominated deformation occurs for n_M = 1, μ_L/μ_M > ～ 50, and H_m/H_l > ～ 1: and (4) folding and overthrusting occur for a power-law viscous matrix and ～0.5 < H_m/H_l < ～2. The power-law stress exponent of the stiff layer has a minor effect on the deformation style. Simulations with layers that contain two weak zones show the formation of a nappe stack. The simulations also show that multilayers can overthrust like a single layer. The change in deformation style as a function of H_m/H_l corroborates field observations from the Helvetic nappe system. The agreement of the numerical results with first-order observations from the Helvetic nappe system suggest that ductile deformation dominated this fold and thrust belt. The results further suggest that overthrusting on an effectively viscous weak layer only occurs if the rheology of the weak layer is power-law viscous.