  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.
