Thermochronology and exhumation rates of granitic intrusions at Mesa Central,Mexico

ABSTRACT

The Mesa Central of Mexico (MC) is an elevated plateau located 2000 m above sea level in central Mexico, where intrusions outcrop that register the history of exhumation-erosion occurring during the Late Cretaceous-Paleogene. The tectonic history of the region records formation of the Late Cretaceous-Paleogene ‘Mexican orogen’; this was followed by extension of the entire region and several plutons were then exhumed. The age and magnitude of the crustal uplift and erosion occurring during exhumation has not been addressed to date. Therefore, this study reports the crystallization and cooling ages of two plutons, the Tesorera Granodiorite and the Comanja Granite, and estimates their emplacement depths. Based on these data, the exhumation age of the Tesorera Granodiorite is estimated to be between ~73 Ma and ~63 Ma at an exhumation rate of ~528 m/m. y. and that of the Comanja Granite is 52 Ma and 48 Ma at an exhumation rate of ~2500 m/m. y. Exhumation-erosion event of the Tesorera Granodiorite was located on the trace of the San Luis-Tepehuanes Fault System and that of the Comanja Granite on the a trace of the El Bajío Fault System. Furthermore, the high exhumation rate in the Comanja Granite suggests that gravitational collapse played an important role during exhumation.     

Spatial Distribution,Scaling and Self-similar Behavior of Fracture Arrays in the Los Planes Fault,Baja California Sur,Mexico

We present a case of detailed analysis of fracture arrays spanning four orders of magnitude in length; all of them measured at a single natural site by acquiring images at progressively larger scales. There is a high dispersion of cumulative-length exponents, box dimensions and fracture densities. However, the fractal analysis supports the fractal nature of fracture arrays. Our data indicate the existence of an upper limit for the density parameters, as similarly reported by other authors. We prove that box dimension is in inverse relation with fracture concentration and in direct relation with fracture density. These relations are also observed in our data and additionally there is an upper limit for the box dimensions. We interpret the dispersion in our results as more fundamental than methodological problems. It represents a truncation in the complete evolution of the fracture systems because in natural cases strain initiates overprinting of previous fracture arrays. Considering that larger fractures accommodate strain more efficiently than small fractures, the generation of small fractures is inhibited in the presence of pre-existing larger fractures. Maximum values of fracture density prevent accommodating an excess of strain in a single or restricted range of scales; we claim this condition produces migration of fracturing to larger scales originating fracture scaling.     

Origin and tectonic interpretation of multiple fault patterns

It is shown that in two-dimensional and three-dimensional deformation accommodated by fracture, the symmetry of the fault patterns is an intrinsic attribute because it reflects the symmetry of either stress or strain tensors. The deformation accommodated by sliding along pre-existing planes, when there is kinematic interaction between that planes, forms multiple fault pattern and multiple slickenline sets during a single deformation event. These fault patterns have no restrictions with respect to symmetry, number of fault sets or fault orientation.

The kinematic analysis developed here shows that an interacting system is formed by two cross cutting faults and three slickenlines. One slickenline must be parallel to the intersection line between the planes. Also, it is demonstrated that the slickenlines generally do not correspond to the shear stress solution on the planes. Thus, the interaction between planes does not satisfy the assumption of parallelism between shear stress and slip vector. We conclude that the inversion methods to calculate paleostress tensors can lead to erroneous interpretations in structurally complex zones with many pre-existing planes of weakness.

We propose four possibilities to form multiple fault patterns: (1) two or more events of faulting obeying Coulomb's law with a change of orientation of the principal stresses in each event; (2) reactivation of non-interacting planes according to the Bott (1959) model; (3) one three-dimensional strain event that obeys the “Slip Model”; this mechanism will form an orthorhombic four-fault pattern and two slickenline sets in a single strain event; and (4) one or more events obeying the interacting block model proposed here, with or without rotation of the principal stresses. We propose the last origin as the most common in continental regions.