Investigation of Iron Ore Mineral Distribution Using Aero-Magnetic Exploration Techniques: Case Study at Pocheon, Korea

: We present our aeromagnetic survey results from an investigation of the iron ore mineral distribution in Pocheon, Korea, in the west-central area of the Korean Peninsula. A manned aeromagnetic system using a helicopter for regional exploration and an unmanned aeromagnetic system using a multicopter for high-resolution exploration were used for the survey. The inversion results of the magnetic data conﬁrmed the possibility of the existence of a new iron ore body. Drilling was carried out based on inversion results and drilling revealed amphibolite including iron ore, as indicated by a strong magnetic response. The position and depth of the iron ore were consistent with the interpretation results of the magnetic data. installed All data be used to calibrate, and analyze magnetic ﬁndings.


Introduction
Securing adequate mineral resources is important for successfully navigating a new era of industrial growth. For some time, we have conducted studies on titanomagnetite deposits located in the west-central portion of the Korean Peninsula. The area has numerous active mines, with magnetite (an iron ore mineral), titanium (Ti), and vanadium (V) as the main deposits. Magnetite is easy to identify with magnetic surveys, as it has a strong susceptibility. Given that susceptibility is a direct indicator of the magnetic mineral content, magnetic exploration is one of the most efficient means of finding new iron ore sources. As such, the identification of iron ore deposits and their distribution throughout the westcentral part of the Korean Peninsula using magnetic surveys has been and continues to be the focus of our research.
Exploration for new mineral resources is challenging in Korea, given that the target areas are mostly located in complex mountainous terrain. In particular, Mount Gonam, the subject area of this study, has an altitude of 848.8 m and very rugged and steep topography. Exploration of harsh environments such as this requires high-resolution measurement devices and efficient data acquisition. In consideration of the conditions of this region, we used an aeromagnetic exploration method to map the iron ore mineral distribution in Pocheon, Korea.
Aeromagnetic surveys are one of the most commonly used techniques in the mineral exploration field [1][2][3][4][5]. These surveys are often conducted using helicopters and airplanes with hanging transmitters and receivers attached. These existing exploration methods can quickly cover large areas. However, there is a safety issue when collecting high-quality data while flying close to the ground. Extensive research on exploration systems based on unmanned aerial vehicles (UAVs) has accelerated in recent years [6][7][8][9][10][11][12][13], in an effort to supplement the limitations of conventional airborne and ground exploration methods. In particular, UAV-based magnetic exploration has been of particular interest, and commercial

Geological Background
The research area for the current study, the Pocheon area, includes the intrusive, igneous body of Mount Gonam. Mount Gonam is located at the boundary between the southeastern part of the Imjingang Belt and the northwestern part of the Gyeonggi Massif (marked with a star in Figure 1a). The Gyeonggi gneiss complex is composed of schist and gneiss. The Paleoproterozoic Gyeonggi gneiss constitutes the bedrock of the study area. The Paleoproterozoic gneiss complex was intruded by Neoproterozoic plutonic rocks, Jurassic Daebo granite, and Cretaceous granite. The intrusive body in Mount Gonam occupies an area of approximately 3 km × 2 km in a north-south direction. Most of the deposits are in contact with the Paleoproterozoic gneiss complex, but in the northeastern part they are in contact with Cretaceous Myeongseongsan granite. According to the surface survey [23], the eastern side of the intrusive rock is melanocratic monzogabbro-monzodiorite, and the western side is leucocratic quartz monzodiorite (Figure 1b [24]). Yellow rectangle indicates the study area.
The deposits in the study area are iron ore that developed within gabbroid rocks intruding into the Precambrian metamorphic sedimentary rock. According to petrochemical studies [25,26], the genesis of the magnetite deposit is interpreted as an orthomagmatic deposit originating from basic magma. The ore production states of mines in this area can be classified into veins, lenticulars, masses, skarn-type ores, and stratified ore bodies, which are closely related to gabbroid rock. In a report presented in 2019 from the Gwanin mine with the location indicated in Figure 1b, the main ore body of the mine is classified as a skarn-type ore.

