Search

During the termination of the last glacial period the western U.S. experienced exceptionally wet conditions, driven by changes in location and strength of the mid-latitude winter storm track. The distribution of modern winter precipitation is frequently characterized by a north-south wet/dry dipole pattern, controlled by interaction of the storm track with ocean-atmosphere conditions over the Pacific and Atlantic Oceans. Here we show that a dipole pattern of similar geographic extent persisted and switched sign during millennial-scale abrupt climate changes of the last deglaciation, based on a new lake level reconstruction for pluvial Lake Chewaucan (northwestern U.S.), and a compilation of regional paleoclimate records. This suggests the dipole pattern is robust, and one mode may be favored for centuries, thereby creating persistent contrasting wet/dry conditions across the western U.S. The TraCE-21k climate model simulation shows an equatorward enhancement of winter storm track activity in the northeastern Pacific, favoring wet conditions in southwestern U.S. during the second half of Heinrich Stadial 1 (16.1–14.6 ka) and consistent with paleoclimate evidence. During the Bølling/Allerød (14.6–12.8 ka), the northeastern Pacific storm track contracted poleward, consistent with wetter conditions concentrated poleward toward the northwest U.S.

Isotopic compositions of lacustrine carbonates are commonly used for dating and paleoenvironmental reconstructions. Here we use carbonate δ 13 C and δ 18 O, clumped (Δ 47 ), and 14 C compositions to better understand the carbonate isotope system in closed-basin lakes and trace the paleohydrologic and temperature evolution in the Chewaucan closed-basin lake system, northern Great Basin, USA, over the Last Glacial/Holocene transition. We focus on shorezone tufas to establish that they form in isotopic equilibrium with lake water and DIC, they can be dated reliably using 14 C, and their clumped isotope composition can be used to reconstruct past lake temperature. Calculations of the DIC budget and reservoir age for the lake indicate residence time is short, and dominated by exchange with atmospheric CO 2 at all past lake levels. Modern lake DIC and shorezone tufas yield δ 13 C and 14 C values consistent with isotopic equilibrium with recent fossil fuel and bomb-influenced atmospheric CO 2 , supporting these calculations. δ 13 C values of fossil tufas are also consistent with isotopic equilibrium with pre-industrial atmospheric CO 2 at all shoreline elevations. This indicates that the 14 C reservoir effect for this material is negligible. Clumped isotope (Δ 47 ) results indicate shorezone tufas record mean annual lake temperature. Modern (average 13 ± 2 °C) and 18 ka BP-age tufas (average 6 ± 2 °C) have significantly different temperatures consistent with mean annual temperature lowering of 7 ± 3 °C (1 SE) under full glacial conditions. For shorezone tufas and other lake carbonates, including spring mounds, mollusk shells, and ostracod tests, overall δ 13 C and δ 18 O values co-vary according to the relative contribution of spring and lacustrine end member DIC and water compositions in the drainage system, but specific isotope values depend strongly upon sample context and are not well correlated with past lake depth. This contrasts with the interpretation that carbonate isotopes in closed-basin lake systems reflect changes in DIC and water budgets connected to higher or lower lake volumes. Instead, a small overlapping range of isotope compositions characterize multiple lake levels, so that none can be identified uniquely by isotope composition alone. Relative to other lake carbonates, δ 13 C and δ 18 O values for ostracods in Ana River Canyon deposits are very strongly influenced by Ana River water, suggesting low lake level and volume characterized Summer Lake for most of the past 100,000 years. Coupled with sedimentologic observations, the Ana River deposits thus suggest dry conditions like today are close to the mean climate state in the northern Great Basin. By contrast, basin-integrating highstands such as that dating to ∼14 ka BP, during the last glacial termination, are hydrologically unique and short-lived. Overall, our results indicate carbonate isotope records must account for the specific geochemical and hydrologic characteristics of lake system in order to provide robust paleoenvironmental reconstructions.

