An original method is presented that allows us to measure the local shear-wave birefringence properties over any depth interval. It requires the acquisition of two shear-wave vertical seismic profiles (VSPs), each with different initial polarizations of the shear wave. The method is based on the estimation of a two by two matrix (called the propagator matrix) that represents a linear operator between two states of polarization. No information is required about layering above the zone of interest (in particular, about the weathering zone). If these two states of polarization correspond to the direct downgoing shear wave at two different depths z 1 and z2, the operator represents the transmission properties between the two depths. Under the previous hypothesis, this operator is independent of the source polarization and can be accurately estimated by a least-squares method in the frequency domain. Physically, this operator is a multicomponent deconvolution, whose column vectors represent the state of polarizations at a depth Z2 for two linear and INTRODUCTION