Aeromagnetic Exploration with a Helicopter
Since 1981, the Korea Institute of Geoscience and Mineral Resources (KIGAM) has been researching aeromagnetic exploration and developing technology to this end, in an attempt to identify magnetic anomalies throughout the Korean Peninsula. There has been great interest among the research community in using high-resolution aerial exploration to identify mineral resources. Figure 2 presents a helicopter-based KIGAM-manned aeromagnetic system. Three booms are attached to the lower part of the helicopter, and a total of four magnetic sensors are installed within the booms. A fluxgate magnetometer is attached to the center of the front boom to compensate for the magnetic field associated with the movement of the helicopter during flight. The other three magnetometers (G-822A, GeoMetrics Inc., San Jose, CA, USA) are scalar magnetometers that are used to measure the total magnetic intensity. One magnetometer is attached to the end of the boom facing the front of the helicopter, and the others are attached to either end of the boom installed in a direction perpendicular to that of the central boom. The three magnetometers have a sampling frequency of 10 Hz. In the helicopter, a DAARC500 data acquisition system (RMS Instruments, Mississauga, ON, Canada) is installed; this system is capable of calibrating against the magnetic response that occurs from the helicopter itself in real time. Additionally, the The deposits in the study area are iron ore that developed within gabbroid rocks intruding into the Precambrian metamorphic sedimentary rock. According to petrochemical studies [25,26], the genesis of the magnetite deposit is interpreted as an orthomagmatic deposit originating from basic magma. The ore production states of mines in this area can be classified into veins, lenticulars, masses, skarn-type ores, and stratified ore bodies, which are closely related to gabbroid rock. In a report presented in 2019 from the Gwanin mine with the location indicated in Figure 1b, the main ore body of the mine is classified as a skarn-type ore.

Aeromagnetic Exploration with a Helicopter
Since 1981, the Korea Institute of Geoscience and Mineral Resources (KIGAM) has been researching aeromagnetic exploration and developing technology to this end, in an attempt to identify magnetic anomalies throughout the Korean Peninsula. There has been great interest among the research community in using high-resolution aerial exploration to identify mineral resources. Figure 2 presents a helicopter-based KIGAM-manned aeromagnetic system. Three booms are attached to the lower part of the helicopter, and a total of four magnetic sensors are installed within the booms. A fluxgate magnetometer is attached to the center of the front boom to compensate for the magnetic field associated with the movement of the helicopter during flight. The other three magnetometers (G-822A, GeoMetrics Inc., San Jose, CA, USA) are scalar magnetometers that are used to measure the total magnetic intensity. One magnetometer is attached to the end of the boom facing the front of the helicopter, and the others are attached to either end of the boom installed in a direction perpendicular to that of the central boom. The three magnetometers have a sampling frequency of 10 Hz. In the helicopter, a DAARC500 data acquisition system (RMS Instruments, Mississauga, ON, Canada) is installed; this system is capable of calibrating against the magnetic response that occurs from the helicopter itself in real time. Additionally, the DAARC500 system collects and stores data from all installed sensors, such as the global positioning system and barometer. All data can be used to match, calibrate, and analyze magnetic findings. inerals 2021, 11, 665 DAARC500 system collects and stores data from all installed sensors, such positioning system and barometer. All data can be used to match, calibrat magnetic findings. In 1975, the National Geological Survey (current KIGAM) conducted s ical and magnetic surveys of the area [23]. Production activities began in e and continue to this day. However, although it was possible to confirm the r netic exploration in some areas with large magnetic anomalies, most of the in general were insufficient for confirming the magnetic response of an entir ing mines. Therefore, in 2019, additional exploration using helicopters was In this study, the flight line spacing and altitude of the surveys were se 60-100 m, respectively, to acquire high-quality data ( Figure 3). Real-tim from navigation systems and the radar altimeter was used to control the hel the specified line spacing and altitude range. In 1975, the National Geological Survey (current KIGAM) conducted surface geological and magnetic surveys of the area [23]. Production activities began in earnest in 1982 and continue to this day. However, although it was possible to confirm the results of magnetic exploration in some areas with large magnetic anomalies, most of the magnetic data in general were insufficient for confirming the magnetic response of an entire area, including mines. Therefore, in 2019, additional exploration using helicopters was planned.
In this study, the flight line spacing and altitude of the surveys were set at 100 m and 60-100 m, respectively, to acquire high-quality data ( Figure 3). Real-time information from navigation systems and the radar altimeter was used to control the helicopter within the specified line spacing and altitude range. DAARC500 system collects and stores data from all installed sensors, such as the global positioning system and barometer. All data can be used to match, calibrate, and analyze magnetic findings. In 1975, the National Geological Survey (current KIGAM) conducted surface geological and magnetic surveys of the area [23]. Production activities began in earnest in 1982 and continue to this day. However, although it was possible to confirm the results of magnetic exploration in some areas with large magnetic anomalies, most of the magnetic data in general were insufficient for confirming the magnetic response of an entire area, including mines. Therefore, in 2019, additional exploration using helicopters was planned.
In this study, the flight line spacing and altitude of the surveys were set at 100 m and 60-100 m, respectively, to acquire high-quality data ( Figure 3). Real-time information from navigation systems and the radar altimeter was used to control the helicopter within the specified line spacing and altitude range.

