Index transforms with the squares of Kelvin functions

New index transforms, involving squares of Kelvin functions, are investigated. Mapping properties and inversion formulas are established for these transforms in Lebesgue spaces. The results are applied to solve a boundary value problem on the wedge for a fourth order partial differential equation.