We study the regularity of solutions to Dirichlet, Neumann and mixed boundary value problems in polyhedral cones K ⊂ R 3 in the specific scale B α τ,τ , 1 τ = α d + 1 p of Besov spaces. The regularity of the solution in these spaces determines the order of approximation that can be achieved by adaptive and nonlinear numerical schemes. We show that the solutions are much smoother in this specific Besov scale compared to the ususal L 2 -Sobolev scale, which justifies the use of adaptive schemes. Our findings considerably extend [12, Thms. 2.1, 3.1]. The proofs are performed by combining Sobolev estimates with mixed weights -which measure the distance to the vertex and the edges of the cone, respectively -with characterizations of Besov spaces by wavelet expansions.