Aeromagnetic Exploration with a Multicopter
In this study, to obtain a more precise reading of the strong magnetic anomalies, we also conducted a multicopter aeromagnetic exploration survey. The multicopter was outfitted with a MagArrow scalar magnetometer (GeoMetrix Inc., Scottsdale, AZ, USA). The MagArrow system weighs only 1 kg and has a sampling rate of 1000 Hz, meaning that data can be acquired every centimeter when flying at a speed of 10 m/s [27]. In this study, the data acquired at 1000 Hz were downsampled to 50 Hz considering the flight speed and data redundancy. Among the various multicopters, the M210RTK (DJI, Shenzhen, China) was selected for this study, based on its maximum payload (<1.23 kg) and reliability. Figure 4 presents the multicopter configuration for this study. Helicopter and multirotor UAVs typically use a single magnetic sensor rigidly mounted on a platform or suspended below the platform [13]. Depending on the mounting method, magnetic signals generated by various UAVs and their on-board components have been investigated, and studies on optimal sensor placement have been conducted [28][29][30][31]. In these papers, it was concluded that the magnetic signal interference between the magnetic sensor and multicopter is minimized when the magnetic sensor is more than 3 m away from the multicopter [28][29][30][31][32]. Based on these studies, we hung the MagArrow 3 m below the multicopter. During the test flight, the magnetic field contributions from the multicopter were not observed in the obtained data and no additional calibration process was required to eliminate the magnetic signal interference. The maximum flight time with the MagArrow was about 20 min. To acquire high-quality data, the flight line spacing and altitude of the surveys were set at 50 m and 60 m, respectively, as shown in Figure 5. In this study, a single drone was used, and it was difficult to secure a line-of-sight distance due to rugged, mountainous terrain covered with trees. Flight data were collected 15 times while moving between the take-off and landing points. It took a total of about 8 h to collect data for an area of 25 km 2 .
In this study, to obtain a more precise reading of the strong magnetic also conducted a multicopter aeromagnetic exploration survey. The multic fitted with a MagArrow scalar magnetometer (GeoMetrix Inc., Scottsdale, A MagArrow system weighs only 1 kg and has a sampling rate of 1000 Hz, data can be acquired every centimeter when flying at a speed of 10 m/s [27] the data acquired at 1000 Hz were downsampled to 50 Hz considering th and data redundancy. Among the various multicopters, the M210RTK (D China) was selected for this study, based on its maximum payload (<1.23 kg ity. Figure 4 presents the multicopter configuration for this study. Helicop rotor UAVs typically use a single magnetic sensor rigidly mounted on a pl pended below the platform [13]. Depending on the mounting method, ma generated by various UAVs and their on-board components have been inv studies on optimal sensor placement have been conducted [28][29][30][31]. In these concluded that the magnetic signal interference between the magnetic sen copter is minimized when the magnetic sensor is more than 3 m away from ter [28][29][30][31][32]. Based on these studies, we hung the MagArrow 3 m below th During the test flight, the magnetic field contributions from the multicopte served in the obtained data and no additional calibration process was req nate the magnetic signal interference. The maximum flight time with the M about 20 min. To acquire high-quality data, the flight line spacing and altit veys were set at 50 m and 60 m, respectively, as shown in Figure 5. In this drone was used, and it was difficult to secure a line-of-sight distance d mountainous terrain covered with trees. Flight data were collected 15 times between the take-off and landing points. It took a total of about 8 h to coll area of 25 km 2 . To increase the amplitude of the sensed magnetic anomalies within the is important to fly as low as possible. However, because the target area to a rugged mountainous terrain covered with trees taller than 20 m, as not dangerous and difficult to perform low-altitude UAV surveying in practice UAV surveying requires advanced UAV operation skills, reliable equipmen ized software for mission planning to safely collect high-quality data while the ground. To increase the amplitude of the sensed magnetic anomalies within the subsurface, it is important to fly as low as possible. However, because the target area to be explored is a rugged mountainous terrain covered with trees taller than 20 m, as noted above, it is dangerous and difficult to perform low-altitude UAV surveying in practice. Lowaltitude UAV surveying requires advanced UAV operation skills, reliable equipment, and specialized software for mission planning to safely collect high-quality data while flying close to the ground. In this study, survey missions were planned and flown in automatic mode with UgCS software version 4.3.82 (SPH Engineering, Riga, Latvia). The altitude of the multicopter was set at a constant altitude based on the digital elevation model (DEM) provided by Google (Mountain View, CA, USA). Because the DEM information cannot reflect changes in vegetation such as the growth of trees over time, the flight altitude was set to 60 m for stable flight. Additionally, to avoid the risk of a flight accident when the multicopter disappeared from view or when the communication with the ground monitoring system was cut off, the survey was conducted by moving the take-off and landing sites as needed. In this study, survey missions were planned and flown in automatic mode with UgCS software version 4.3.82 (SPH Engineering, Riga, Latvia). The altitude of the multicopter was set at a constant altitude based on the digital elevation model (DEM) provided by Google (Mountain View, CA, USA). Because the DEM information cannot reflect changes in vegetation such as the growth of trees over time, the flight altitude was set to 60 m for stable flight. Additionally, to avoid the risk of a flight accident when the multicopter disappeared from view or when the communication with the ground monitoring system was cut off, the survey was conducted by moving the take-off and landing sites as needed.

Interpretation Methods of Aeromagnetic Data
The goal of interpretation is to determine the target sources based on measured aeromagnetic data. Inversion is a well-established tool that supports the interpretation of geophysical data by generating the subsurface distribution of a physical property from observed data. The obtained magnetic data must undergo a basic correction process before inversion is performed. Preprocessing for inversion of the aeromagnetic data in this study was as follows. First, the diurnal variation correction was performed based on the data obtained at the Cheongyang geomagnetic observatory in Korea. The correction was carried out to separate sudden temporal variation in the Earth's magnetic field, caused by the rotation of the Earth and solar wind from geographical variations. International Geomagnetic Reference Field (IGRF) correction was then applied to correct for secular variation. Given that the published IGRF coefficients represent the field originating from within the core, subtracting the IGRF field from the Earth's field should only leave anomalies caused by sources within the Earth's crust [14]. Lastly, interpolation and resampling were applied considering the inversion block size.

Interpretation Methods of Aeromagnetic Data
The goal of interpretation is to determine the target sources based on measured aeromagnetic data. Inversion is a well-established tool that supports the interpretation of geophysical data by generating the subsurface distribution of a physical property from observed data. The obtained magnetic data must undergo a basic correction process before inversion is performed. Preprocessing for inversion of the aeromagnetic data in this study was as follows. First, the diurnal variation correction was performed based on the data obtained at the Cheongyang geomagnetic observatory in Korea. The correction was carried out to separate sudden temporal variation in the Earth's magnetic field, caused by the rotation of the Earth and solar wind from geographical variations. International Geomagnetic Reference Field (IGRF) correction was then applied to correct for secular variation. Given that the published IGRF coefficients represent the field originating from within the core, subtracting the IGRF field from the Earth's field should only leave anomalies caused by sources within the Earth's crust [14]. Lastly, interpolation and resampling were applied considering the inversion block size.
We used VPmg, which is a magnetic 3D modeling and inversion program developed by Fullagar Geophysics Pty Ltd. (Blackmans Bay, Tasmania, Australia) [33][34][35][36]. VPmg permits a wide variety of inversion styles to accommodate homogeneous unit properties, geometry, VP basement properties, and heterogeneous properties. In this research, due to a lack of known information about the subsurface, heterogeneous property inversion was selected without constraints. We first established that the starting model was a two-layer model. The first layer extended from the topography to a flat subsurface at 0 m of elevation, and was split into a heterogeneous layer. The 1:5000-scale digital topographic map obtained from the National Geographic Information Institute was used as the topography data. The elevation of the flat subsurface was determined considering the depth of the iron ore deposit already confirmed. The inversion cell size of the layer was assigned based on the line spacing. A scaling factor of 1.2% was applied to the vertical cell size with increasing depth. The initial susceptibility was set to 0 SI (Magnetic susceptibility). The magnetic field parameters used for inversion were an inclination of 55.059 • , a declination of −8.5731 • , and an amplitude of 51,427.8 nT. One important problem in potential field inversion is the lack of depth resolution. To enhance the depth resolution, depth-based weighting was used. The detailed inversion results are presented in the next section.

Results and Discussion
This section presents the inversion results obtained using the aeromagnetic data introduced above and the inversion technique. Figure 6 shows the inversion results of the aeromagnetic survey using a helicopter system. Figure 6a Figure 6a. The cross section shows the predicted distribution of potential iron ore (pink region) based on the position information of iron ore confirmed through drilling in the mine and the region developed so far. Area B is the area where iron ore is exposed on the surface as the pit that had been developed in the past collapsed. In the inversion results, the high susceptibility distribution is consistent with area B. In the case of area A, contrary to expectations, the anomalous body present directly under the surface and the susceptibility is relatively smaller than that of area B.
For a more detailed investigation, as described above, we also carried out an aeromagnetic exploration using a multicopter. Figure 7 shows the inversion result of the aeromagnetic survey using the multicopter system. Figure 7a displays the plan view of the inversion result superimposed on the geological map of Gwanin mine. Figure 7b,c displays the 3D inversion result superimposed on the geological cross section interpreted in Gwanin mine. Similar to the previous inversion results, Figure 7b,c only displays the region of the 3D inverted model with a susceptibility greater than 0.28 and 0.25 SI, respectively. Considering the line spacing (50 m), the size of inversion block was chosen (50 m × 50 m × 50 m). Compared with the previous inversion results using the aeromagnetic data from the helicopter, the location of a highly anomalous body (area B) was identified similarly. Additionally, a newly discovered anomalous body in area C was confirmed. However, the inversion result for area A, where there is an active mine, was still different from the location of magnetite that had already been identified (pink region in area A in Figure 7c). Therefore, to analyze the cause of the inconsistency in the inversion result in area A, the magnetic field data was simulated using the model constructed based on geological cross sections. Even when there was an anomaly having a higher magnetic susceptibility than the actual ore body (0.81 SI) and having a larger size, the magnetic intensity was close to zero at the investigation altitude. Based on these results, in the case of area A, we determined that there is a possibility that another anomalous body exists in the upper part and that through the surface geological survey, we could find the magnetite on the surface.
By comprehensively analyzing the inversion using the two sets of aeromagnetic data and the results of the geological analysis, two directional drilling and one vertical drilling project were designed (indicated in Figures 6a and 7a) for verification. In the case of two-directional drilling, to obtain drilling data passing through area C, drilling was carried out at a location as close to area C as possible using the previously used pit. Well 1 was drilled to 201.2 m. The strike was N75E, and the slope was 10 • . Well 2 was drilled to 180.5 m, with a strike of N80E and a slope of 10 • . Well 3 (vertical drilling) was drilled to 271 m. Figure 8 shows the facies classification profiles of three wells. The results are shown based on the measured depth without considering the strike and slope. For clear distinction, the area that reacts with the magnet is marked with color, and the other areas are marked with achromatic color. In the case of well 1, iron-containing minerals appeared from 92 to 201.2 m. In well 2, iron-containing minerals appeared from 136 to 180.5 m. Comparing the two sets of drilling data, well 1 showed a relatively strong magnetic response; the strongest magnetic spot was located in the vicinity of 100 m of excavation depth. In the case of well 3, the weak magnetic response appeared at the vicinity of 250 m. troduced above and the inversion technique. Figure 6 shows the inversion results of the aeromagnetic survey using a helicopter system. Figure 6a displays the inversion results superimposed on the geological map of Gwanin mine, and Figure 6b,c displays the 3D inversion results superimposed on the geological cross section interpreted from Gwanin mine. Figure 6a displays the region of the 3D inverted model with a susceptibility greater than 0.2 SI. Figure 6b,c only displays the region of the 3D inverted model with a susceptibility greater than 0.28 SI and 0.25 SI, respectively. Considering the line spacing (100 m), the size of inversion block was chosen (100 m × 100 m × 100 m). The location of the geological cross section shown in Figure 6b,c is indicated by dashed lines DE and EF in Figure  6a. The cross section shows the predicted distribution of potential iron ore (pink region) based on the position information of iron ore confirmed through drilling in the mine and the region developed so far. Area B is the area where iron ore is exposed on the surface as the pit that had been developed in the past collapsed. In the inversion results, the high susceptibility distribution is consistent with area B. In the case of area A, contrary to expectations, the anomalous body present directly under the surface and the susceptibility is relatively smaller than that of area B. For a more detailed investigation, as described above, we also carried out an aeromagnetic exploration using a multicopter. Figure 7 shows the inversion result of the aeromagnetic survey using the multicopter system. Figure 7a displays the plan view of the inversion result superimposed on the geological map of Gwanin mine. Figure 7b,c displays the 3D inversion result superimposed on the geological cross section interpreted in Gwanin mine. Similar to the previous inversion results, Figure 7b,c only displays the region of the 3D inverted model with a susceptibility greater than 0.28 and 0.25 SI, respectively. Considering the line spacing (50 m), the size of inversion block was chosen (50 m × 50 m × 50 m). Compared with the previous inversion results using the aeromagnetic data from the helicopter, the location of a highly anomalous body (area B) was identified similarly. Additionally, a newly discovered anomalous body in area C was confirmed. However, the inversion result for area A, where there is an active mine, was   Figure 6a. The cross section shows the predicted distribution of potential iron ore (pink region) based on the position information of iron ore confirmed through drilling in the mine and the region developed so far.
Additionally, to interpret the interior of the 3D magnetic inversion results quantitatively, we used the measured susceptibility results obtained from laboratory experiments of a core sample obtained from the main active mine (area A) and three new boreholes. Figure 9 shows box plots of the ranges of measured susceptibility according to rock type. In the case of area A, where iron ore development is in progress, it shows a relatively strong susceptibility distribution in the range of 0.78 to 0.84 SI. Among the drilling data of area C, the SI range of 0.3 to 0.4 for rock samples is classified as relatively high-grade ore bodies (red in Figure 8). The SI range of 0.11 to 0.19 for rock samples is classified as low-grade ore bodies (orange area in Figure 8). In the case of well 3, the region that reacted with the magnet shows an SI distribution of 0.013 to 0.16. Additionally, to interpret the interior of the 3D magnetic inversion results quantitatively, we used the measured susceptibility results obtained from laboratory experiments of a core sample obtained from the main active mine (area A) and three new boreholes. Figure 9 shows box plots of the ranges of measured susceptibility according to rock type. In the case of area A, where iron ore development is in progress, it shows a relatively strong susceptibility distribution in the range of 0.78 to 0.84 SI. Among the drilling data of area C, the SI range of 0.3 to 0.4 for rock samples is classified as relatively high-grade ore bodies (red in Figure 8). The SI range of 0.11 to 0.19 for rock samples is classified as low-grade ore bodies (orange area in Figure 8). In the case of well 3, the region that reacted with the magnet shows an SI distribution of 0.013 to 0.16.
By performing a comparative analysis by superimposing the drilling data on the inversion result of the aeromagnetic data acquired using the multicopter system (Figure 10), well 1 was designed to penetrate the inside of the anomalous body in area C, and well 2 was configured to penetrate the boundary of the anomalous body. This is consistent with the observation that the intensity of the magnetic strength decreased from well 1 to well 2. Additionally, the appearance of iron-containing amphibolite after 92 m excavation in well 1 was consistent with the analysis results that area C was separated from area B. Additionally, when only values greater than 0.2 SI are shown among the 3D susceptibility By performing a comparative analysis by superimposing the drilling data on the inversion result of the aeromagnetic data acquired using the multicopter system (Figure 10), well 1 was designed to penetrate the inside of the anomalous body in area C, and well 2 was configured to penetrate the boundary of the anomalous body. This is consistent with the observation that the intensity of the magnetic strength decreased from well 1 to well 2. Additionally, the appearance of iron-containing amphibolite after 92 m excavation in well 1 was consistent with the analysis results that area C was separated from area B. Additionally, when only values greater than 0.2 SI are shown among the 3D susceptibility models obtained through the inversion, the depth of the anomaly decreases from region B toward the drilling position of Well 3. The inversion results are consistent with the tendency to have susceptibility less than 0.2 SI from a depth of 250 m in well 3.       As a result of analyzing the distribution of ore bodies by applying two aeromagnetic methods, it was difficult for the manned aerial survey using helicopter to fly at low altitude along narrow line spacing due to safety concerns. On the other hand, in the case of using a multicopter, because automatic flight was performed based on sensors, it was possible to obtain data of narrow line spacing and relatively close to the ground. For this reason, the area that appeared to be a single ore body in the inversion result of the aeromagnetic data obtained using a helicopter was clearly distinguished in the inversion result of the aeromagnetic data obtained using the multicopter.
However, since UAV multicopters need to fly within line-of-sight, data must be collected while moving to the take-off and landing points for safety in the mountainous terrain, which can be an inefficient method to acquire a wide range of data. Therefore, it is recommended that these two methods are used as complements to each other for studying the distribution of ore bodies in complex mountainous terrain covered with trees.
Minerals 2021, 11, 665 Figure 10. Drilling results shown overlaid on the inversion result of the aeromagnetic data quired using the multicopter system.
As a result of analyzing the distribution of ore bodies by applying two aerom methods, it was difficult for the manned aerial survey using helicopter to fly at tude along narrow line spacing due to safety concerns. On the other hand, in th using a multicopter, because automatic flight was performed based on sensors, it w sible to obtain data of narrow line spacing and relatively close to the ground. reason, the area that appeared to be a single ore body in the inversion result of t magnetic data obtained using a helicopter was clearly distinguished in the inversio of the aeromagnetic data obtained using the multicopter.
However, since UAV multicopters need to fly within line-of-sight, data mus lected while moving to the take-off and landing points for safety in the mountain rain, which can be an inefficient method to acquire a wide range of data. Theref recommended that these two methods are used as complements to each other for s the distribution of ore bodies in complex mountainous terrain covered with trees

Conclusions
We describe a case study in which we identified the distribution of iron ore in the west-central portion of the Korean Peninsula using aeromagnetic techniq obtained regional information based on manned aeromagnetic data from a helico vey. Next, we acquired magnetic data with higher resolution from a UAV mu focusing on the target area. The advantage of using a helicopter was that it could data quickly over a wide area. The advantage of using the multicopter was that perform precise automatic flight control, so it was suitable for acquiring high-re magnetic data. Through the inversion result of the higher-resolution magnetic

Conclusions
We describe a case study in which we identified the distribution of iron ore mineral in the west-central portion of the Korean Peninsula using aeromagnetic techniques. We obtained regional information based on manned aeromagnetic data from a helicopter survey. Next, we acquired magnetic data with higher resolution from a UAV multicopter focusing on the target area. The advantage of using a helicopter was that it could acquire data quickly over a wide area. The advantage of using the multicopter was that it could perform precise automatic flight control, so it was suitable for acquiring high-resolution magnetic data. Through the inversion result of the higher-resolution magnetic data, we confirmed the possibility of the existence of a new ore body. Finally, drilling was performed based on the analysis results, and new ore bodies were successfully identified.