No abstract

No abstract

Abstract. New results are presented from the teleseismic component of the Jemez Tomography Experiment conducted across Valles caldera in northern New Mexico. We invert 4872 relative P wave arrival times recorded on 50 portable stations to determine velocity structure to depths of 40 km. The three principle features of our model for Valles caldera are: (1) near-surface low velocities of-17% beneath the Toledo embayment and the Valle Grande, (2) midcrustal low velocities of-23% in an ellipsoidal volume underneath the northwest quadrant of the caldera, and (3) a broad zone of low velocities (-15%) in the lower crust or upper mantle. Crust shallower than 20 km is generally fast to the northwest of the caldera and slow to the southeast. Near-surface low velocities are interpreted as thick deposits of Bandelier tuff and postcaldera volcaniclastic rocks. Lateral variation in the thickness of these deposits supports increased caldera collapse to the southeast, beneath the Valle Grande. We interpret the midcrustal low-velocity zone to contain a minimum melt fraction of 10%. While we cannot rule out the possibility that this zone is the remnant 1.2 Ma Bandelier magma chamber, the eruption history and geochemistry of the volcanic rocks erupted in Valles caldera following the Bandelier tuff make it more likely that magma results from a new pulse of intrusion, indicating that melt flux into the upper crust beneath Valles caldera continues. The low-velocity zone near the crust-mantle boundary is consistent with either partial melt in the lower crust or mafic rocks without partial melt in the upper mantle. In either case, this low-velocity anomaly indicates that underplating by mantle-derived melts has occurred.

A new approximate algorithm for two-point ray tracing is proposed and tested in a variety of laterally heterogeneous velocity models. An initial path estimate is perturbed using a geometric interpretation of the ray equations, and the travel time along the path is minimized in a piecewise fashion. This perturbation is iteratively performed until the travel time converges within a specified limit. Test results show that this algorithm successfully finds the correct travel time within typical observational error much faster than existing three-dimensional ray tracing programs. The method finds an accurate ray path in a fully three-dimensional form even where lateral variations in velocity are severe. Because our algorithm utilizes direct minimization of the travel time instead of solving the ray equations, a simple linear interpolation scheme can be employed to compute velocity as a function of position, providing an added computational advantage.