Jonathan Saul Caine

Brittle structures and their role in controlling porosity and permeability in a complex Precambrian crystalline-rock aquifer system in the Colorado Rocky Mountain Front Range

Expansion of the Denver metropolitan area has resulted in substantial residential development in the foothills of the Rocky Mountain Front Range. This type of suburban growth, characteristic of much of the semiarid intermountain west, often relies on groundwater from individual domestic wells and is exemplified in the Turkey Creek watershed. The watershed is underlain by complexly deformed and fractured crystalline bedrock in which groundwater resources are poorly understood, and concerns regarding groundwater mining and degradation have arisen. As part of a pilot project to establish quantitative bounds on the groundwater resource, an outcrop-based geologic characterization and numerical modeling study of the brittle structures and their controls on the flow system was initiated. Existing data suggest that groundwater storage, flow, and contaminant transport are primarily controlled by a heterogeneous array of fracture networks. Inspections of well-permit data and field observations led to a conceptual model in which three dominant lithologic groups underlying sparse surface deposits form the aquifer system—metamorphic rocks, a complex array of granitic intrusive rocks, and major brittle fault zones. Pervasive but variable jointing of each lithologic group forms the “background” permeability structure and is an important component of the bulk storage capacity. This “background” is cut by brittle fault zones of varying structural styles and by pegmatite dikes, both with much higher fracture intensities relative to “background” that likely make them spatially complex conduits. Probabilistic, discrete-fracture-network and finite-element modeling was used to estimate porosity and permeability at the outcrop scale using fracture network data collected in the field. The models were conditioned to limited aquifer test and borehole geophysical data and give insight into the relative hydraulic properties between locations and geologic controls on storage and flow. Results from this study reveal a complex aquifer system in which the upper limits on estimated hydraulic properties suggest limited storage capacity and permeability as compared with many sedimentary-rock and surficial-deposit aquifers.

Structural and geochemical characteristics of faulted sediments and inferences on the role of water in deformation, Rio Grande Rift, New Mexico

The San Ysidro fault is a spectacularly exposed normal fault located in the northwestern Albuquerque Basin of the Rio Grande Rift. This intrabasin fault is representative of many faults that formed in poorly lithified sediments throughout the rift. The fault is exposed over nearly 10 km and accommodates nearly 700 m of dip slip in subhorizontal, siliciclastic sediments. The extent of the exposure facilitates study of along-strike variations in deformation mechanisms, architecture, geochemistry, and permeability. The fault is composed of structural and hydrogeologic components that include a clay-rich fault core, a calcite-cemented mixed zone, and a poorly developed damage zone primarily consisting of deformation bands. Structural textures suggest that initial deformation in the fault occurred at low temperature and pressure, was within the paleosaturated zone of the evolving Rio Grande Rift, and was dominated by particulate flow. Little geochemical change is apparent across the fault zone other than due to secondary processes. The lack of fault-related geochemical change is interpreted to reflect the fundamental nature of water-saturated, particulate flow. Early mechanical entrainment of low-permeability clays into the fault core likely caused damming of groundwater flow on the up-gradient, footwall side of the fault. This may have caused a pressure gradient and flow of calcite-saturated waters in higher-permeability, fault-entrained siliciclastic sediments, ultimately promoting their cementation by sparry calcite. Once developed, the cemented and clay-rich fault has likely been, and continues to be, a partial barrier to cross-fault groundwater flow, as suggested by petrophysical measurements. Aeromagnetic data indicate that there may be many more unmapped faults with similar lengths to the San Ysidro fault buried within Rio Grande basins. If these buried faults formed by the same processes that formed the San Ysidro fault and have persistent low-permeability cores and cemented mixed zones, they could compartmentalize the basin-fill aquifers more than is currently realized, particularly if pumping stresses continue to increase in response to population growth.

Controls on groundwater flow in a semiarid folded and faulted intermountain basin

The major processes controlling groundwater flow in intermountain basins are poorly understood, particularly in basins underlain by folded and faulted bedrock and under regionally realistic hydrogeologic heterogeneity. To explore the role of hydrogeologic heterogeneity and poorly constrained mountain hydrologic conditions on regional groundwater flow in contracted intermountain basins, a series of 3-D numerical groundwater flow models were developed using the South Park basin, Colorado, USA as a proxy. The models were used to identify the relative importance of different recharge processes to major aquifers, to estimate typical groundwater circulation depths, and to explore hydrogeologic communication between mountain and valley hydrogeologic landscapes. Modeling results show that mountain landscapes develop topographically controlled and predominantly local-scale to intermediate-scale flow systems. Permeability heterogeneity of the fold and fault belt and decreased topographic roughness led to permeability controlled flow systems in the valley. The structural position of major aquifers in the valley fold and fault belt was found to control the relative importance of different recharge mechanisms. Alternative mountain recharge model scenarios showed that higher mountain recharge rates led to higher mountain water table elevations and increasingly prominent local flow systems, primarily resulting in increased seepage within the mountain landscape and nonlinear increases in mountain block recharge to the valley. Valley aquifers were found to be relatively insensitive to changing mountain water tables, particularly in structurally isolated aquifers inside the fold and fault belt.

