Abstract: The discontinuities developed in grotto rock mass are multiple in scales and complex in distribution, which pose a serious threat to the long-term stability of grotto temples. A set of integrated detection methods formed by UAV, 3D laser scanner, ground penetrating radar, camera and manual survey statistics were applied in Anyue Yuanjue Cave, and the parameters and spatial distribution characteristics of multi-scale discontinuities of the grotto temple were systematically analyzed. Results show that large scale unloading fractures and bedding planes as well as medium and small scale weathering fractures are mainly developed in grotto temple rock mass. And the discontinuities on both the surface and interior of grotto roof present a spatially differentiated distribution characteristic of more on the north side than on the south. In addition, the zoning statistics reveal that the occurrence of discontinuities of grotto roof conforms to Fisher distribution, and the discontinuities at different scales develop a trend of nearly NW330°, with a consistent developmental rule. While the trace length, spacing, and depth of discontinuities mainly obey the negative exponential or log-normal distribution. On this basis, integrating Agisoft Metashape-Cloudcompare-Geomagic Studio-HyperMesh-3DEC multi-software processing, a 3D modeling technology for multi-source discontinuities information of grotto temple is proposed. Then, the digital distribution information of discontinuities of grotto temple rock mass is obtained. This paper can provide a basis for analyzing key blocks and overall stability of grotto temples, as well as for developing of reinforcement schemes.