ALBERT

Description: New biostratigraphic zonations, core descriptions, sandstone petrography, facies analysis, and seismic information are compared with published detrital and bedrock geo- and thermochronology to build a Cenozoic paleogeographic reconstruction of the Andean retroarc region of Colombia, encompassing the ancestral Central Cordillera, Middle Magdalena Valley, Eastern Cordillera, and Llanos basin. We identify uplifted sediment source areas, provenance domains, depositional environments, and thickness changes to propose a refined paleogeographic evolution of eastern Colombia. We conclude that Cenozoic evolution of the northernmost Andes includes (1) a period of contractional deformation focused in the Central Cordillera and Middle Magdalena Valley that may have started by the Late Cretaceous, although thermochronological data points to maximum shortening and exhumation during the late Paleocene; (2) a period of slower deformation rates or even tectonic quiescence during the middle Eocene; and (3) a renewed phase of contractional deformation from the late Eocene to the Pleistocene/Holocene expressed in provenance, bedrock thermochronology, and increased subsidence rates in the Llanos foreland. The sedimentary response in the Llanos foreland basin is controlled by source area proximity, exhumation and shortening rates, relationships between accommodation and sediment supply, as well as potential paleoclimate forcing. This new reconstruction changes the picture of Cenozoic basin evolution offered by previous reconstructions, providing an updated chronology of deformation, which is tied to a more precise understanding of basin evolution.

Description: Surface uplift of the Garzón Massif in the northern Andes formed a critical orographic barrier (2500–3000 m elevation) that generated a deep rain shadow and strongly influenced the evolution of the largest river systems draining northern South America. This basement massif and its corresponding foreland basement high define the headwaters and drainage divides of the Amazon, Orinoco, and Magdalena Rivers. Despite its pivotal role, the exhumation history of the Garzón Massif and its relationships to the structural evolution of the broader Eastern Cordillera fold-thrust belt remain unclear. The northern Andes underwent major Cenozoic shortening, with considerable thin-skinned and thick-skinned deformation and topographic development in the Eastern Cordillera focused during late Miocene time. On the basis of widespread coarse-grained nonmarine sedimentation, previous studies have inferred that uplift of the Garzón Massif began during the late Miocene, coincident with rapid elevation gain elsewhere in the Eastern Cordillera. We take an integrated, multiproxy approach to better reconstruct Andean topographic growth and distinguish between exhumation and surface uplift of the Garzón Massif. We present new U-Pb detrital zircon provenance data, sandstone petrographic data, and paleoprecipitation data from upper Miocene clastic fill of the Neiva Basin within the adjacent Upper Magdalena Valley of the modern hinterland. In addition, six new apatite fission track (AFT) ages from the central segment of the northeast-trending Garzón Massif (Jurassic granite and Proterozoic gneiss and schist) directly constrain its Neogene exhumation history. The results indicate that early exhumation may have initiated by ca. 12.5 Ma, but a substantial orographic barrier was not fully established until ca. 6–3 Ma, when 〉1 km/m.y. of material was exhumed. Thermal history modeling of the AFT data suggests diminished exhumation thereafter (3–0 Ma), during latest Cenozoic oblique Nazca–South America convergence. This exhumation history is consistent with paleontological data suggesting late Miocene divergence of the three river systems, with associated transcontinental drainage of the Amazon River.

