Mesozoic transtensional basin history of the Eastern Cordillera, Colombian Andes: Inferences from tectonic models

Backstripping analysis and forward modeling of 162 stratigraphic columns and wells of the Eastern Cordillera (EC), Llanos, and Magdalena Valley shows the Mesozoic Colombian Basin is marked by five lithosphere stretching pulses. Three stretching events are suggested during the Triassic–Jurassic, but additional biostratigraphical data are needed to identify them precisely. The spatial distribution of lithosphere stretching values suggests that small, narrow (<150 km), asymmetric graben basins were located on opposite sides of the paleo-Magdalena–La Salina fault system, which probably was active as a master transtensional or strike-slip fault system. Paleomagnetic data suggesting a significant (at least 10°) northward translation of terranes west of the Bucaramanga fault during the Early Jurassic, and the similarity between the early Mesozoic stratigraphy and tectonic setting of the Payandé terrane with the Late Permian transtensional rift of the Eastern Cordillera of Peru and Bolivia indicate that the areas were adjacent in early Mesozoic times. New geochronological, petrological, stratigraphic, and structural research is necessary to test this hypothesis, including additional paleomagnetic investigations to determine the paleolatitudinal position of the Central Cordillera and adjacent tectonic terranes during the Triassic–Jurassic. Two stretching events are suggested for the Cretaceous: Berriasian–Hauterivian (144–127 Ma) and Aptian–Albian (121–102 Ma). During the Early Cretaceous, marine facies accumulated on an extensional basin system. Shallow-marine sedimentation ended at the end of the Cretaceous due to the accretion of oceanic terranes of the Western Cordillera. In Berriasian–Hauterivian subsidence curves, isopach maps and paleomagnetic data imply a (>180 km) wide, asymmetrical, transtensional half-rift basin existed, divided by the Santander Floresta horst or high. The location of small mafic intrusions coincides with areas of thin crust (crustal stretching factors >1.4) and maximum stretching of the subcrustal lithosphere. During the Aptian–early Albian, the basin extended toward the south in the Upper Magdalena Valley. Differences between crustal and subcrustal stretching values suggest some lowermost crustal decoupling between the crust and subcrustal lithosphere or that increased thermal thinning affected the mantle lithosphere. Late Cretaceous subsidence was mainly driven by lithospheric cooling, water loading, and horizontal compressional stresses generated by collision of oceanic terranes in western Colombia. Triassic transtensional basins were narrow and increased in width during the Triassic and Jurassic. Cretaceous transtensional basins were wider than Triassic–Jurassic basins. During the Mesozoic, the strike-slip component gradually decreased at the expense of the increase of the extensional component, as suggested by paleomagnetic data and lithosphere stretching values. During the Berriasian–Hauterivian, the eastern side of the extensional basin may have developed by reactivation of an older Paleozoic rift system associated with the Guaicáramo fault system. The western side probably developed through reactivation of an earlier normal fault system developed during Triassic–Jurassic transtension. Alternatively, the eastern and western margins of the graben may have developed along older strike-slip faults, which were the boundaries of the accretion of terranes west of the Guaicáramo fault during the Late Triassic and Jurassic. The increasing width of the graben system likely was the result of progressive tensional reactivation of preexisting upper crustal weakness zones. Lateral changes in Mesozoic sediment thickness suggest the reverse or thrust faults that now define the eastern and western borders of the EC were originally normal faults with a strike-slip component that inverted during the Cenozoic Andean orogeny. Thus, the Guaicáramo, La Salina, Bitúima, Magdalena, and Boyacá originally were transtensional faults. Their oblique orientation relative to the Mesozoic magmatic arc of the Central Cordillera may be the result of oblique slip extension during the Cretaceous or inherited from the pre-Mesozoic structural grains. However, not all Mesozoic transtensional faults were inverted.     

A sequence stratigraphic model for an intensely bioturbated shallow-marine sandstone: the Bridport Sand Formation, Wessex Basin, UK