The arid northwest Great Basin underwent substantial hydroclimate changes in the past 15,000 years, greatly affecting its desert ecosystems and prehistoric people. There are conflicting interpretations of the timing of hydrologic changes in this region, requiring more records to resolve the dominant climatic drivers. The Paisley Caves archaeological site, located near former pluvial Lake Chewaucan, contains well-dated, stratified sediments best known for evidence of early human occupation in North America. We present a novel paleohydrologic record for the Chewaucan basin based on the frequency of fish remains (Salmonidae and Cypriniformes, likely tui chub) and their carbon, oxygen, and strontium isotope compositions, from the Paisley Caves. Cypriniformes abundance peaks first at the start of the Bolling/Allerod warm interval (∼14.7 ka) and again during the early Younger Dryas (∼12.8 ka). Isotope compositions indicate tui chub were derived from an expansive Lake Chewaucan throughout the Bolling/Allerod, but mainly from spring- or stream-influenced sources by the late Younger Dryas to the present. Fish abundance dropped sharply through the Younger Dryas and early Holocene, when isotope compositions indicate a mix of habitats. Isotope compositions indicate the driest conditions during the middle Holocene, followed by slightly wetter conditions up to the present. This record agrees with recent pluvial lake reconstructions, supporting the hypothesis that a northward shift in the winter storm track supported deep lakes throughout the Bolling/Allerod in the northwest Great Basin. Lake level decline during the Younger Dryas suggests drying climate, differing from more southerly records. During the Holocene, however, shifts in Chewaucan basin hydrology are consistent with the rest of the western U.S. This highlights the need for region-specific records to inform predictions of the hydrologic impact of climate change on arid regions.

The nature of multidecadal to centennial variability of the Asian monsoon is largely unknown, hindering our understanding of how the modern precipitation regime compares to preindustrial variability, as well as how the Asian monsoon may change in coming decades. Here we use the sediment record from a closed-basin lake in southern Tibet, Ngamring Tso, along with other published paleoclimatic records, to assess summer monsoon precipitation variability over the 20th century and during the late Holocene (4100 cal yr BP to present). The first principal component of the Ngamring Tso grain size record, which is highly correlated with median grain size, correlates significantly with observed summer (June–September) precipitation. From CE 1940 to 2007, grain size decreased with increasing summer precipitation and increased with decreasing summer precipitation. Satellite images of Ngamring Tso reveal larger lake area in the 1970s and 2000s and smaller lake area in the 1980s and 1990s, suggesting monsoon precipitation-induced changes in lake depth or area likely govern grain size variability. In the context of the last 4100 yr, prolonged periods of relatively weak summer monsoon precipitation occurred from 2800–2600 cal yr BP, 2500–2300 cal yr BP, and 1600–400 cal yr BP. A trend toward increased summer precipitation began around 1000 cal yr BP, with above-average summer precipitation from 400 cal yr BP to present, peaking between 200–100 cal yr BP (CE 1750–1850). Dry and wet periods are coincident with dry and wet periods in other south-central Tibetan lake sediment records that reflect precipitation or moisture balance, and with regional proxies of the ISM and EASM, indicating south-central Tibet is influenced by both monsoon subsystems. However, western Tibet lake sediment records do not covary with the Ngamring Tso record, suggesting the spatial variability in Tibet precipitation that occurs on interannual timescales today may also characterize precipitation on multidecadal and centennial timescales during the late Holocene. 20th century precipitation variability in southern Tibet falls within the range of natural variability in the last 4100 yr, and does not show a clear trend of increasing precipitation as projected by models. Instead, it appears that poorly understood multidecadal to centennial internal modes of monsoon variability remained influential throughout the last 4100 yr. Thus, substantial multidecadal to centennial-scale hydroclimatic variability will complicate our ability to project and prepare for anthropogenic changes in the region's monsoon precipitation.

