From Pliocene to middle Pleistocene time, a large lake occupied most of the San Luis Valley above 2300 m elevation (7550 ft) in southern Colorado. This ancient lake accumulated sediments of the Alamosa Formation (Siebenthal, 1910), for which the lake is herein named. The existence of this lake was first postulated in 1822 and proven in 1910 from well logs. At its maximum extent of nearly 4000 km 2 , it was one of the largest high-altitude lakes in North America, similar to but larger than Lake Texcoco in the Valley of Mexico. Lake Alamosa persisted for ~3 m.y., expanding and contracting and filling the valley with sediment until ca. 430 ka, when it overtopped a low sill and cut a deep gorge through Oligocene volcanic rocks in the San Luis Hills and drained to the south. As the lake drained, nearly 100 km 3 (81 × 10 6 acre-ft or more) of water coursed southward and flowed into the Rio Grande, entering at what is now the mouth of the Red River. The key to this new interpretation is the discovery of ancient shoreline deposits, including spits, barrier bars, and lagoon deposits nestled among bays and in backwater positions on the northern margin of the San Luis Hills, southeast of Alamosa, Colorado. Alluvial and lacustrine sediment nearly filled the basin prior to the lake's overflow, which occurred ca. 430 ka as estimated from 3 He surface-exposure ages of 431 ± 6 ka and 439 ± 6 ka on a shoreline basalt boulder, and from strongly developed relict calcic soils on barrier bars and spits at 2330–2340 m (7645–7676 ft), which is the lake's highest shoreline elevation. Overtopping of the lake's hydrologic sill was probably driven by high lake levels at the close of marine oxygen-isotope stage (OIS) 12 (452–427 ka), one of the most extensive middle Pleistocene glacial episodes on the North American continent. Hydrologic modeling of stream inflow during full-glacial-maximum conditions suggests that Lake Alamosa could fill at modern precipitation amounts if the mean annual temperature were just 5 °C (10 °F) cooler, or could fill at modern temperatures with 1.5 times current mean annual precipitation. Thus, during pluvial epochs the lake would rise to successively higher levels owing to sedimentation; finally during OIS 12, the lake overflowed and spilled to the south. The integration of the upper (Colorado) and lower (New Mexico) reaches of the Rio Grande expanded the river's drainage basin by nearly 18,000 km 2 and added recharge areas in the high-altitude, glaciated San Juan Mountains, southern Sawatch Range, and northern Sangre de Cristo Mountains. This large increase in mountainous drainage influenced the river's dynamics downstream in New Mexico through down-cutting and lowering of water tables in the southern part of the San Luis Valley.

Volcanic clasts incorporated in the lower portion of the Tertiary Santa Fe Group sedimentary rocks of the Culebra graben, San Luis Basin, Colorado, provide constraints on the timing of regional tectonic events by provenance determination. Based on currently exposed volcanic terrains, possible clast sources include Spanish Peaks and Mount Mestas to the east, the San Juan volcanic field to the west, and the Thirtynine Mile volcanic field, a remnant of the Central Colorado volcanic field, to the north and east of the San Luis Basin. Provenance was determined by a variety of geochemical, mineral chemical, and geochronologic data. Large porphyritic Santa Fe Group volcanic clasts are potassic with a wide compositional range from potassic trachybasalt to rhyolite. The whole-rock chemistry of the Culebra graben clasts is similar to that of the Thirtynine Mile and San Juan volcanic fields. Culebra graben amphibole and biotite chemistry is generally consistent with that of rocks of the San Juan volcanic field, but not with Spanish Peaks samples. Trace-element data of Culebra graben volcanic clasts overlap with those of the San Juan and Thirtynine Mile volcanic fields, but differ from those of the Mount Mestas. Thermobarometric calculations using mineral chemistry suggest that many Culebra graben rocks underwent a three-stage crystallization history: ~1120 °C at 7–10 kbar, ~1100 °C at 2.3–4.6 kbar, and hornblende formation ~800 °C at 3 kbar. Within the Culebra graben clasts, zircon rim U-Pb geochronologic systematics as well as amphibole and biotite 40 Ar/ 39 Ar plateau data yield ages ranging from 36 to 29 Ma. These ages are consistent with ages of the Thirtynine Mile volcanic field (36–27 Ma) and the Conejos Formation of the San Juan volcanic field (35–29 Ma), but predate Spanish Peaks (ca. 27–21 Ma) and Mount Mestas (ca. 25 Ma). Based on these data, Spanish Peaks and Mount Mestas are excluded as potential source areas for the Santa Fe Group volcanic clasts in the Culebra graben. The San Juan volcanic field is also an unlikely source due to the distance from the depositional site, the inconsistent paleo-current directions, and the pressure-temperature conditions of the rocks. The most likely scenario is that the Central Colorado volcanic field originally extended proximal to the current location of the Culebra graben and local delivery of volcanic clasts was from the north and northeast prior to the uplift of the Culebra Range and Sangre de Cristo Mountains.

