To calibrate airborne gravity gradiometers currently in development in China, it is urgent to build an airborne gravity gradiometer test site. The site’s selection depends on the preknowledge of high-resolution gravity and gradient structures. The residual terrain modelling (RTM) technique is generally applied to recover the short-scale gravity field signals. However, due to limitations in the quality and resolution of density models, RTM terrain generally assumes a constant density. This assumption can introduce significant errors in areas with substantial density anomalies and of reggued terrain, such as volcano areas. In this study, we promote a method to determine a high-resolution gravity field by integrating long-wavelength signals generated by EGM2008 with short-wavelength signals from terrain relief and shallow density anomalies. These short wavelength signals are recovered using the RTM technique with both constant density and density anomalies obtained through the equivalent source layer (ESL) method, utilizing sparse terrestrial gravity measurements. Compared to the recovery rate of 54.62% using the classical RTM method, the recovery rate increases to 86.22% after involving density anomalies. With this method, we investigate the gravity field signals over the Wudalianchi Volcano Field (WVF) both on the Earth’s surface and at a flight height of 100 m above the terrain. The contribution of each part and their attenuation characters are studied. In particular, the 5 km × 5 km area surrounding Bijiashan (BJS) and Wohushan (WHS) volcanos shows a strong gravity signature, making it a good candidate for the test site location. This study gives the location of the airborne gravity gradiometer test site which is an essential step in the instruments’ development. Furthermore, the method presented in this study offers a foundational framework for future data processing within the test site.