The Huancane II moraines deposited by the Quelccaya Ice Cap in southern Peru were selected by the CRONUS-Earth Project as a primary site for evaluating cosmogenic-nuclide scaling methods and for calibrating production rates. The CRONUS-Earth Project is an effort to improve the state of the art for applications of cosmogenic nuclides to earth-surface chronology and processes. The Huancane II moraines are situated in the southern Peruvian Andes at about 4850 m and ∼13.9°S, 70.9°W. They are favorable for cosmogenic-nuclide calibration because of their low-latitude and high-elevation setting, because their age is very well constrained to 12.3 ± 0.1 ka by 34 radiocarbon ages on peat bracketing the moraines, and because boulder coverage by snow or soil is thought to be very unlikely. However, boulder-surface erosion by granular disintegration is observed and a ∼4% correction was applied to measured concentrations to compensate. Samples from 10 boulders were analyzed for 10Be, 26Al, and 36Cl. Interlaboratory bias at the ∼5% level was the largest contributor to variability of the 10Be samples, which were prepared by three laboratories (the other two nuclides were only prepared by one laboratory). Other than this issue, variability for all three nuclides was very low, with standard deviations of the analyses only slightly larger than the analytical uncertainties. The site production rates (corrected for topographic shielding, erosion, and radionuclide decay) at the mean site elevation of 4857 m were 45.5 ± 1.6 atoms 10Be (g quartz)−1 yr−1, 303 ± 15 atoms 26Al (g quartz)−1 yr−1, and 1690 ± 100 atoms 36Cl (g K)−1 yr−1. The nuclide data from this site, along with data from other primary sites, were used to calibrate the production rates of these three nuclides using seven global scaling methods. The traditional Lal formulation and the new Lifton-Sato-Dunai calibrations yield average ages for the Huancane samples that are in excellent-to-good agreement with the radiocarbon age control (within 0.7% for 10Be and 36Cl and 6% for 26Al). However, all of the neutron-monitor-based methods yielded ages that were too young by about 20%. The nuclide production ratios at this site are 6.74 ± 0.34 for 26Al/10Be in quartz and 37.8 ± 2.3 (atoms 36Cl (g K)−1) (atom 10Be (g SiO2)−1)−1 for 36Cl/10Be, in sanidine and quartz, respectively.

The CRONUS-Earth and CRONUS-EU Projects have recently examined the reproducibility of cosmogenic nuclide measurements and the ability of cosmogenic-production models to fit high-quality global calibration data sets. In both cases, although results are adequate for present geochronological needs, they must be improved to meet future demands. We recommend that the cosmogenic-nuclide community embark on a continuing effort to understand the sources of the remaining model-versus-data discrepancy, and that in the meanwhile all cosmogenic ages be benchmarked against reference ages calculated in a consistent fashion. We also recommend that interlaboratory reference materials be routinely analyzed along with samples and blanks, that the results be tracked, and that methods of achieving consistency and precision in sample preparation and analysis be shared within the community. We also recommend additional organization of the community to facilitate the use of intercomparison materials, foster analytical improvements, oversee the evolution of the community age calculator, and synthesize ongoing calibration efforts.