The Sunshine Valley–Costilla Plain, a structural subbasin of the greater San Luis Basin of the northern Rio Grande rift, is bounded to the north and south by the San Luis Hills and the Red River fault zone, respectively. Surficial mapping, neotectonic investigations, geochronology, and geophysics demonstrate that the structural, volcanic, and geomorphic evolution of the basin involves the intermingling of climatic cycles and spatially and temporally varying tectonic activity of the Rio Grande rift system. Tectonic activity has transferred between range-bounding and intrabasin faults creating relict landforms of higher tectonic-activity rates along the mountain-piedmont junction. Pliocene–Pleistocene average long-term slip rates along the southern Sangre de Cristo fault zone range between 0.1 and 0.2 mm/year with late Pleistocene slip rates approximately half (0.06 mm/year) of the longer Quaternary slip rate. During the late Pleistocene, climatic influences have been dominant over tectonic influences on mountain-front geomorphic processes. Geomorphic evidence suggests that this once-closed subbasin was integrated into the Rio Grande prior to the integration of the once-closed northern San Luis Basin, north of the San Luis Hills, Colorado; however, deep canyon incision, north of the Red River and south of the San Luis Hills, initiated relatively coeval to the integration of the northern San Luis Basin. Long-term projections of slip rates applied to a 1.6 km basin depth defined from geophysical modeling suggests that rifting initiated within this subbasin between 20 and 10 Ma. Geologic mapping and geophysical interpretations reveal a complex network of northwest-, northeast-, and north-south–trending faults. Northwest- and northeast-trending faults show dual polarity and are crosscut by north-south– trending faults. This structural model possibly provides an analog for how some intracontinental rift structures evolve through time.

Interpretation of gravity, aeromagnetic, and magnetotelluric (MT) data reveals patterns of rifting, rift-sediment thicknesses, distribution of pre-rift volcanic and sedimentary rocks, and distribution of syn-rift volcanic rocks in the central San Luis Basin, one of the northernmost major basins that make up the Rio Grande rift. Rift-sediment thicknesses for the central San Luis Basin determined from a three-dimensional gravity inversion indicate that syn-rift Santa Fe Group sediments have a maximum thickness of ~2 km in the Sanchez graben near the eastern margin of the basin along the central Sangre de Cristo fault zone, and reach nearly 1 km within the Monte Vista graben near the western basin margin along the San Juan Mountains. In between, Santa Fe Group thickness is negligible under the San Luis Hills and estimated to reach ~1.1 km under the Costilla Plains (although no independent thickness constraints exist, and a range of thicknesses of 600 m to 2 km is geophysically reasonable). From combined geophysical and geologic considerations, pre-rift, dominantly sedimentary rocks appear to increase in thickness from none in the Sanchez graben on the east to perhaps 800 m under the San Luis Hills on the west. The pre-rift rocks are most likely early Tertiary in age, but the presence of Mesozoic and Paleozoic sedimentary rocks cannot be ruled out. Geophysical data provide new evidence that an isolated exposure of Proterozoic rocks on San Pedro Mesa is rooted in the Precambrian basement. This narrow, north-south–trending basement high has ~2 km of positive relief with respect to the base of the Sanchez graben, and separates the graben from the structural depression beneath the Costilla Plains. A structural high composed of pre-rift rocks, long inferred to extend from under the San Luis Hills to the Taos Plateau, is confirmed and found to be denser than previously believed, with little or no overlying Santa Fe Group sediments. Major faults in the study area are delineated by geophysical data and models; these faults include significant vertical offsets (≥1 km) of Precambrian rocks along the central and southern zones of the Sangre de Cristo fault system. Other faults with similarly large offsets of the Santa Fe Group include a fault bounding the western margin of San Pedro Mesa, and other faults that bound the Monte Vista graben in an area previously assumed to be a simple hinge zone at the western edge of the San Luis Basin. A major north-south–trending structure with expression in gravity and MT data occurs at the boundary between the Costilla Plains and the San Luis Hills structural high. Although it has been interpreted as a down-to-the-east normal fault or fault zone, our modeling suggests that it also is likely related to pre-rift tectonics. Aeromagnetic anomalies over much of the area are interpreted to mainly reflect variations of remanent magnetic polarity and burial depth of the 5.3–3.7 Ma Servilleta Basalt of the Taos Plateau volcanic field. Magnetic-source depth estimates are interpreted to indicate patterns of subsidence following eruption of the basalt, with maximum subsidence in the Sanchez graben.