The Bridport Sand Formation is an intensely bioturbated sandstone that represents part of a mixed siliciclastic‐carbonate shallow‐marine depositional system. At outcrop and in subsurface cores, conventional facies analysis was combined with ichnofabric analysis to identify facies successions bounded by a hierarchy of key stratigraphic surfaces. The geometry of these surfaces and the lateral relationships between the facies successions that they bound have been constrained locally using 3D seismic data. Facies analysis suggests that the Bridport Sand Formation represents progradation of a low‐energy, siliciclastic shoreface dominated by storm‐event beds reworked by bioturbation. The shoreface sandstones form the upper part of a thick (up to 200 m), steep (2–3°), mud‐dominated slope that extends into the underlying Down Cliff Clay. Clinoform surfaces representing the shoreface‐slope system are grouped into progradational sets. Each set contains clinoform surfaces arranged in a downstepping, offlapping manner that indicates forced‐regressive progradation, which was punctuated by flooding surfaces that are expressed in core and well‐log data. In proximal locations, progradational shoreface sandstones (corresponding to a clinoform set) are truncated by conglomerate lags containing clasts of bored, reworked shoreface sandstones, which are interpreted as marking sequence boundaries. In medial locations, progradational clinoform sets are overlain across an erosion surface by thin (<5 m) bioclastic limestones that record siliciclastic‐sediment starvation during transgression. Near the basin margins, these limestones are locally thick (>10 m) and overlie conglomerate lags at sequence boundaries. Sequence boundaries are thus interpreted as being amalgamated with overlying transgressive surfaces, to form composite erosion surfaces. In distal locations, oolitic ironstones that formed under conditions of extended physical reworking overlie composite sequence boundaries and transgressive surfaces. Over most of the Wessex Basin, clinoform sets (corresponding to high‐frequency sequences) are laterally offset, thus defining a low‐frequency sequence architecture characterized by high net siliciclastic sediment input and low net accommodation. Aggradational stacking of high‐frequency sequences occurs in fault‐bounded depocentres which had higher rates of localized tectonic subsidence.     

Controls on low-pressure anatexis

Low-pressure anatexis, whereby rocks melt in place after passing through the andalusite stability field, develops under more restricted conditions than does low-pressure metamorphism. Our thermal modelling and review of published work indicate that the following mechanisms, operating alone, may induce anatexis in typical pelitic rocks without inducing wholesale melting in the lower crust: (i) magmatic advection by pervasive flow; (ii) crustal-scale detachment faulting; and (iii) the presence of a high heat-producing layer. Of these, only magmatic advection by pervasive flow and crustal-scale detachment faulting have been shown quantitatively to provide sufficient heat to cause widespread melting. Combinations of the above mechanisms with pluton-scale magmatic advection, shear heating, removal of the lithospheric mantle, or with each other provide additional means of developing suitable high temperatures at shallow crustal levels to generate low-pressure anatexis.     

3-D seismic structure of the source area of the 1993 Latur, India, earthquake and its implications for rupture nucleations

The Latur earthquake (M_w 6.1) of 29 September 1993 is a rare stable continental region (SCR) earthquake that occurred on a previously unknown blind fault. In this study, we determined detailed three-dimensional (3-D) P- and S-wave velocity (V_p, V_s) and Poisson's ratio (σ) structures by inverting the first P- and S-wave high-quality arrival time data from 142 aftershocks that were recorded by a network of temporary seismic stations. The source zone of the Latur earthquake shows strong lateral heterogeneities in V_p, V_s and σ structures, extending in a volume of about 90 × 90 × 15 km³. The mainshock occurred within, but near the boundary, of a low-V_p, high-V_s and low-σ zone. This suggests that the structural asperities at the mainshock hypocenter are associated with a partially fluid-saturated fractured rock in a previously unknown source zone with intersecting fault surfaces. This might have triggered the 1993 Latur mainshock and its aftershock sequence. Our results are in good agreement with other geophysical studies that suggest high conductivity and high concentration of radiogenic helium gas beneath the source zone of the Latur earthquake. Our study provides an additional evidence for the presence of fluid related anomaly at the hidden source zone of the Latur earthquake in the SCR and helps us understand the genesis of damaging earthquakes in the SCR of the world.