Geological surface-exposure dating using cosmogenic-nuclide accumulation became a practical geochronological endeavor in 1986, when the utility of 10Be, 26Al, 36Cl, and 3He were all demonstrated. In response to the lack of a common basis for quantifying analytical consistency and calibrating cosmogenic-nuclide production, the CRONUS-Earth Project in the U.S. was started in 2005, along with a European partner project, CRONUS-EU. The goal of the CRONUS-Earth Project was to improve the accuracy and precision of terrestrial cosmogenic nuclide dating in general, focusing especially on nuclide production rates and their variation with altitude, latitude, and time, and to attempt to move from empirically based methods to ones with a stronger basis in physics. The CRONUS-Earth Project conducted extensive intercomparisons of reference materials to attempt to quantify analytical reproducibility at the community level. We found that stated analytical uncertainties nearly always underestimate the actual degree of variability, as quantified by the over-all coefficient of variation of the intercalibration data. The average amount by which the actual coefficient of variation exceeded the analytical uncertainty was a factor of two (100%), but ranged from 15% to 300% depending on the nuclide and material. Coefficients of variation ranged from 3−4% for 10Be to 6–8% for 36Cl, 14C, and 21Ne, to 5–11% for 26Al. Both interlaboratory bias and within-laboratory excess spread of the data played a role in increasing variability above the stated analytical uncertainties. The physical basis for cosmogenic nuclide production was investigated through numerical modeling and the measurement of energy-dependent neutron cross sections for nuclide interactions. We formulated new, physically based, scaling models, denoted LSD and LSDn, by generalizing global numerical simulations of cosmic-ray processes. The CRONUS-Earth Project identified new geological calibration sites, including one at low latitude and high elevation (Huancane, Peru), and replicated nuclide measurement at numerous laboratories. At many sites multiple nuclides were measured, providing much more confidence in the equivalence of surface-exposure ages calculated from differing nuclides. The data were interpreted using an original cosmogenic-nuclide calculator, CRONUScalc, that incorporates the new physically based scaling. The new data and model produced significantly better fits than previous efforts, but do not fully resolve apparent spatial variations in production rates. The CRONUS-Earth and CRONUS-EU Projects have provided a firm foundation for assessing the strengths and weaknesses of cosmogenic-nuclide analytical methods, adjusted the AMS standards for 10Be and consequently revised the half-life, and have provided improved calibration data sets and interpretative tools.

Well-dated records of tropical glacier fluctuations are essential for developing hypotheses and testing proposed mechanisms for past climate changes. Since organic material for radiocarbon dating is typically scarce in low-latitude, high-altitude environments, surface exposure-age dating, based on the measurement of in situ produced cosmogenic nuclides, provides much of the chronologic information on tropical glacier moraines. Here, we present a locally calibrated Be-10 production rate for a low-latitude, high-altitude site near Quelccaya Ice Cap (similar to 13.95 degrees S, 70.89 degrees W, 4857 m asl) in the southeastern Peruvian Andes. Using an independent age (12.35 +0.2, -0.02 ka) of the late glacial Huancane Ha moraines based on thirty-four bracketing radiocarbon ages and twelve Be-10 concentrations of boulders on the moraines, we determine a local production rate of 43.28 +/- 2.69 atoms gram(-1) year (at g(-1) yr(-1)). Reference Be-10 production rates (i.e., production rates by neutron spallation appropriate for sea-level, high-latitude sites) range from 3.97 +/- 0.09 to 3.78 +/- 0.09 at g(-1) yr(-1), determined using scaling after Lal (1991) and Stone (2000) and depending on our assumed boulder surface erosion rate. Since our boulder surface erosion rate estimate is a minimum value, these reference production rates are also minimum values. A secondary control site on the Huancane Mb moraines suggests that the Be-10 production rates are at least as low as, or possibly lower than, those derived from the Huancane Ha moraines. These sea-level, high-latitude production rates are at least 11-15% lower than values derived using the traditional global calibration dataset, and they are also lower than those derived from the late glacial Breque moraine in the Cordillera Blanca of Peru. However, our sea-level, high-latitude production rates agree well with recently published, locally calibrated production rates from the Arctic, New Zealand, and Patagonia. The production rates presented here should be used to calculate Be-10 exposure ages in low-latitude, high-altitude locations, particularly in the tropical Andes, and should improve the ability to compare the results of studies using Be-10 exposure-age dating with other chronological data. (C) 2013 Elsevier B.V. All rights reserved.