Two Oligocene conglomeratic units, one primarily nonvolcaniclastic and the other volcaniclastic, are preserved on the west side of the Jemez Mountains beneath the 14 Ma to 40 ka lavas and tuffs of the Jemez Mountains volcanic field. Thickness changes in these conglomeratic units across major normal fault zones, particularly in the southwestern Jemez Mountains, suggest that the western margin of the Rio Grande rift was active in this area during Oligocene time. Furthermore, soft-sediment deformation and stratal thickening in the overlying Abiquiu Formation adjacent to the western boundary faults are indicative of syndepositional normal-fault activity during late Oligocene–early Miocene time. The primarily nonvolcaniclastic Oligocene conglomerate, which was derived from erosion of Proterozoic basement-cored Laramide highlands, is exposed in the northwestern Jemez Mountains, southern Tusas Mountains, and northern Sierra Nacimiento. This conglomerate, formerly called, in part, the lower member of the Abiquiu Formation, is herein assigned to the Ritito Conglomerate in the Jemez Mountains and Sierra Nacimiento. The clast content of the Ritito Conglomerate varies systematically from northeast to southwest, ranging from Proterozoic basement clasts with a few Cenozoic volcanic pebbles, to purely Proterozoic clasts, to a mix of Proterozoic basement and Paleozoic limestone clasts. Paleocurrent directions indicate flow mainly to the south. A stratigraphically equivalent volcaniclastic conglomerate is present along the Jemez fault zone in the southwestern Jemez Mountains. Here, thickness variations, paleocurrent indicators, and grain-size trends suggest north-directed flow, opposite that of the Ritito Conglomerate, implying the existence of a previously unrecognized Oligocene volcanic center buried beneath the northern Albuquerque Basin. We propose the name Gilman Conglomerate for this deposit. The distinct clast composition and restricted geographic nature of each conglomerate suggests the presence of two separate fluvial systems, one flowing south and the other flowing north, separated by a west-striking topographic barrier in the vicinity of Fenton Hill and the East Fork Jemez River in the western Jemez Mountains during Oligocene time. In contrast, the Upper Oligocene–Lower Miocene Abiquiu Formation overtopped this barrier and was deposited as far south as the southern Jemez Mountains. The Abiquiu Formation, which is derived mainly from the Latir volcanic field, commonly contains clasts of dacite lava and Amalia Tuff in the northern and southeastern Jemez Mountains, but conglomerates are rare in the southwestern Jemez Mountains.

Geologic mapping, age determinations, and geochemistry of rocks exposed in the Abiquiu area of the Abiquiu embayment of the Rio Grande rift, north-central New Mexico, provide data to determine fault-slip and incision rates. Vertical-slip rates for faults in the area range from 16 m/m.y. to 42 m/m.y., and generally appear to decrease from the eastern edge of the Colorado Plateau to the Abiquiu embayment. Incision rates calculated for the period ca. 10 to ca. 3 Ma indicate rapid incision with rates that range from 139 m/m.y. on the eastern edge of the Colorado Plateau to 41 m/m.y. on the western part of the Abiquiu embayment. The Abiquiu area is located along the margin of the Colorado Plateau–Rio Grande rift and lies within the Abiquiu embayment, a shallow, early extensional basin of the Rio Grande rift. Cenozoic rocks include the Eocene El Rito Formation, Oligocene Ritito Conglomerate, Oligocene–Miocene Abiquiu Formation, and Miocene Chama–El Rito and Ojo Caliente Sandstone Members of the Tesuque Formation (Santa Fe Group). Volcanic rocks include the Lobato Basalt (Miocene; ca. 15–8 Ma), El Alto Basalt (Pliocene; ca. 3 Ma), and dacite of the Tschicoma Formation (Pliocene; ca. 2 Ma). Quaternary deposits consist of inset axial and side-stream deposits of the ancestral Rio Chama (Pleistocene in age), landslide and pediment alluvium and colluvium, and Holocene main and side-stream channel and floodplain deposits of the modern Rio Chama. The predominant faults are Tertiary normal high-angle faults that displace rocks basinward. A low-angle fault, referred to as the Abiquiu fault, locally separates an upper plate composed of the transitional zone of the Ojo Caliente Sandstone and Chama–El Rito Members from a lower plate consisting of the Abiquiu Formation or the Ritito Conglomerate. The upper plate is distended into blocks that range from about 0.1 km to 3.5 km long that may represent a larger sheet that has been broken up and partly eroded. Geochronology ( 40 Ar/ 39 Ar) from fifteen volcanic and intrusive rocks resolves discrete volcanic episodes in the Abiquiu area: (1) emplacement of Early and Late Miocene basaltic dikes at 20 Ma and ca. 10 Ma; (2) extensive Late Miocene–age lava flows at 9.5 Ma, 7.9 Ma, and 5.6 Ma; and (3) extensive basaltic eruptions during the early Pliocene at 2.9 Ma and 2.4 Ma. Clasts of biotite- and hornblende-rich trachyandesites and trachydacites from the base of the Abiquiu Formation are dated at ca. 27 Ma, possibly derived from the Latir volcanic field. The most-mafic magmas are interpreted to be generated from a similar lithospheric mantle during rifting, but variations in composition are correlated with partial melting at different depths, which is correlated with thinning of the crust due to extensional processes.

