Two-dimensional numerical modeling of tectonic and metamorphic histories at active continental margins

The evolution of an active continental margin is simulated in two dimensions, using a finite difference thermomechanical code with half-staggered grid and marker-in-cell technique. The effect of mechanical properties, changing as a function of P and T, assigned to different crustal layers and mantle materials in the simple starting structure is discussed for a set of numerical models. For each model, representative P–T paths are displayed for selected markers. Both the intensity of subduction erosion and the size of the frontal accretionary wedge are strongly dependent on the rheology chosen for the overriding continental crust. Tectonically eroded upper and lower continental crust is carried down to form a broad orogenic wedge, intermingling with detached oceanic crust and sediments from the subducted plate and hydrated mantle material from the overriding plate. A small portion of the continental crust and trench sediments is carried further down into a narrow subduction channel, intermingling with oceanic crust and hydrated mantle material, and to some extent extruded to the rear of the orogenic wedge underplating the overriding continental crust. The exhumation rates for (ultra)high pressure rocks can exceed subduction and burial rates by a factor of 1.5–3, when forced return flow in the hanging wall portion of the self-organizing subduction channel is focused. The simulations suggest that a minimum rate of subduction is required for the formation of a subduction channel, because buoyancy forces may outweigh drag forces for slow subduction. For a weak upper continental crust, simulated by a high pore pressure coefficient in the brittle regime, the orogenic wedge and megascale melange reach a mid- to upper-crustal position within 10–20 Myr (after 400–600 km of subduction). For a strong upper crust, a continental lid persists over the entire time span covered by the simulation. The structural pattern is similar in all cases, with four zones from trench toward arc: (a) an accretionary complex of low-grade metamorphic sedimentary material; (b) a wedge of mainly continental crust, with medium-grade HP metamorphic overprint, wound up and stretched in a marble cake fashion to appear as nappes with alternating upper and lower crustal provenance, and minor oceanic or hydrated mantle interleaved material; (c) a megascale melange composed of high-pressure and ultrahigh-pressure metamorphic oceanic and continental crust, and hydrated mantle, all extruded from the subduction channel; (d) zone represents the upward tilted frontal part of the remaining upper plate lid in the case of a weak upper crust. The shape of the P–T paths and the time scales correspond to those typically recorded in orogenic belts. Comparison of the numerical results with the European Alps reveals some similarities in their gross structural and metamorphic pattern exposed after collision. A similar structure may be developed at depth beneath the forearc of the Andes, where the importance of subduction erosion is well documented, and where a strong upper crust forms a stable lid.