14C dating of shoreline deposits of closed-basin lake Baqan Tso in the western Tibetan Plateau shows that lake level regressed from the undated highstand (46 m above modern, 4.3 × modern surface area) of likely earliest Holocene age by 11.5 ka, and remained larger than modern until at least ≈ 5.0 ka. The shoreline record broadly matches other regional climate records, with lake level closely following Northern Hemisphere summer insolation overprinted by sub-millennial lake-level oscillations. A model coupling modern land runoff and lake surface heat closely reproduces estimated modern precipitation of ≈ 240 mm/yr. We estimate that the Baqan Tso basin required ≈ 380 mm/yr precipitation to sustain the maximum early Holocene lake area, a 55% increase over modern. Precipitation increases, not glacial meltwater, drove lake-level changes, as Baqan Tso basin was not glaciated during the Holocene. Our estimate assumes early Holocene insolation (≈ 1.3% overall increase), and mean annual increases of 2°C in temperature, and 37% in relative humidity. We additionally developed a Holocene precipitation history for Baqan Tso using dated paleolake areas. Using the modern and early Holocene model results as end-members, we estimate precipitation in the western Tibetan Plateau which was 300–380 mm/yr between 5.0 and 11.5 ka, with error of ± 29–57 mm/yr (± 12–15%).

We present a shoreline-based, millennial-scale record of lake-level changes spanning 12.8–2.3 ka for a large closed-basin lake system on the southwestern Tibetan Plateau. Fifty-three radiocarbon and eight U–Th series ages of tufa and beach cement provide age control on paleoshorelines ringing the basin, supplemented by nineteen ages from shell and aquatic plant material from natural exposures generally recording lake regressions. Our results show that paleo-Ngangla Ring Tso exceeded modern lake level (4727 m asl) continuously between ~ 12.8 and 2.3 ka. The lake was at its highstand 135 m (4862 m asl) above the modern lake from 10.3 ka to 8.6 ka. This is similar to other closed-basin lakes in western Tibet, and coincides with peak Northern Hemisphere summer insolation and peak Indian Summer Monsoon intensity. The lake experienced a series of millennial-scale oscillations centered on 11.5, 10.8, 8.3, 5.9 and 3.6 ka, consistent with weak monsoon events in proxy records of the Indian Summer Monsoon. It is unclear whether these events were forced by North Atlantic or Indian Ocean conditions, but based on the abrupt lake-level regressions recorded for Ngangla Ring Tso, they resulted in significant periodic reductions in rainfall over the western Tibetan Plateau throughout the Holocene.

The timing and pattern of Tibetan Pla- teau rise provide a critical test of possible mechanisms for the development and sup- port of high topography, yet views range widely on the history of surface uplift to modern elevations of ~4.5 km. To address this issue we present clumped isotope ther- mometry data from two well-studied ba- sins in central and southwestern Tibet, for which previous carbonate δ 18 O data have been used to reconstruct high paleoeleva- tions from late Oligocene to Pliocene time. Clumped isotope thermometry uses mea- surements of the 13 C- 18 O bond ordering in carbonates to constrain the temperature (T(Δ47)) and δ 18 O value of the water from which the carbonate grew. These data can be used to infer paleoelevation by exploiting the systematic decrease of surface tempera- ture and the δ 18 O-based paleoaltimetry and with paleontological and isotopic data indicating the presence of cold- adapted mammals living in a cold, high-ele- vation climate. We suggest that late Neogene elevation loss across the Zhada Basin catch- ment probably related to local expression of east-west extension across much of the south- ern Tibetan Plateau at this time.