The late Cenozoic extension in the Rio Grande rift of north-central New Mexico was predominantly accommodated by the north-south–trending Pajarito and Sangre de Cristo normal faults and the intervening east-northeast–striking predominantly strike-slip Embudo fault. Using this segment of the rift as our primary example, we have analyzed a series of three-dimensional nonlinear elastic-plastic finite-element models to assess the role of mechanical interactions between pairs of en echelon rift-scale listric normal faults in the evolution of intervening relay zones. The model results demonstrate that under orthogonal extension and an overall plane-strain deformation, relay zones may evolve in a three-dimensional strain field and along non-coaxial strain paths. The extent of non-plane strain and non-coaxial deformation depends on the fault overlap to spacing ratio, the relative orientations of the bounding faults, and the structural position within the relay zone. The model-derived minimum compressive stress vectors within the relay zone are oblique to the regional extension direction throughout the deformation. Within the Rio Grande rift of north-central New Mexico, the occurrence of northerly striking Neogene faults suggestive of east-west extension in the Española and the San Luis Basins, geographic variations in the vertical-axis rotations from paleomagnetic studies, and secondary fault patterns are consistent with the near-surface variations in the strain field predicted by the model. The model suggests that interaction between the Pajarito and the Sangre de Cristo faults may have played a major role in the evolution of this segment of the rift.

We describe the structure of the eastern Española Basin and use stratigraphic and stratal attitude data to interpret its tectonic development. This area consists of a west-dipping half graben in the northern Rio Grande rift that includes several intrabasinal grabens, faults, and folds. The Embudo–Santa Clara–Pajarito fault system, a collection of northeast- and north-striking faults in the center of the Española Basin, defines the western boundary of the half graben and was active throughout rifting. Throw rates near the middle of the fault system (i.e., the Santa Clara and north Pajarito faults) and associated hanging-wall tilt rates progressively increased during the middle Miocene. East of Española, hanging-wall tilt rates decreased after 10–12 Ma, coinciding with increased throw rates on the Cañada del Almagre fault. This fault may have temporarily shunted slip from the north Pajarito fault during ca. 8–11 Ma, resulting in lower strain rates on the Santa Clara fault. East of the Embudo–Santa Clara–Pajarito fault system, deformation of the southern Barrancos monocline and the Cañada Ancha graben peaked during the early–middle Miocene and effectively ceased by the late Pliocene. The north-striking Gabeldon faulted monocline lies at the base of the Sangre de Cristo Mountains, where stratal dip relations indicate late Oligocene and Miocene tilting. Shifting of strain toward the Embudo–Santa Clara–Pajarito fault system culminated during the late Pliocene–Quaternary. Collectively, our data suggest that extensional tectonism in the eastern Española Basin increased in the early Miocene and probably peaked between 14–15 Ma and 9–10 Ma, preceding and partly accompanying major volcanism, and decreased in the Plio-Pleistocene.

