Shari A. Kelley

Extension and basin formation in the southern Andes caused by increased convergence rate: A mid‐Cenozoic trigger for the Andes

The southern Andes between 33° and 45°S latitude are characterized by a series of complex basins that spanned the contemporaneous active continental margin, which itself was characterized by volcanic activity. The basins are filled with thick (up to 3000 m) accumulations of interbedded sedimentary and volcanic strata of late Oligocene-early Miocene age. We interpret that these basins developed during a phase of moderate extension within the plate margin system, triggered by an increased rate of convergence of the Farallon (Nazca) and South American plates between 28 and 26 Ma. This history is inconsistent with models of convergence that link high rates of convergence of a continental margin and an oceanic plate to strong compressional coupling. Although extensional basins of this age are only well-described in the southern Andes, the convergence history and volcanic chronology are similar farther north in the central Andes (18°–33°S), leading to the speculation that extension may have characterized the late Oligocene-early Miocene interval regionally. We hypothesize that this extension was a necessary condition to subsequent building of the modern Andes Mountains.

Heat production in an Archean crustal profile and implications for heat flow and mobilization of heat-producing elements

We have measured concentrations of heat producing elements (Th, U, and K) in 58 samples representative of the main lithologies in a 100 km transect of the Superior Province of the Canadian Shield, from the Michipicoten (Wawa) greenstone belt, near Wawa, Ontario, through a domal gneiss terrane of amphibolite grade, to the granulite belt of the Kapuskasing Structural Zone, near Foleyet. This transect has been interpreted as an oblique cross section through some 25 km of crust, uplifted along a major thrust fault, and thus provides an opportunity to examine in detail a continuous profile into deep continental crust of Archean age. Mean heat production values for these terranes, based on aereal distribution of major rock types and calculated from their Th, U, and K concentrations are: Michipicoten greenstone belt = 0.72 μW m−3; Wawa domal gneiss terrane (amphibolite grade) = 1.37 μW m−3; Kapuskasing granulites = 0.44 μW m−3. Among the silicic plutonic rocks (tonalites, granites, and their derivative gneisses), the relatively large variation in heat production correlates with modal abundances of accessory minerals including allanite, sphene, zircon, and apatite. We interpret these variations as primary (pre-metamorphic). The relatively high weighted mean heat production of the domal gneiss terrane can be accounted for by the larger proportion there of late-stage Th-, U-, and K-rich granitoid plutons. These may have been derived from the underlying Kapuskasing granulite terrane, leaving it slightly depleted in heat producing elements. Transport of Th, U, and K, therefore, could have taken place in silicate melts rather than in aqueous or carbonic metamorphic fluids. This conclusion is supported by the lack of a statistically significant difference in heat production between tonalites, tonalite gneisses and mafic rocks of amphibolite versus granulite grade. The pre-metamorphic radioactivity profile for this crustal section is likely to have been uniformly low, with a mean heat production value less than 1 μW m−3. This result is distinctly different from measured profiles in more silicic terranes, which show decreasing heat production with depth. This implies fundamental differences in crustal radioactivity distributions between granitic and more mafic terranes, and may be an important factor in selective reactivation of lithologically different terranes, possibly resulting in preferential stabilization of basic terranes in the geological record. Our results indicate that a previously determined apparently linear heat flow-heat production relationship for the Kapuskasing area does not relate to the distribution of heat production with depth. Low, but significant heat production, 0.4–0.5 μW m−3, continues to lower crustal depths with no correlation to the depth parameter from the linear relationship. This low heat production may be a minimum average granulite heat production and suggests that, in general, heat flow through the Moho is 8–10 mW m−2 lower than the reduced heat flow calculated from the heat flow-heat production regression.

Mantle-driven dynamic uplift of the Rocky Mountains and Colorado Plateau and its surface response: Toward a unified hypothesis

The correspondence between seismic velocity anomalies in the crust and mantle and the differential incision of the continental-scale Colorado River system suggests that significant mantle-to-surface interactions can take place deep within continental interiors. The Colorado Rocky Mountain region exhibits low-seismic-velocity crust and mantle associated with atypically high (and rough) topography, steep normalized river segments, and areas of greatest differential river incision. Thermochronologic and geologic data show that regional exhumation accelerated starting ca. 6–10 Ma, especially in regions underlain by low-velocity mantle. Integration and synthesis of diverse geologic and geophysical data sets support the provocative hypothesis that Neogene mantle convection has driven long-wavelength surface deformation and tilting over the past 10 Ma. Attendant surface uplift on the order of 500–1000 m may account for ∼25%–50% of the current elevation of the region, with the rest achieved during Laramide and mid-Tertiary uplift episodes. This hypothesis highlights the importance of continued multidisciplinary tests of the nature and magnitude of surface responses to mantle dynamics in intraplate settings.