The Rio Grande is the fourth longest river in North America extending over 3,000 km from the Rocky Mountains to the Gulf of Mexico. The Pleistocene evolution of this river from individual subbasins into a coalesced fluvial system has been long debated. Herein, we constrain the middle Pleistocene evolution of the northernmost and largest Rio Grande basin, the San Luis basin, and the timing of incision of the Rio Grande Gorge, based on new geologic mapping, 3He surface exposure dating, and U-series dating on pedogenic carbonates. 3He dating of shoreline deposits of closed basin Lake Alamosa and fluvially scoured bedrock along the Rio Grande gorge rim indicate the San Luis basin was connected with southerly Rio Grande basins and progressively incised since ∼400 ka. Waning tectonic activity, coupled with the hydrologic response to intensified middle Pleistocene glaciation, drove southward spillover of Lake Alamosa ∼400 ka, expanding the river's drainage basin by nearly 24,000 km2, and adding recharge areas in the high-altitude, glaciated San Juan, Tusas, and Sangre de Cristo Mountains. Geologic barriers to streamflow on the Taos Plateau were breached and basin tributary drainages were integrated and inset into their present canyons ∼250-200 ka, with approximately 120 m of incision during ∼400-200 ka. Previous work demonstrates that southerly streamflow through basins south of the SLB terminated into a large, rift-related bolson complex until the middle Pleistocene, and most axial Rio Grande river incision has occurred since 630 ka. We propose that middle Pleistocene integration of the SLB was the major geomorphic event that eroded Plio-Pleistocene tectonic and volcanic barriers and re-established flow of the Rio Grande from the Rocky Mountains to the Gulf of Mexico.

Microlithic artifacts, some found in situ, are abundant in the Zhongba archaeological site in southwestern Tibet. The site environment consists of extant wetlands and paleo-wetland deposits found in depressions between sand dunes derived from the Yarlung Tsangpo floodplain. Constraining 14C dates from wetland vegetation and shell from one site fall between ca. 6600–2600 cal. yr B.P., while a second site is dated 3400–1200 cal. yr B.P. A significant and variable 14C reservoir effect—up to 1400 14C yr—limits these ranges to terminus post quem constraints. The in situ artifacts are supplemented by surface collections fully characterizing raw material and typological variability for each site. Raw material found at Zhongba is chert and chalcedony likely sourced from Cretaceous bedrock near the site. Typologically, microblades are nongeometric and are derived from conical and wedge-shaped cores similar to those identified in the Qinghai Lake Basin and the Chang Tang Nature Reserve of similar or greater age. The later occupation period at Zhongba is broadly contemporaneous with sites on the Qinghai-Tibet Plateau containing bronze and iron artifacts, indicating microlithic technology remained an important tool-making strategy in western Tibet late into the protohistoric period.

International audience; Paleoaltimetric studies have characterized in detail the relationship between carbonate oxygen isotope ratios (d18Oc) and elevation in orogens with simple, single-moisture-source hydrological systems, and applied this relationship to ancient continental carbonates to provide constraints on their past elevation. However, mixing of different atmospheric moisture sources in low-elevation orogens should affect d18Oc values, but this effect has not yet been confirmed unequivocally. In the American Southwest, summer monsoonal moisture, sourced in the Equatorial Pacific and the Gulf of Mexico, and winter moisture, sourced in the East Pacific, both contribute to annual rainfall. We present stable isotope results from Quaternary carbonates within the American Southwest to characterize the regional d18Oc-elevation relationship. We then provide stable isotope results from local Eocene carbonates to reconstruct late Laramide paleoelevations. The Quaternary d18Oc-elevation relationship in the American Southwest is not as straightforward as in more simple hydrological systems. d18Oc changes with altitude are non-linear, scattered, and display an apparent isotopic lapse rate inversion above 1200 m of elevation. We speculate that decreasing surface temperatures at high altitudes limit the duration of carbonate growth to the summer months, biasing d18Oc values toward higher values typical of the summer monsoon and leading to lapse rate inversion. d18Oc-elevation relationships based on modern water isotope data or distillation models predict paleoelevations that range up to as much as 2 km higher than the modern elevations of 2000 to 2400 m for our late Eocene sites located at the southern edge of the Colorado Plateau. By contrast, our d18Oc-elevation relationship for the American Southwest yields lower paleoelevation estimates. These alternate estimates nonetheless suggest that significant elevation (at least ϳ1 km) had already been attained by the Eocene, but are also compatible with < 1 km of uplift by post-Laramide mechanisms. Our results show the limitations of standard d18Oc-elevation models in complex hydrological systems and suggest that similar mechanisms may have led to summer-biased paleoaltimetry estimates for the initial stages of other orogenies -in the American Southwest and elsewhere.

