A computer-implemented method to calculate an orientation of a multicomponent seismic sensor (216, 290, 450), which is deployed on a seabed at a fluid-solid interface between a first and a second elastic medium, wherein a subsea formation below the seabed (208, 288, 303) forms the first elastic medium and seawater forms the second elastic medium; the method comprising: receiving information associated with the interface between the first elastic medium and the second elastic medium, wherein the information comprises sensor data (510, 520, 530, 550) acquired by the multicomponent seismic sensor (216, 290, 450); and based on at least a portion of the information, estimating seismic wave properties that comprise elastic properties, depth-dependent properties and horizontal slowness; wherein the estimation of the depth-dependent properties comprises: calculating vertical slowness in the first elastic medium; estimating a polarization angle of an incident seismic wave in the first elastic medium corresponding to the calculated vertical slowness; calculating an angle of incidence in the second elastic medium; and inducing a computing system (350, 354) to calculate the orientation of the multicomponent seismic sensor (216, 290, 450) according to the equation: tan θ g − θ h − θ n = p q α 2 − q α 1 p 2 + q α 1 q α 2 where θ g denotes the angle of incidence of a seismic wave, θ h denotes the polarization angle of the incident seismic wave, θ n denotes the orientation of the multicomponent seismic sensor when rotated around an axis perpendicular to the propagation plane, q α 1 is the vertical slowness just above the seabed, q α 2 denotes the vertical slowness just below the seabed and p denotes the horizontal slowness. receiving information associated with the interface between the first elastic medium and the second elastic medium, wherein the information comprises sensor data (510, 520, 530, 550) acquired by the multicomponent seismic sensor (216, 290, 450); and based on at least a portion of the information, estimating seismic wave properties that comprise elastic properties, depth-dependent properties and horizontal slowness; wherein the estimation of the depth-dependent properties comprises: calculating vertical slowness in the first elastic medium; estimating a polarization angle of an incident seismic wave in the first elastic medium corresponding to the calculated vertical slowness; calculating an angle of incidence in the second elastic medium; and inducing a computing system (350, 354) to calculate the orientation of the multicomponent seismic sensor (216, 290, 450) according to the equation: tan θ g − θ h − θ n = p q α 2 − q α 1 p 2 + q α 1 q α 2 where θ g denotes the angle of incidence of a seismic wave, θ h denotes the polarization angle of the incident seismic wave, θ n denotes the orientation of the multicomponent seismic sensor when rotated around an axis perpendicular to the propagation plane, q α 1 is the vertical slowness just above the seabed, q α 2 denotes the vertical slowness just below the seabed and p denotes the horizontal slowness. calculating vertical slowness in the first elastic medium; estimating a polarization angle of an incident seismic wave in the first elastic medium corresponding to the calculated vertical slowness; calculating an angle of incidence in the second elastic medium; and inducing a computing system (350, 354) to calculate the orientation of the multicomponent seismic sensor (216, 290, 450) according to the equation: tan θ g − θ h − θ n = p q α 2 − q α 1 p 2 + q α 1 q α 2 where θ g denotes the angle of incidence of a seismic wave, θ h denotes the polarization angle of the incident seismic wave, θ n denotes the orientation of the multicomponent seismic sensor when rotated around an axis perpendicular to the propagation plane, q α 1 is the vertical slowness just above the seabed, q α 2 denotes the vertical slowness just below the seabed and p denotes the horizontal slowness.