New 40 Ar/ 39 Ar results from drill-hole cuttings of basaltic and basaltic andesite flows from the Guaje well field of the Pajarito Plateau along the western part of the Española Basin in north-central New Mexico yielded Middle Miocene ages (11.5–13.2 Ma). The volcanic eruptions were closely associated with intense faulting, subsidence, and sedimentation, and the results provide age constraints for the volcanic and tectonic processes along the western margin of the Española Basin. The Middle Miocene volcanic rocks are interbedded within the Santa Fe Group, which is divided into the Hernandez and Vallito Members of the Chamita Formation and the Chama–El Rito Member of the Tesuque Formation, in descending stratigraphic order. New and published geochemical results from the Guaje well field and from other surface and subsurface mafic and intermediate lava flows within the Pajarito Plateau suggest that the volcanic rocks erupted from different magmatic sources and centers close to the Pajarito fault zone. Multiple pulses of volcanic eruptions mostly confined to the hanging wall of the Pajarito fault zone, which represents the current western boundary fault of the Española Basin, suggest that the Pajarito fault system has been sporadically reactivated several times, beginning at least in the Middle Miocene and continuing to the Plio-Pleistocene. Moreover, the volcanic, tectonic, and sedimentary records in the Pajarito Plateau suggest that there is no evidence for eastward migration of tectonic and volcanic activities from the Cañada de Cochiti fault zone in the southern part of the Jemez Mountains to the Pajarito fault zone during the early Pliocene (4–5 Ma).

In the Galisteo drainage basin south of Santa Fe, a fold and several faults related to the Rio Grande rift deform late Eocene–Oligocene dikes, laccoliths, and the Espinaso Formation. The largest rift-related feature, a northerly plunging syncline, comprises the south end of the Santa Fe embayment of the Española Basin and the northern end of the Estancia Basin. The east limb of the syncline is cut by northerly trending, graben-forming, normal faults of the Agua Fria fault system in the Santa Fe embayment. East of the Tijeras-Cañoncito fault system, the east limb of the Estancia Basin is disrupted by down-to-the-west, normal faults of the Glorieta Mesa boundary fault and the Hub Mesa fault system. The fold is offset by down-to-the-northwest movement, and a small component of right slip, on the Tijeras-Cañoncito fault system, which separates the two basins. The above-mentioned rift-related fold and north-trending faults are superimposed across the southeastern margin of the San Luis uplift and the younger Galisteo Basin. Geologic maps and drill data reveal four, and possibly five, phases of Laramide deformation associated with recurrent movement along the Tijeras-Cañoncito fault system: (1) a possible Late Cretaceous, cryptic phase of strike slip associated with elevation of the highest portions of the Santa Fe Range uplift to the north-northeast; (2) the early Paleocene San Luis uplift that formed a southwest plunging, V-shaped anticlinal nose whose southeast limb is the Lamy monocline, which extends 25 km southwest from Precambrian basement at the margin of the Santa Fe Range at Cañoncito to the Cretaceous Menefee Formation; (3) following erosional beveling, the collapse of the southern shoulder of the San Luis uplift, forming a portion of the north-northeast–trending, latest Paleocene–Eocene Diamond Tail subbasin, the axial portion of which lies along the trend of the Tijeras-Cañoncito fault zone; (4) minor Eocene uplift which interrupted deposition in the basin; and (5) Eocene subsidence across the broader Galisteo subbasin and deposition of the Galisteo Formation and latest Eocene–Oligocene Espinaso Formation. Late Eocene–Oligocene intrusions in Los Cerrillos and the Ortiz Mountains deformed the Cretaceous and Tertiary host rocks. Across the Tijeras-Cañoncito fault system, the northwest-trending erosional edge of the Campanian Point Lookout Sandstone displays 400 m of pre–Diamond Tail Formation, right-lateral separation, and the Diamond Tail Formation shows no lateral offset between the overlapping San Lazarus and Los Angeles faults. Although the axis of the Galisteo Basin parallels the fault system, and the basin has been proposed to have formed in a releasing bend of a strike-slip fault along the Tijeras-Cañoncito fault system, any major Laramide strike-slip movement pre-dates the deposition of the Diamond Tail Formation and the formation of the Lamy monocline. The faulted core of the pre–Diamond Tail Lamy monocline, initially up ~800 m on the northwest, was reactivated during rift development and downdropped on the northwest by ~250 m of dip slip. An earlier period of movement (either early Laramide or older strike slip or down-to-the southeast Pennsylvanian movement) is suggested by contrasting thicknesses of Paleozoic formations across the fault zone.