A GIS and statistical approach to identify variables that control water quality in hydrothermally altered and mineralized watersheds, Silverton, Colorado, USA

Hydrothermally altered bedrock in the Silverton mining area, southwest Colorado, USA, contains sulfide minerals that weather to produce acidic and metal-rich leachate that is toxic to aquatic life. This study utilized a geographic information system (GIS) and statistical approach to identify watershed-scale geologic variables in the Silverton area that influence water quality. GIS analysis of mineral maps produced using remote sensing datasets including Landsat Thematic Mapper, advanced spaceborne thermal emission and reflection radiometer, and a hybrid airborne visible infrared imaging spectrometer and field-based product enabled areas of alteration to be quantified. Correlations between water quality signatures determined at watershed outlets, and alteration types intersecting both total watershed areas and GIS-buffered areas along streams were tested using linear regression analysis. Despite remote sensing datasets having varying watershed area coverage due to vegetation cover and differing mineral mapping capabilities, each dataset was useful for delineating acid-generating bedrock. Areas of quartz–sericite–pyrite mapped by AVIRIS have the highest correlations with acidic surface water and elevated iron and aluminum concentrations. Alkalinity was only correlated with area of acid neutralizing, propylitically altered bedrock containing calcite and chlorite mapped by AVIRIS. Total watershed area of acid-generating bedrock is more significantly correlated with acidic and metal-rich surface water when compared with acid-generating bedrock intersected by GIS-buffered areas along streams. This methodology could be useful in assessing the possible effects that alteration type area has in either generating or neutralizing acidity in unmined watersheds and in areas where new mining is planned.

A comprehensive survey of faults, breccias, and fractures in and flanking the eastern Española Basin, Rio Grande rift, New Mexico

A comprehensive survey of geologic structures formed in the Earth’s brittle regime in the eastern Española Basin and flank of the Rio Grande rift, New Mexico, reveals a complex and protracted record of multiple tectonic events. Data and analyses from this representative rift flank-basin pair include measurements from 53 individual fault zones and 22 other brittle structures, such as breccia zones, joints, and veins, investigated at a total of just over 100 sites. Structures were examined and compared in poorly lithified Tertiary sedi ments, as well as in Paleozoic sedimentary and Proterozoic crystalline rocks. Data and analyses include geologic maps; field observations and measurements; orientation, kinematic, and paleostress analyses; statistical examination of fault trace lengths derived from aeromagnetic data; mineralogy and chemistry of host and fault rocks; and investigation of fault versus bolideimpact hypotheses for the origin of enigmatic breccias found in the Proterozoic basement rocks. Fault kinematic and paleostress analyses suggest a record of transitional, and perhaps partitioned, strains from the Laramide orogeny through Rio Grande rifting. Normal faults within Tertiary basin-fill sediments are consistent with more typical WNW-ESE Rio Grande rift extension, perhaps de coupled from bedrock structures due to strength contrasts favoring the formation of new faults in the relatively weak sediments. Analyses of the fault-length data indicate power-law length distributions similar to those reported from many geologic settings globally. Mineralogy and chemistry in Proterozoic fault-related rocks reveal geochemical changes tied to hydro thermal alteration and nearly isochemical transformation of feldspars to clay minerals. In sediments, faulted minerals are characterized by mechanical entrainment with minor secondary chemical changes. Enigmatic breccias in rift-flanking Protero zoic rocks are autoclastic and isochemical with respect to their protoliths and exist near shatter cones believed to be related to a previously reported pre-Pennsylvanian impact event. A weak iridium anomaly is associated with the breccias as well as adjacent protoliths, thus an impact shock wave cannot be ruled out for their origin. Major fault zones along the eastern rift-flank mountain front are discontinuous and unlikely to impede regional groundwater flow into Española Basin aquifers. The breccia bodies are not large enough to constitute aquifers, and no faultor breccia-related geochemical anomalies were identified as potential contamination sources for ground or surface waters. The results of this work provide a broad picture of structural diversity and tectonic evolution along the eastern flank of the central Rio Grande rift and the adjacent Española Basin representative of the rift as a whole and many rifts worldwide.