We prove a new upper bound on the second moment of Maass form symmetric square L-functions defined over Gaussian integers. Combining this estimate with the recent result of Balog–Biro–Cherubini–Laaksonen, we improve the error term in the prime geodesic theorem for the Picard manifold.