Controls on architecture of the Late Cretaceous to Cenozoic southern Middle Magdalena Valley Basin, Colombia

The Maastrichtian-Cenozoic southern Middle Magdalena Valley Basin of Colombia contains a unique record of unconformities, strata, and structure, from which we extract the histories of exhumation of the Central Cordillera, to the west, and evolution of the Eastern Cordillera fold-and-thrust belt, to the east. This study integrates field-based analyses of stratigraphy, laboratory analyses of provenance, fission-track thermochronology, vitrinite-reflectance data, volcanic-ash geochronology, and studies of synorogenic geometries and structure displayed in seismic data. A major unconformity, the Middle Magdalena Valley unconformity, formed by eastward migration of Central Cordillera uplift during Late Cretaceous to early Eocene time, which transformed a Maastrichtian marine basin into a Paleocene depositional piedmont area This transformation is recorded by a coarsening-upward sequence of marine shales to alluvial-fan conglomerates, which was partly eroded during further early Eocene propagation of Central Cordillera deformation. Cessation of this phase of uplift led to formation of a pediment surface, the Middle Magdalena Valley unconformity, which was buried by westward-onlapping middle Eocene to Neogene alluvial deposits. Middle Eocene to Neogene facies, paleoflow, and unconformities were controlled by Eastern Cordillera deformation. In the Eastern Cordillera foot-hills, growth strata and thermal history reveal two phases of folding of middle Eocene-Oligocene and late Miocene ages, prior to intense Pliocene-Pleistocene uplift. Two unconformities of early late Miocene and Pliocene-Pleistocene ages occur to the west of the Eastern Cordillera and record flexural tilting related to episodes of Eastern Cordillera loading.

Syntectonic Cenozoic sedimentation in the northern middle Magdalena Valley Basin of Colombia and implications for exhumation of the Northern Andes

Development of the Colombian Middle Magdalena Valley Basin (MMVB) was determined by late Cretaceous-early Eocene uplift of the Central Cordillera to the west, and subsequent transferal of deformation to the Eastern Cordillera to the east. These phases are separated in the tectono-stratigraphic record by a major unconformity, the Middle Magdalena Valley unconformity (MMVU). Paleocene coastal to alluvial facies underneath the MMVU were deposited in a foreland basin coupled to Central Cordillera kilometer-scale uplift. The middle Eocene to Neogene continental strata that onlap the MMVU document transformation of the MMVB into an interior basin due to Eastern Cordillera deformation recorded by growth strata in seismic lines and changing provenance and paleoflow patterns. Exhumation histories of the MMVB bounding ranges are further constrained by apatit-fission-track and vitrinite reflectance thermochronology and volcanic ash chronology. This study indicates that the MMVB is fundamentally related to evolution of the entire Andean margin of South America.

Tectonic controls on the hydrogeology of the Rio Grande Rift, New Mexico

Mathematical modeling is used in this study to assess how tectonic movement of fault blocks and fault permeability influence the present-day and paleohydrogeology of the Rio Grande Rift near Socorro, New Mexico. Our analysis focuses on active and ancient groundwater flow patterns and hot spring development within the southern La Jencia and Socorro subbasins. The best agreement between model results and present-day and paleoheat flow data was achieved by representing faults as passive surfaces and incorporating 2 km of moderately permeable (10−14.0 m2) fractured crystalline rocks into the hydrogeologic model. Quantitative results indicate that changes in groundwater flow patterns across the basin are primarily generated by the truncation/reconnection of aquifers and confining units. Calculated flow patterns help to explain the annealing of apatite fission tracks within Eocene Baca Formation clasts to the east of Socorro, potassium metasomatism mass balance constraints within Oligocene volcanics and overlying Santa Fe Group deposits, and the timing of barite/fluorite ore mineralization within the Gonzales prospect on the eastern edge of the Rio Grande Rift. We estimate that about 5% of mountain front recharge penetrates to a depth of 2.8 km below the sedimentary pile. This may have implications for water resource planners who have historically focused on groundwater resource development within the shallow alluvial deposits along the Rio Grande Rift valley.

