10.07.2015 Views

Quantitative structural analyses and numerical modelling of ...

Quantitative structural analyses and numerical modelling of ...

Quantitative structural analyses and numerical modelling of ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

LEXA ET AL.: COLLISION IN WEST CARPATHIANS 5 - 11Figure 7. (a) Initial geometry <strong>and</strong> boundary conditions <strong>of</strong> the <strong>numerical</strong> model. The arrow indicates thevelocity <strong>and</strong> trajectory <strong>of</strong> the indenter northern margin, with v = 0; zero Dirichlet boundary condition forvelocity; outflow, zero Neumann boundary condition for velocity. (b) Finite strain pattern developed inweak zone after 1 Myr <strong>of</strong> shortening. Distribution <strong>of</strong> strain intensity expressed in D value <strong>and</strong> orientationpattern <strong>of</strong> XY plane <strong>of</strong> finite strain ellipsoid. The foliations trajectories are shown by lines <strong>and</strong> theorientation <strong>of</strong> lineation is expressed by the color <strong>of</strong> foliation trace. White line corresponds to horizontallineation <strong>and</strong> black line to vertical one. (c) Finite strain pattern developed in weak zone after 3 Myr <strong>of</strong>shortening. Distribution <strong>of</strong> strain intensity expressed in D value <strong>and</strong> orientation pattern <strong>of</strong> XY plane <strong>of</strong>finite strain ellipsoid. (d) Distribution <strong>of</strong> finite strain symmetry expressed in K value.ments within the Vepor-Gemer system are responsible for itsinternal deformation. Important observation is that thedeformation intensity in central part <strong>of</strong> the Gemer Unitvanishes to the north. In this area, the Vepor basement ispresent (covered by Tertiary sediments) but no increase inCretaceous deformation intensity has been observed. Thismeans that this Vepor segment does not creating importantdeformation resulting from possible movement to the south.Therefore we suggest that the Vepor promontories did notmove actively to the south, <strong>and</strong> extreme deformation alongwestern <strong>and</strong> eastern Vepor promontories was imposed bygenerally northward flow <strong>of</strong> weak material. In conclusion,the only differentially moving rigid block is the northwardthrusted part <strong>of</strong> sub-Gemer basement. All other basementunits can be further considered kinematically fixed.[34] Our field studies showed that apart from GCF, anintense deformation was concentrated along southeasternedge <strong>of</strong> the western Vepor promontory producing TGSZ <strong>and</strong>also along southwestern edge <strong>of</strong> the eastern Vepor promontoryresponsible for the origin <strong>of</strong> EGT (Figure 6b). Thedevelopment <strong>of</strong> TGSZ probably results from a major changein mutual translation direction <strong>of</strong> southern sub-Gemer block<strong>and</strong> western Vepor promontory due to their oblique collisionat deeper crustal levels. Localized transpressional deformationin upper crustal levels is a typical expression <strong>of</strong> obliqueconvergence in many active regions, e.g., San Andreas faultzone [Teyssier <strong>and</strong> Tik<strong>of</strong>f, 1997], Sumatra [Tik<strong>of</strong>f <strong>and</strong>Teyssier, 1994], or Alpine fault in New Zeal<strong>and</strong> [Teyssieret al., 1995].4.2. Numerical Modeling <strong>of</strong> ProgressiveDeformation <strong>of</strong> the Gemer Unit[35] The presented <strong>numerical</strong> approach enables to modelthe deformation in a weak zone surrounded by rigid blocksor free boundaries. The approach is based on the thinviscous sheet approximation being similar to that one usedby Engl<strong>and</strong> et al. [1985] for modeling the deformation <strong>of</strong>the whole lithosphere. We assume a horizontal weak tabulardomain to have been subjected to flow with no tractions attop <strong>and</strong> bottom surface. We consider vertical gradients <strong>of</strong>the horizontal velocity to be negligible, which allows us tointegrate the equations <strong>of</strong> motion over the vertical dimension<strong>and</strong> to work with vertical averages <strong>of</strong> stress <strong>and</strong> strainrates. When linear constitutive relation between stress <strong>and</strong>strain rates is considered, the procedure leads to a system <strong>of</strong>elliptic partial differential equations for two horizontalvelocity components (see Appendix A). The system issolved by the finite element method, with the Dirichlet<strong>and</strong> Neumann boundary conditions applied to segments <strong>of</strong>the domain boundaries corresponding to the describedgeological settings (rigid indenter, free inflow or outflow95

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!