Method and system to identify, verify and remediate subterranean anomalies. Infrared (IR) thermographic scanning [10] of a selected surface area obtains image area data [12, 14] inferring such an anomaly. Ground penetrating radar (GPR) is used [22] at predetermined surface locations penetrates subterraneously to a depth including the anomaly, creating vertical dimension radar data [120] showing anomaly depth. Surface image IR data is correlated with GPR data to verify the anomaly and its vertical dimension and finds [124] a central location in the anomaly. Anomaly volume is predetermined from the area data and vertical-dimension data. Grout injected [26, 124] into the central location at first pressure secures the anomaly by surrounding the central location. After verifying centrally securement, more grout is injected [32, 126] into the anomaly region at second pressure at least as great as the first pressure until the total grout injected approximates the predetermined anomaly volume.