Description: New results from two novel lipid biomarker–based proxies help to constrain the late Cenozoic surface elevation history of the Eastern Cordillera in the tropical northern Andes of Colombia. Previous well-known studies have suggested rapid latest Miocene–Pliocene (6–3 Ma) uplift on the basis of an abrupt shift in pollen species assemblages within sedimentary basin fill of the elevated Bogotá Plateau. From resampling of these original study localities, we provide a revised chronology based on magnetic polarity stratigraphy, and we evaluate paleotemperatures using the MBT'/CBT indices (methylation of branched tetraethers/cyclization of branched tetraethers) and hydrogen isotopic composition (D) of leaf waxes as two independent proxies of past surface elevation. Reconstructed paleotemperatures from the MBT'/CBT proxy show a more gradual cooling trend from ca. 7.6 Ma to present, consistent with less than 1000 m of elevation gain since latest Miocene–Pliocene time and in agreement with geologic evidence for accelerated shortening and exhumation at this time. The leaf wax isotopic data, on the other hand, lack a systematic trend, potentially due to fractionation changes associated with uplift-induced turnover in floral populations. Such changes could obscure isotopic variations in meteoric water, suggesting that leaf-wax isotopic compositions may not provide a direct proxy for elevation in this particular situation involving uplift in a heterogeneous tropical environment. More promising in this case is the MBT'/CBT proxy, which may offer a suitable alternative to carbonate-based paleoelevation proxies in tropical regions where year-round high precipitation rates prevent formation of soil carbonates.

Description: Sedimentologic and provenance analyses for the Qaidam Basin in the northern Tibetan Plateau help to elucidate the stratigraphic signatures of initial deformation and exhumation in basin-bounding ranges. The basin recorded sedimentary transitions in response to uplift and unroofing of several distinctive source regions. Along the NE basin margin, a detrital record of exhumation and basin isolation is preserved in the 6200-m-thick Cenozoic succession at the Dahonggou anticline. An up-section shift from axial fluvial and marginal lacustrine deposition to transverse fluvial sedimentation suggests progradation and increasingly proximal sediment sources, reflecting activation and advance of crustal deformation. Provenance results from sandstone petrology, U-Pb geochronology, and heavy mineral analyses indicate initial late Paleocene–early Eocene derivation from igneous, metamorphic, and sedimentary sources, consistent with Permian–Triassic arc rocks dominating the southern (Kunlun Shan) or southwestern (Qimen Tagh) basin margins. Up-section variations in sediment composition and detrital zircon U-Pb age distributions are attributed to Eocene–Oligocene derivation from lower Paleozoic and Mesozoic igneous and metamorphic rocks of the central to northern Qilian Shan–Nan Shan. Disappearance of igneous sources and persistence of metamorphic sources are consistent with derivation from the southern Qilian Shan–Nan Shan during early–middle Miocene shortening along the frontal Nan Shan–North Qaidam thrust belt. These results are supported by paleocurrent analyses revealing an Eocene shift from roughly E-directed (axial) to SW-directed (transverse) dispersal of sediment. Variations in lithofacies, composition, U-Pb ages, and paleoflow are consistent with late Paleocene–early Eocene exhumation in the Kunlun Shan followed by middle Eocene–middle Miocene exhumation in the Qilian Shan–Nan Shan. The up-section disappearance and reappearance of diagnostic U-Pb age populations can be associated with progressive unroofing of multiple thrust sheets, successive input of sedimentary and magmatic sources, and southward encroachment of Qilian Shan–Nan Shan shortening into the Qaidam Basin. The sedimentary record presented here indicates that during the Paleogene, the unified Qaidam-Tarim Basin was partitioned and uplifted as it was incorporated into the growing Tibetan Plateau. Comparison with basins on and surrounding the Tibetan Plateau suggests that basement strength and lateral homogeneity, and formation of syndepositional structural dams are among the primary controls on formation of giant sedimentary basins.