We have developed a conceptual model for the Tesuque aquifer system in the southeastern Española Basin near Santa Fe, New Mexico, based on measurements of chemical, isotopic, and thermal properties of groundwater from 120 wells. This study concentrates on a single groundwater-flow unit (GFU) of the Tesuque aquifer associated with the Santa Fe River drainage, where groundwater flows east to west across north-trending rift structures. We examine links between groundwater flow, temperature, water chemistry, and major fault structures. Hydrologic and hydrochemical processes are assessed through spatial mapping of temperature and chemical composition (Ca:Na ratios, F, As, B, Li, δ 2 H, and δ 18 O), Piper and bivariate plots, Spearman rank-order correlations, and flow-line modeling of mineral saturation (PHREEQC software). Results help delineate recharge and discharge areas and demonstrate spatial correspondence of major rift structures with changes in chemical and thermal data. Thermal wells with anomalous discharge temperatures and regional thermal gradients exceeding 40 °C/km align with structural boundaries of the Cañada Ancha graben and Caja del Rio horst. Mg-Li geothermometry characterizes temperatures associated with deep circulating groundwater. Important features of the conceptual model are (1) a forced convection system in the Tesuque aquifer associated with the Caja del Rio horst drives upward flow and discharge of warm, Na-rich groundwater in the western half of the Cañada Ancha graben; and (2) major horst-graben structures concentrate upward flow of deep, NaSO 4 thermal waters from underlying bedrock. Both features likely contribute to chemical anomalies and thermal disturbances in the shallow Tesuque aquifer.

We used tephrochronology for upper Neogene deposits in the Española Basin and the adjoining Jemez Mountains volcanic field in the Rio Grande rift, northern New Mexico, to correlate key tephra strata in the study area, identify the sources for many of these tephra, and refine the maximum age of an important stratigraphic unit. Electron-microprobe analyses on volcanic glass separated from 146 pumice-fall, ash-fall, and ash-flow tephra units and layers show that they are mainly rhyolites and dacites. Jemez Mountains tephra units range in age from Miocene to Quaternary. From oldest to youngest these are: (1) the Canovas Canyon Rhyolite and the Paliza Canyon Formation of the lower Keres Group (ca. <12.4–7.4 Ma); (2) the Peralta Tuff Member of the Bearhead Rhyolite of the upper Keres Group (ca. 6.96–6.76 Ma); (3) Puye Formation tephra layers (ca. 5.3–1.75 Ma); (4) the informal San Diego Canyon ignimbrites (ca. 1.87–1.84 Ma); (5) the Otowi Member of the Bandelier Tuff, including the basal Guaje Pumice Bed (both ca. 1.68–1.61 Ma); (6) the Cerro Toledo Rhyolite (ca. 1.59–1.22 Ma); (7) the Tshirege Member of the Bandelier Tuff, including the basal Tsankawi Pumice Bed (both ca. 1.25–1.21 Ma); and (8) the El Cajete Member of the Valles Rhyolite (ca. 60–50 ka). The Paliza Canyon volcaniclastic rocks are chemically variable; they range in composition from dacite to dacitic andesite and differ in chemical composition from the younger units. The Bearhead Rhyolite is highly evolved and can be readily distinguished from the younger units. Tuffs in the Puye Formation are dacitic rather than rhyolitic in composition, and their glasses contain significantly higher Fe, Ca, Mg, and Ti, and lower contents of Si, Na, and K. We conclude that the Puye is entirely younger than the Bearhead Rhyolite and that its minimum age is ca. 1.75 Ma. The San Diego Canyon ignimbrites can be distinguished from all members of the overlying Bandelier Tuff on the basis of Fe and Ca. The Cerro Toledo tephra layers are readily distinguishable from the overlying and underlying units of the Bandelier Tuff primarily by lower Fe and Ca contents. The Tshirege and Otowi Members of the Bandelier Tuff are difficult to distinguish from each other on the basis of electron-microprobe analysis of the volcanic glass; the Tshirege Member contains on average more Fe than the Otowi Member. Tephra layers in the Española Basin that correlate to the Lava Creek B ash bed (ca. 640 ka) and the Nomlaki Tuff (Member of the Tuscan and Tehama Formations, ca. 3.3 Ma) indicate how far tephra from these eruptions traveled (the Yellowstone caldera of northwestern Wyoming and the southern Cascade Range of northern California, respectively). Tephra layers of Miocene age (16–10 Ma) sampled from the Tesuque Formation of the Santa Fe Group in the Española Basin correlate to sources associated with the southern Nevada volcanic field (Timber Mountain, Black Mountain, and Oasis Valley calderas) and the Snake River Plain–Yellowstone hot spot track in Idaho and northwestern Wyoming. Correlations of these tephra layers across the Santa Clara fault provide timelines through various stratigraphic sections despite differences in stratigraphy and lithology. We use tephra correlations to constrain the age of the base of the Ojo Caliente Sandstone Member of the Tesuque Formation to 13.5–13.3 Ma.

