Tracing short-lived hydrothermal circulation systems and water-rock interactions around small-scale intrusions

Hydrothermal circulation around magmatic intrusions affects hazards, critical metal deposition, and heat transfer, however it is often difficult to directly observe. We combine stable isotope and thermochronologic measurements, hydrothermal circulation modeling, and forward thermal history modeling to understand the short-lived hydrothermal circulation system that developed around the Little Devils Postpile, a small-scale mafic intrusion in the eastern Sierra Nevada Mountains of western North America. Bedrock samples collected in transects orthogonal to the intrusion reveal reset and partially reset apatite (U-Th-Sm)/He and zircon (U-Th)/He ages up to similar to 13m from the contact. Measurements of hydrogen isotopes in apatite crystals from the same samples show deuterium depletion near the contact; delta D values vary between -120 parts per thousand and -132 parts per thousand within 3m of the contact, -112 parts per thousand to -104 parts per thousand at 4-6m from the contact, and -85 parts per thousand to -95 parts per thousand at distances of 12m to 40m from the contact. In contrast, oxygen isotopes do not vary spatially with a delta O-18 of +6.4 +/- 0.4 parts per thousand. This depletion of delta D and consistency of delta O-18 suggests that short-lived water-rock interaction occurred following the intrusion emplacement, facilitating the exchange of hydrogen isotopes via diffusion rather than dissolution-precipitation processes. Conductive cooling calculations predict an alteration aureole ten times larger than is observed by the low-temperature thermochronometric ages. Numerical models of hydrothermal circulation, for a range of host bedrock permeabilities and magmatic emplacement temperatures, reveal that the geometry of alteration is most consistent with buoyancy-driven hydrothermal circulation and a relatively high crustal permeability. With this case study, we show that combining thermochronology and hydrothermal modeling allows for the opportunity to model hydrogen isotope exchange and predict plausible apatite-water hydrogen fractionation. These new constraints on the geometry of water-rock interactions coupled with hydrothermal circulation modeling provide (1) a rare window into the pressure-temperature evolution of a small magmatic intrusion, (2) constraints on crustal permeability, and (3) a useful tool for interpreting low-temperature thermochronometric data in volcanic settings.