Description: Detrital zircon U-Pb geochronology has been used extensively to develop provenance histories for surface outcrops of key stratigraphic localities within sedimentary basins. However, many basins lack sufficiently continuous and widespread exposures of complete successions to evaluate proposed long-term tectonic histories, stratigraphic correlations, and paleodrainage patterns within individual basins. Here, we demonstrate the utility of subsurface detrital zircon U-Pb analysis by integrating ages from three key wells (21 subsurface samples) with previously reported data from six exposed intervals (90 surface samples) within a single basin. Samples from the 5–10-km-thick clastic successions span several structural blocks over an ~300 x 50 km swath of the Middle Magdalena Valley Basin, a north-trending intermontane basin in the northern Andes of Colombia. Available U-Pb age distributions for modern rivers highlight the distinctive signatures of several competing sediment sources, including two major contiguous ranges (Central Cordillera and Eastern Cordillera) and two localized block uplifts (Santander Massif and San Lucas range). U-Pb results from Jurassic through Neogene stratigraphic units spanning the nine surface and subsurface sites, including several type localities, enable comparisons of provenance shifts at specific sites and spatial variations among key stratigraphic intervals across multiple sites. Distinctive age populations for the Andean magmatic arc, retroarc fold-thrust belt, and South American craton facilitate correlation of stratigraphic units and reconstruction of the long-term provenance and tectonic evolution of the Middle Magdalena Valley Basin. Nearly all surface and subsurface localities show up-section changes in age spectra consistent with (1) Jurassic growth of extensional subbasins fed by local igneous sources, (2) Cretaceous deposition in an extensive postrift setting, and (3) protracted Cenozoic growth of basin-bounding ranges during Andean crustal shortening. Subsurface samples augment surface samples, highlighting their utility in developing regional source-to-sink relationships, the timing of paleodrainage integration, and tectonic reconstructions. Provenance shifts of mid-Paleocene and latest Eocene–earliest Oligocene age are consistent with incipient uplift of the flanking Central Cordillera and Eastern Cordillera, respectively. However, a well-documented phase of latest Paleocene–middle Eocene beveling of basement uplifts in the Middle Magdalena Valley Basin appears to be largely aliased in the detrital record. Moreover, despite the proximity of the magmatic arc, there is insufficient syndepositional evidence for a proposed Paleogene pulse of magmatism and, in this case, limited utility of U-Pb ages in pinpointing precise depositional (stratigraphic) ages. U-Pb age spectra for Oligocene through Pliocene basin fill underscore complex along-strike (north-south) and cross-strike (east-west) variations reflective of compartmentalized transverse deposystems demarcated by point-source contributions from the Central Cordillera and Eastern Cordillera. The late Miocene appearance of 100–0 Ma grains and a regional switch to broad, multimodal age distributions suggest the initial integration of the longitudinal proto–Magdalena River, linking the Middle Magdalena Valley Basin with southern headwaters in the Upper Magdalena Valley and likely driving increased sedimentation rates farther north in the offshore Magdalena submarine fan of the southern Caribbean margin.

Description: The increase in detrital geochronological data presents challenges to existing approaches to data visualization and comparison, and highlights the need for quantitative techniques able to evaluate and compare multiple large data sets. We test five metrics commonly used as quantitative descriptors of sample similarity in detrital geochronology: the Kolmogorov-Smirnov (K-S) and Kuiper tests, as well as Cross-correlation, Likeness, and Similarity coefficients of probability density plots (PDPs), kernel density estimates (KDEs), and locally adaptive, variable-bandwidth KDEs (LA-KDEs). We assess these metrics by applying them to 20 large synthetic data sets and one large empirical data set, and evaluate their utility in terms of sample similarity based on the following three criteria. (1) Similarity of samples from the same population should systematically increase with increasing sample size. (2) Metrics should maximize sensitivity by using the full range of possible coefficients. (3) Metrics should minimize artifacts resulting from sample-specific complexity. K-S and Kuiper test p-values passed only one criterion, indicating that they are poorly suited as quantitative descriptors of sample similarity. Likeness and Similarity coefficients of PDPs, as well as K-S and Kuiper test D and V values, performed better by passing two of the criteria. Cross-correlation of PDPs passed all three criteria. All coefficients calculated from KDEs and LA-KDEs failed at least two of the criteria. As hypothesis tests of derivation from a common source, individual K-S and Kuiper p-values too frequently reject the null hypothesis that samples come from a common source when they are identical. However, mean p-values calculated by repeated subsampling and comparison (minimum of 4 trials) consistently yield a binary discrimination of identical versus different source populations. Cross-correlation and Likeness of PDPs and Cross-correlation of KDEs yield the widest divergence in coefficients and thus a consistent discrimination between identical and different source populations, with Cross-correlation of PDPs requiring the smallest sample size. In light of this, we recommend acquisition of large detrital geochronology data sets for quantitative comparison. We also recommend repeated subsampling of detrital geochronology data sets and calculation of the mean and standard deviation of the comparison metric in order to capture the variability inherent in sampling a multimodal population. These statistical tools are implemented using DZstats, a MATLAB-based code that can be accessed via an executable file graphical user interface. It implements all of the statistical tests discussed in this paper, and exports the results both as spreadsheets and as graphic files.