Two- and three-dimensional electrical resistivity models derived from the magnetotelluric method were interpreted to provide more accurate hydrogeologic parameters for the Albuquerque and Española Basins. Analysis and interpretation of the resistivity models are aided by regional borehole resistivity data. Examination of the magnetotelluric response of hypothetical stratigraphic cases using resistivity characterizations from the borehole data elucidates two scenarios where the magnetotelluric method provides the strongest constraints. In the first scenario, the magnetotelluric method constrains the thickness of extensive volcanic cover, the underlying thickness of coarser-grained facies of buried Santa Fe Group sediments, and the depth to Precambrian basement or overlying Pennsylvanian limestones. In the second scenario, in the absence of volcanic cover, the magnetotelluric method constrains the thickness of coarser-grained facies of buried Santa Fe Group sediments and the depth to Precambrian basement or overlying Pennsylvanian limestones. Magnetotelluric surveys provide additional constraints on the relative positions of basement rocks and the thicknesses of Paleozoic, Mesozoic, and Tertiary sedimentary rocks in the region of the Albuquerque and Española Basins. The northern extent of a basement high beneath the Cerros del Rio volcanic field is delineated. Our results also reveal that the largest offset of the Hubbell Spring fault zone is located 5 km west of the exposed scarp. By correlating our resistivity models with surface geology and the deeper stratigraphic horizons using deep well log data, we are able to identify which of the resistivity variations in the upper 2 km belong to the upper Santa Fe Group sediments

The structural geometry of transfer and accommodation zones that relay strain between extensional domains in rifted crust has been addressed in many studies over the past 30 years. However, details of the kinematics of deformation and related stress changes within these zones have received relatively little attention. In this study we conduct the first-ever systematic, multi-basin fault-slip measurement campaign within the late Cenozoic Rio Grande rift of northern New Mexico to address the mechanisms and causes of extensional strain transfer associated with a broad accommodation zone. Numerous (562) kinematic measurements were collected at fault exposures within and adjacent to the NE-trending Santo Domingo Basin accommodation zone, or relay, which structurally links the N-trending, right-stepping en echelon Albuquerque and Española rift basins. The following observations are made based on these fault measurements and paleostresses computed from them. (1) Compared to the typical northerly striking normal to normal-oblique faults in the rift basins to the north and south, normal-oblique faults are broadly distributed within two merging, NE-trending zones on the northwest and southeast sides of the Santo Domingo Basin. (2) Faults in these zones have greater dispersion of rake values and fault strikes, greater dextral strike-slip components over a wide northerly strike range, and small to moderate clockwise deflections of their tips. (3) Relative-age relations among fault surfaces and slickenlines used to compute reduced stress tensors suggest that far-field, ~E-W–trending σ 3 stress trajectories were perturbed 45° to 90° clockwise into NW to N trends within the Santo Domingo zones. (4) Fault-stratigraphic age relations constrain the stress perturbations to the later stages of rifting, possibly as late as 2.7–1.1 Ma. Our fault observations and previous paleomagnetic evidence of post–2.7 Ma counterclockwise vertical-axis rotations are consistent with increased bulk sinistral-normal oblique shear along the Santo Domingo rift segment in Pliocene and later time. Regional geologic evidence suggests that the width of active rift faulting became increasingly confined to the Santo Domingo Basin and axial parts of the adjoining basins beginning in the late Miocene. We infer that the Santo Domingo clockwise stress perturbations developed coevally with the oblique rift segment mainly due to mechanical interactions of large faults propagating toward each other from the adjoining basins as the rift narrowed. Our results suggest that negligible bulk strike-slip displacement has been accommodated along the north-trending rift during much of its development, but uncertainties in the maximum ages of fault slip do not allow us to fully evaluate and discriminate between earlier models that invoked northward or southward rotation and translation of the Colorado Plateau during early (Miocene) rifting.