The last glacial-interglacial transition (ca. 18–11 ka) was interrupted by abrupt climate events that differed in each hemisphere. During the Antarctic Cold Reversal (ca. 14.5–12.9 ka), the Southern Hemisphere high and mid latitudes cooled, while the Northern Hemisphere warmed. The pattern of change then reversed during the Younger Dryas (ca. 12.9–11.7 ka), which was characterized by cold conditions in much of the Northern Hemisphere. Well-dated paleoclimate records serve to reveal the possible mechanisms for these events. Here we present a reconstruction of the late glacial fluctuations of Quelccaya Ice Cap, located in the southern tropics, based on 38 new radiocarbon ages. Quelccaya was retreating from its Last Glacial Maximum (ca. 21 ka) extent by ca. 17.2 ka, and was located upvalley from its late glacial moraines by 13.6–12.8 ka. Quelccaya experienced a significant readvance that culminated at 12.5–12.4 ka, and then receded several kilometers to near, or within, its late Holocene extent by ca. 11.6 ka. This record provides the most detailed evidence yet of glacier fluctuations in the southern tropics during late glacial time.

Fluvial channels in metamorphic core complexes are preferentially oriented parallel and perpendicular to the direction of tectonic extension. This pattern has been variably attributed to such causes as tectonic tilting during extension, channel elongation by slip along the range-bounding detachment fault, and the exploitation of extension-related joint sets during channel incision. In this paper we use field measurements, digital elevation model analyses, and numerical modeling to test hypotheses for the tectonic and structural control of fluvial channels in metamorphic core complexes, using the Catalina-Rincon core complex in southern Arizona, USA, as a type example. Field measurements and aerial photographic analyses indicate that channels of all sizes exploit steeply dipping joint sets during fluvial incision. As a consequence, channels become preferentially aligned along those joint sets. First and second Strahler-order channels preferentially exploit a joint set oriented perpendicular to the extension direction, while higher-order channels preferentially exploit a joint set oriented parallel to the extension direction. While these observations support the joint-exploitation hypothesis for structural control of drainage architecture, numerical modeling indicates that the spatial distribution of rock uplift during the initial phase of extension plays a crucial role by determining which joint set is preferentially exploited by channels of which Strahler orders. Numerical models indicate that higher-order channels exploit the joint set that is most closely aligned with the direction of initial tectonic tilting, even if that tilting is active for only a short period of time following the initiation of uplift. We conclude that the drainage architecture in the Catalina-Rincon core complex is the result of a combination of joint exploitation and tectonic tilting mechanisms. Structure also plays an important role in controlling the longitudinal profiles of channels in metamorphic core complexes. Channels in the Catalina-Rincon core complex are characterized by structurally controlled knickpoints with a wide distribution of heights and spacings. Field observations indicate that the occurrence of structurally controlled knickpoints and the resulting variability in longitudinal profile form is related to spatial variations in joint density. Numerical models that incorporate spatial variations in joint density using a stochastic bedrock erodibility coefficient are capable of reproducing the statistical properties of longitudinal profiles in the Catalina-Rincon core complex, including the power spectrum of longitudinal profiles and the frequency size distribution of structurally controlled knickpoints. The results of this study illustrate the important roles played by both jointing and the spatial distribution of rock uplift on the geomorphic evolution of metamorphic core complexes. More broadly, the study provides a recipe for how to incorporate joint-related structural controls into landscape evolution models.