Middle Tertiary buoyancy modification and its relationship to rock exhumation, cooling, and subsequent extension at the eastern margin of the Colorado Plateau

In the southern Rocky Mountains and Rio Grande rift, rock cooling patterns from apatite fission-track (AFT) data spatially correlate with areas of voluminous middle Tertiary caldera-complex magmatism. We use thermochronology and gravity data to explore lithospheric modification by voluminous middle Tertiary magmatism. These data are not traditionally used to constrain magmatic processes, but in our study area they provide first-order constraints on the degree of mantle dedensification by basalt removal. We show that thermal isostatic responses to middle Tertiary magmatism drove spatially variable rock uplift and thermal perturbations that, coupled with variable exhumation, can explain AFT cooling and rock preservation patterns. We further argue that, if rock uplift exceeded exhumation, then surface uplift combined with magmatic weakening of the lithosphere could have influenced subsequent Neogene extension.

Diachronous episodes of Cenozoic erosion in southwestern North America and their relationship to surface uplift, paleoclimate, paleodrainage, and paleoaltimetry

The history of erosion of southwestern North America and its relationship to surface uplift is a long-standing topic of debate. We use geologic and thermochronometric data to reconstruct the erosion history of southwestern North America. We infer that erosion events occurred mostly in response to surface uplift by contemporaneous tectonism, and were not long-delayed responses to surface uplift caused by later climate change or drainage reorganization. Rock uplift in response to isostatic compensation of exhumation occurred during each erosion event, but has been quantified only for parts of the late Miocene–Holocene erosion episode. We recognize four episodes of erosion and associated tectonic uplift: (1) the Laramide orogeny (ca. 75–45 Ma), during which individual uplifts were deeply eroded as a result of uplift by thrust faults, but Laramide basins and the Great Plains region remained near sea level, as shown by the lack of significant Laramide exhumation in these areas; (2) late middle Eocene erosion (ca. 42–37 Ma) in Wyoming, Montana, and Colorado, which probably occurred in response to epeirogenic uplift from lithospheric rebound that followed the cessation of Laramide dynamic subsidence; (3) late Oligocene–early Miocene deep erosion (ca. 27–15 Ma) in a broad region of the southern Cordillera (including the southern Colorado Plateau, southern Great Plains, trans-Pecos Texas, and northeastern Mexico), which was uplifted in response to increased mantle buoyancy associated with major concurrent volcanism in the Sierra Madre Occidental of Mexico and in the Southern Rocky Mountains; (4) Late Miocene–Holocene erosion (ca. 6–0 Ma) in a broad area of southwestern North America, with loci of deep erosion in the western Colorado–eastern Utah region and in the western Sierra Madre Occidental. Erosion in western Colorado–eastern Utah reflects mantle-related rock uplift as well as an important isostatic component caused by compensation of deep fluvial erosion in the upper Colorado River drainage following its integration to the Gulf of California. Erosion in the western Sierra Madre Occidental occurred in response to rift-shoulder uplift and the proximity of oceanic base level following the late Miocene opening of the Gulf of California. We cannot estimate the amount of rock or surface uplift associated with each erosion episode, but the maximum depths of exhumation for each were broadly similar (typically ∼1–3 km). Only the most recent erosion episode is temporally correlated with climate change. Paleoaltimetric studies, except for those based on leaf physiognomy, are generally compatible with the uplift chronology we propose here. Physiognomy-based paleo elevation data suggest that near-modern elevations were attained during the Paleogene, but are the only data that uniquely support such interpretations. High Paleogene elevations require a complex late Paleogene–Neogene uplift and subsidence history for the Front Range and western Great Plains of Colorado that is not compatible with the regional sedimentation and erosion events we describe here. Our results suggest that near-modern surface elevations in southwestern North America were generally not attained until the Neogene, and that these high elevations are the cumulative result of four major episodes of Cenozoic rock uplift of diverse origin, geographic distribution, and timing.

New thermochronometric constraints on the Tertiary landscape evolution of the central and eastern Grand Canyon, Arizona

