Partial melt in the upper-middle crust of the northwest Himalaya revealed by Rayleigh wave dispersion

Abstract

Seismic shear-wave velocities are sensitive to the partial melts that should be present in the Himalayan orogen if low-viscosity channel flow is active at the present day. We analyzed regional earthquakes in the western Himalaya and Tibet recorded on 16 broadband seismometers deployed across the NW Indian Himalaya, from the Indian platform to the Karakoram Range. We used a multiple filter technique to calculate the group velocity dispersion of fundamental-mode Rayleigh waves, and then inverted the dispersion records to obtain separate one-dimensional shear-wave velocity models for five geologic provinces: the Tibetan plateau, Ladakh arc complex, Indus Tsangpo suture zone, Tethyan Himalaya, and Himalayan thrust belt. Our velocity models show a low-velocity layer (LVL) with 7–17% velocity reduction centered at ~ 30 km depth and apparently continuous from the Tethyan Himalaya to the Tibetan plateau. This LVL shows good spatial correspondence with observations of low resistivity from magnetotelluric studies along the same profile. Of the possible explanations for low velocity and low resistivity in the mid-crust, only the presence of melts or aqueous fluids (or both) satisfactorily explains both sets of observations. Elevated heat flow observed in the NW Himalaya implies that if aqueous fluids are present in the mid-crust, then the mid-crust is well above its solidus. Comparison of our results with laboratory measurements and theoretical models suggests 3–7% melt is present in a channel in the upper-middle crust of the NW Himalaya at the present day, and the physical conditions to enable active channel flow may be present.