We investigated a Plio-Pleistocene alluvial succession in the Albuquerque Basin of the Rio Grande rift in New Mexico using geomorphic, stratigraphic, sedimentologic, geochronologic, and magnetostratigraphic data. New 40 Ar/ 39 Ar age determinations and magnetic-polarity stratigraphy refine the ages of the synrift Santa Fe Group. The Pliocene Ceja Formation lies on the distal hanging-wall ramp across much of the Albuquerque Basin. The Ceja onlapped and buried a widespread, Upper Miocene erosional paleosurface by 3.0 Ma. Sediment accumulation rates in the Ceja Formation decreased after 3.0 Ma and the Ceja formed broad sheets of amalgamated channel deposits that prograded into the basin after ca. 2.6 Ma. Ceja deposition ceased shortly after 1.8 Ma, forming the Llano de Albuquerque surface. Deposition of the Sierra Ladrones Formation by the ancestral Rio Grande was focused near the eastern master fault system before piedmont deposits (Sierra Ladrones Formation) began prograding away from the border faults between 1.8 and 1.6 Ma. Widespread basin filling ceased when the Rio Grande began cutting its valley, shortly after 0.78 Ma. Although the Albuquerque Basin is tectonically active, the development of through-going drainage of the ancestral Rio Grande, burial of Miocene unconformities, and coarsening of upper Santa Fe Group synrift basin fill were likely driven by climatic changes. Valley incision was approximately coeval with increased northern- hemisphere climatic cyclicity and magnitude and was also likely related to climatic changes. Asynchronous progradation of coarse-grained, margin-sourced detritus may be a consequence of basin shape, where the basinward tilting of the hanging wall promoted extensive sediment bypass of coarse-grained, margin-sourced sediment across the basin.

Discrepancies among previous models of the geometry of the Albuquerque Basin motivated us to develop a new model using a comprehensive approach. Capitalizing on a natural separation between the densities of mainly Neogene basin fill (Santa Fe Group) and those of older rocks, we developed a three-dimensional (3D) geophysical model of syn-rift basin-fill thickness that incorporates well data, seismic-reflection data, geologic cross sections, and other geophysical data in a constrained gravity inversion. Although the resulting model does not show structures directly, it elucidates important aspects of basin geometry. The main features are three, 3–5-km-deep, interconnected structural depressions, which increase in size, complexity, and segmentation from north to south: the Santo Domingo, Calabacillas, and Belen subbasins. The increase in segmentation and complexity may reflect a transition of the Rio Grande rift from well-defined structural depressions in the north to multiple, segmented basins within a broader region of crustal extension to the south. The modeled geometry of the subbasins and their connections differs from a widely accepted structural model based primarily on seismic-reflection interpretations. Key elements of the previous model are an east-tilted half-graben block on the north separated from a west-tilted half-graben block on the south by a southwest-trending, scissor-like transfer zone. Instead, we find multiple subbasins with predominantly easterly tilts for much of the Albuquerque Basin, a restricted region of westward tilting in the southwestern part of the basin, and a northwesterly trending antiform dividing subbasins in the center of the basin instead of a major scissor-like transfer zone. The overall eastward tilt indicated by the 3D geophysical model generally conforms to stratal tilts observed for the syn-rift succession, implying a prolonged eastward tilting of the basin during Miocene time. An extensive north-south synform in the central part of the Belen subbasin suggests a possible path for the ancestral Rio Grande during late Miocene or early Pliocene time. Variations in rift-fill thickness correspond to pre-rift structures in several places, suggesting that a better understanding of pre-rift history may shed light on debates about structural inheritance within the rift.

A P-wave seismic velocity model derived from analysis of data from a seismic refraction/reflection survey provides the first regional-scale view of the subsurface structure of the upper crust of the southern Rio Grande rift. The seismic survey was conducted along a 205-km-long east-west transect that crosses a portion of the Basin and Range province and the Rio Grande rift in southernmost New Mexico and west Texas. In the upper few kilometers of the crust, the model shows a series of low- and high-velocity zones that correlate with the mid-Tertiary to Holocene Basin and Range structure at the surface. Typically, basins are 10–20 km wide and ~1–3 km deep. Beneath the mid-Tertiary to Holocene structures are velocity anomalies interpreted to be the result of tectonic activity of Paleozoic to early Tertiary age. The geometry of high-velocity zones at 3–10 km depth in the eastern half of the model correlates well with Laramide block uplifts mapped at the surface and suggests that the thrusts that bound these blocks may sole into a regional detachment at 10–15 km depth. In the western half of the model, a low-velocity zone that reaches depths as great as 11 km has geometry highly suggestive of a deep basin. This previously unrecognized feature may represent a combination of great thicknesses of the Upper Jurassic to Lower Cretaceous rocks of the Chihuahua Trough and Bisbee Basin stacked above rocks of the older Paleozoic Pedregosa Basin. Alternatively, portions of this region may have been thickened by thrust stacking during the Laramide orogeny.