Critical zone (CZ) structure, including the spatial distribution of minerals, elements, and fluid-filled pores, evolves on geologic time scales resulting from both top-down climatic forcing and bottom-up geologic controls. Climate and lithology may be imprinted in subsurface structure as depth-dependent trends in geophysical, geochemical, mineralogical, and biological datasets. As the weathering profile is as much (or more) a product of past environmental conditions, development of predictive models requires understanding the relative roles of climatic forcing and the geologic template on which CZ processes evolve. Doing so in complex volcanic terrains with high initial bedrock porosity and distinct depositional and hydrothermal alteration histories is particularly challenging. To resolve CZ structure in a rhyolitic catchment in the Valles Caldera National Preserve (NM, USA), this study combined geophysics, drilling, and laboratory analyses to produce depth-resolved porosity, geochemistry, and mineralogy datasets to >40 m in depth. Quantitative X-ray diffraction analysis showed that local mineral transformations control complex chemical enrichment/depletion (τ) patterns. Using linear discriminant analysis, key variables enabled separation of complex-layered geology into discrete zones. Contemporary, matrix-dominated weathering processes and modern hydrologic fluxes occur dominantly within the top 15 m of the weathering profile. This zone is convoluted by incomplete primary mineral weathering and overprinted by post-eruption weathering and metasomatism. Matrix weathering transitions to fracture surface weathering driven by deep percolation of slower moving, longer residence time meteoric waters at depth. By altering initial conditions and weathering trajectory, geologic legacy is a critical factor in how this subsurface landscape evolved and functions.

Coda Q -1 was calculated at Nevado del Ruiz Volcano (NRV) before and after two phreatomagmatic eruptions (November 1985, September 1989 and for a period of stability (May 1987) using a functional form for coda derived from a single scattering model (Sato 1977). Substantial changes were found before and after the eruptions. The highest value of Q -1 was found during the November 1985 eruption, an intermediate value for the September 1989 eruption, and the lowest value for May 1987. It seems that the changes in coda Q -1 at NRV have a still-unknown relationship with the volume or magnitude of the eruption. A relatively strong frequency dependence was found for all periods. Also Q -1 clearly changed with time, suggesting that the scattering was strong for the eruption of November 1985 and decreased for the eruption of September 1989, and that the intrinsic absorption probably increased. This suggests the possibility that crystallization is taking place in the NRV magma. The clear change of coda Q -1 before and after the eruptions at NRV also suggests the possibility that coda Q -1 is a premonitory tool of activity at this volcano.