Thermal histories are modeled from new apatite (U-Th)/He and apatite fission-track data in order to quantitatively constrain the landscape evolution of the Grand Canyon region. Fifty new samples and their associated thermochronometric ages are presented here. Samples span from Lee’s Ferry in the east to Quartermaster Canyon in the west and include four age-elevation transects into Grand Canyon and borehole samples from the Coconino Plateau. Twenty-seven samples are inversely modeled to provide continuous thermal histories. This represents the most extensive and complete dataset on patterns of long-term exhumation in the Grand Canyon region, and it enables us to constrain the timing and magnitude of erosion and also discriminate between canyon incision and broader planation. The new data suggest that the early Cenozoic landscape in eastern Grand Canyon was low in relief and does not indicate the presence of an early Cenozoic precursor to the modern Grand Canyon. However, there is evidence for the incision of a smaller-scale canyon across the Kaibab Uplift at 28–20 Ma. This middle-Cenozoic denudation event was accompanied by the removal of a majority of remaining Mesozoic strata west of the Kaibab Uplift. In contrast, just upstream in the area of Lee’s Ferry, ∼2 km of Mesozoic strata remained over the middle Cenozoic and were removed after 10 Ma.

Denudation and flexural isostatic response of the Colorado Plateau and southern Rocky Mountains region since 10 Ma

Over the past 10 Ma, the high-relief landscapes of the Colorado Plateau–southern Rocky Mountains region have been shaped by erosional processes. Incision rates have increased in the southern Rocky Mountains, the Colorado River system has been superimposed across buried Laramide structures as it was integrated from the Rocky Mountains to the Gulf of California, the modern Grand Canyon formed, and there has been widespread denudation of the Canyonlands region of the Colorado Plateau. We examine the spatial and temporal distribution of erosion and its associated isostatic rebound since 10 Ma. Erosion estimates come from apatite fission track (AFT) and apatite (U-Th)/He (AHe) thermochronometric studies at 14 sites across the region, including recent AHe data with ages younger than 12 Ma, and from ca. 10 Ma 40 Ar/ 39 Ar dated basalt paleosurfaces at 55 locations on the perimeter of the Colorado Plateau and in the southern Rocky Mountains. Estimated eroded thickness is added to modern topography above numerous control points to reconstruct a 10 Ma paleosurface across the region (referenced to modern elevations); this also yields an eroded thickness volume. Erosion has been spatially variable since 10 Ma: we find widespread denudation with as much as 2 km of incision along rivers in the Canyonlands region of Utah, 1–1.5 km of incision along rivers exiting the Rocky Mountains onto the eastern piedmont since 6 Ma, ∼1 km removed across the high peaks of the southern Rocky Mountains since 10 Ma, and little net erosion in the Basin and Range. Post–10 Ma flexural isostatic response to the eroded volume is calculated using known variable elastic thickness. This rebound caused much of the Colorado Plateau region to undergo more than 800 m of rock uplift, exceeding 1 km in local areas in the Canyonlands and southwestern Colorado. The Lees Ferry and Glen Canyon areas have been isostatically uplifted >500 m relative to the eastern Grand Canyon and the Tavaputs Plateau has been isostatically uplifted 400 m relative to Browns Park. This differential rock uplift driven by erosional isostasy has created or accentuated many of the features of the modern landscape. This component of rock uplift is “removed” by adding the eroded thickness onto modern topography, then subtracting the calculated rebound. The resulting (pre-erosion and pre-rebound) map provides a model of the 10 Ma landscape, neglecting any tectonic uplift contribution to regional elevations. This model suggests the presence of internal drainages on the Colorado Plateau, that the elevation of the Green River Basin and the Tavaputs Plateau were subequivalent, allowing the Green River to flow southward, and shows high topography in the Rocky Mountains that mimicked modern topography, but with potentially lower relief. Future refinements of both the timing and magnitude of differential erosion and rebound models provide an avenue for improved models for Cenozoic landscape evolution of the region. This paper is an advance over previous studies that focused just on the Colorado Plateau. Here we evaluate isostatic response to erosion in an extended region that includes parts of the Basin and Range, Colorado Plateau, southern Rocky Mountains, and eastern piedmont of the Rocky Mountains. We find that erosion of the southern Rocky Mountains and eastern piedmont is comparable to that of the Colorado Plateau and that the flexural isostatic rebounds of all these regions are coupled and cannot be considered in isolation. Furthermore, we focus on the 10 Ma time frame, rather than the 30 or 70 Ma period of previous researchers, as the key time frame during which the modern landscape rapidly evolved. In addition, the use of AFT and AHe thermochronometric constraints on thicknesses and ages of now-eroded sediments has solved key problems that hampered previous erosion studies. Data and analyses of regional post–10 Ma differential erosion and its resulting differential isostatic rebound provide essential constraints for any viable models for landscape evolution in this classic region.