Description: Sedimentologic, provenance, geochronologic, and magnetostratigraphic results from clastic nonmarine deposits in the northern Altiplano Plateau of southern Peru (14–15°S) demonstrate late Eocene–Oligocene (37–26 Ma) accumulation of the 〉4-km-thick San Jerónimo (Puno) Group within a retroarc foreland basin related to early Andean shortening and crustal thickening. Punctuated Oligocene (29–26 Ma) displacement along deep-seated contractional structures, as revealed by growth stratal relationships, abruptly partitioned this regional flexural basin and established the structural boundaries of the smaller intermontane Ayaviri Basin, which continued to evolve in a hinterland setting during late Oligocene–Miocene shortening. This brief episode of shortening along the Altiplano–Eastern Cordillera boundary is correlated with exceptionally rapid sediment accumulation (〉1100–1800 m/m.y.), tightly constrained to 30–28 Ma on the basis of U-Pb geochronology and magnetic polarity stratigraphy. Provenance data from detrital zircon U-Pb age populations and sandstone compositions indicate derivation from a complex belt of Paleogene shortening and probable basin inversion in the Western Cordillera that was subsequently overprinted by Andean arc magmatism. This early Andean zone is interpreted as the along-strike continuation of the better-exposed Marañon fold-thrust belt to the north (5–13°S) and a proposed belt of shortening to the south along the Chilean Precordillera and Western Cordillera of Bolivia and northern Argentina (17–25°S). Subsequent focusing of late Oligocene shortening along the Eastern Cordillera–Altiplano boundary may have been linked to shallowing of the subducting slab and potential reactivation of crustal anisotropies.

Description: Evidence for a long-lived, high-elevation plateau contradicts the hypothesis that uplift of the southern Tibetan Plateau was a driver of late Miocene–Pliocene changes in ecology and monsoon strength. We illuminate the mechanisms underlying late Miocene–Pliocene environmental changes using a multi-proxy record from the Zhada Basin, southwestern Tibetan Plateau, between ca. 9.2 and 2.3 Ma. An increase in mean carbonate 18 O c and 13 C c values, decrease in grain size, and onset of lacustrine deposition at 6.0 Ma is attributed to local tectonic damming and transition from a through-flowing fluvial system to a terminal lake. This is followed at 3.5 Ma by a decrease in lake size indicated by synchronous increase in grain size, progradation of lake-margin depositional systems, increase in 18 O c values at lake-central locations, and divergence in 18 O c values between inflowing and lake water. We attribute lake shrinking to decreasing Indian summer monsoon precipitation due to the combined effects of decreased southern Tibetan Plateau elevations and cooler Indian Ocean temperatures. We attribute increased grain size and 18 O c value variability to changes in fluvial discharge due to increased sensitivity to orbital forcing, itself possibly coupled with onset of glacial conditions at high elevations. These mechanisms link regional tectonic events to Pliocene changes in climate and environmental conditions, including erosion, granularity, sediment accumulation rates, and potentially biological turnover on and around the Tibetan Plateau.