The evolution of metamorphic core complexes is closely related to gold formation. Due to the thick cover and orebodies’ lying depth, exploration regularities, and metallogenic prediction have faced challenges. Therefore, new data were obtained through controlled-source audio magnetotelluric (MT) and broadband MT methods to predict the ore-bearing properties at depth and decipher the spatial coupling relationship, respectively. The results of broadband MT surveys allow us to identify the thick cover (apparent resistivity of 3,000 Ω m), the strongly deformed metamorphic rock (apparent resistivity of 18,000 Ω m), and two low-angle detachment faults (apparent resistivity from 5,000 to 8,000 Ω m). The fault-dip value on the left side is roughly 30°, while on the other side, the values significantly change from 30° to 75°, although they are in the same structure. Moreover, the controllable source audio MT method provides more detailed examinations of the right side fault. The proved ore-bearing gold orebodies were all located in the transition zone where the inclination angle changes from steep to gentle or gentle to steep. The discovery of this mineralization pattern allows us to evaluate the mining prospecting potential and predict the deep-seated metallogenic location. Based on the geotectonic environment and the geophysical profiles, we conclude that the coeval exhumation of the metamorphic core complex with non-symmetrical northwest-southeast shear senses may have resulted from a decratonization event during the retreat of the Paleo-Pacific Plate.

Three-dimensional (3D) geological modeling is an essential tool for assessing digital mine safety and detecting underground cavities. This study systematically reviews the capabilities, advantages, and limitations of mainstream 3D geological modeling software, integrating orthophoto imagery and ground-penetrating radar (GPR) data with the GOCAD platform to construct a three-dimensional model of shallow cavity-related anomalies in Qinglong Manchu Autonomous County, Hebei Province. The GPR survey included 20 reconnaissance profiles and eight detailed profiles in the priority exploration area. Two principal target zones were identified from the measured GPR data: Target Zone 1 near C0071 on Line C at an interpreted depth of approximately 5–25 m, and Target Zone 2 near D0090 on Line D at approximately 3–12 m depth. A laterally continuous reflective interface was also identified at approximately 10 m depth and interpreted as the boundary between Quaternary unconsolidated sediments and the underlying bedrock. Orthophoto imagery provides surface constraints, while GPR data provides subsurface structural information. Their integration in GOCAD provides a common three-dimensional spatial framework for relating surface features, measured GPR anomalies, and field-verification results, rather than a quantitatively demonstrated improvement in detection accuracy. Field investigation and excavation at Target Zone 2 revealed a dolomite-hosted cave-fissure space filled with calcareous clay, angular dolostone breccia, and mammalian fossils, supporting the geological interpretation of the GPR anomaly. The detailed inter-profile geometry remains dependent on survey-line distribution and interpolation.