• Journal of the Korean earth science society
• /
• v.32 no.6
• /
• pp.548-559
• /
• 2011

3-D Multi-geophysical surveys were carried out around the Hwasan caldera at the Euisung Sub-basin. To overcome the limitations of resolutions in previous studies, dense gravity data and magnetotelluric (MT) data were obtained and analyzed. In this study, the independent inversion models from gravity and MT data were integrated using correlation and classification approaches for 3-D imaging of the geologic structures. A Structure Index (SI) method was proposed and applied to the integration and classification analyses. This method consists of Type Angle (TA) and Type Intensity (TI) values, which are estimated by the spatial correlation and abnormality of the physical properties. The SI method allowed the classification analysis to be effectively performed. Major findings are as follows: 1) pyroclastic rocks around the central area of the Hwasan caldera with lower density and resistivity than those of neighboring regions extended to a depth of around 1 km, 2) intrusive igneous rocks with high resistivity and density were imaged around the ring fault boundary, and 3) a basement structure with low resistivity and high density, at a depth of 3-5 km, was inferred by the SI analysis.

• Geophysics and Geophysical Exploration
• /
• v.15 no.4
• /
• pp.177-189
• /
• 2012

Evaluation of three-dimensional ore body modeling has been performed by applying the geostatistical integration technique to multiple geophysical (electrical resistivity, MT) and geological (borehole data, physical properties of core) information. It was available to analyze the resistivity range in borehole and other area through multiple geophysical data. A correlation between resistivity and density from physical properties test of core was also analyzed. In the case study results, the resistivity value of ore body is decreased contrast to increase of the density, which seems to be related to a reason that the ore body (magnetite) includes heavy conductive component (Fe) in itself. Based on the lab test of physical properties in iron mine region, various geophysical, geological and borehole data were used to provide ore body modeling, that is electrical resistivity, MT, physical properties data, borehole data and grade data obtained from borehole data. Of the various geostatistical techniques for the integrated data analysis, in this study, the SGS (sequential Gaussian simulation) method was applied to describe the varying non-homogeneity depending on region through the realization that maintains the mean and variance. With the geostatistical simulation results of geophysical, geological and grade data, the location of residual ore body and ore body which is previously reported was confirmed. In addition, another highly probable region of iron ore bodies was estimated deeper depth in study area through integrated modeling.

• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.361-366
• /
• 2009

Layered-earth Green's functions in electormagnetic (EM) surveys play a key role in modeling the response of exploration targets. They are computed through the Hankel transforms of analytic kernels. Computational precision depends upon the choice of algebraically equivalent forms by which these kemels are expressed. Since three-dimensional (3D) modeling can require a huge number of Green's function evaluations, total computational time can be influenced by computational time for the Hankel transform evaluations. Linear digital filters have proven to be a fast and accurate method of computing these Hankel transforms. In EM modeling for 3D inversion, electric fields are generally evaluated by the secondary field formulation to avoid the singularity problem. In this study, three components of electric fields for five different sources on the surface of homogeneous half-space were derived as primary field solutions. Moreover, reflection coefficients in TE and TM modes were produced to calculate EM responses accurately for a two-layered model having a sea layer. Accurate primary fields should substantially improve accuracy and decrease computation times for Green's function-based problems like MT problems and marine EM surveys.

• Journal of the Korean Geophysical Society
• /
• v.5 no.4
• /
• pp.305-319
• /
• 2002

The three-dimensional MT(Magneto-Telluric) modeling is performed to examine the validity of difference arrow of GDS(Geomagnetic Depth Sounding) survey, In this paper, we investigate the validity of the difference arrow on three configurations of conductors; which is located 1) at surface, 2) at the deep part and 3) vertically extended f개m surface to the deep part, respectively, For conductors located at surface, the validity of difference arrows is certified in our numerical model when long periods over 40 minutes are used or the distance between sea and conductor is over 150 km. However, for conductors located at the deep part, the validity of difference arrow is dependent on the size of conductors. Further, if the size of conductor is adequately larger than that of our model, we recognize the possibility that the mutual coupling of them influences up to longer periods, Moreover, in case of conductors which is vertically extended from surface to the deer part, the mutual coupling of them is reinforced for all periods, especially for longer periods, so that the validity of difference arrow is considerably in doubt. Therefore, to remove the known conductor effect such as the sea effect from the observed induction arrow, the mutual coupling between them must be examined. The difference arrow that certifies the validity in this way can only provide the Subsurface information based on physical supports.

• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.145-155
• /
• 2005

A two-dimensional (2-D) interpretation of MT data has been performed for the purpose of fracture detection for geothermal development. Remote stations have been operated in Kyushu, Japan (480 km apart) as well as in Korea (60 km and 165 km apart in 2002 and 2003 data set, respectively). Apparent resistivity and phase curves calculated by remote processing with the Japan remote data showed enough quality for 2-D inversion for the whole frequency range. Remote reference processing with Korea remote reference data also showed quite good continuity in apparent resistivity and phase curves except some noisy frequency bands; around the power frequency, 60 Hz, and around the dead band $10^{-1}Hz\;Hz\;\~1\;Hz$, where the natural EM signal is known to be very weak. Even though the subsurface showed severe three-dimensional (3-D) characteristics in the survey area so that 2-D inversion by itself could not give enough information for deep geological structures, the 2-D inversion for the 5 survey lines showed several common features. The conductive semi-consolidate mudstone layer is dipping from north to south (about 500 m depth on the south and 200 m on the north most part of the survey area). The boundary between the low (L-2) and high (H-2) resistivity anomalies can be thought as a major fault with strike $N15^{\circ}E$, passing through the sites 206, 112 and 414. The shallow (< 1 km) conductive anomalies (L-4) seem to be fracture zones having strike E-W (at site 105) and $N60^{\circ}W$ (at site 434). And there exists a conductive layer in the western and west-southern part of the survey area in the depth below $2\~3\;km$, for which further investigation is to be needed.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.147-153
• /
• 2007

The conjugate gradient (CG) method is one of the most efficient algorithms for solving a linear system of equations. In addition to being used as a linear equation solver, it can be applied to a least-squares problem. When the CG method is applied to large-scale three-dimensional inversion of magnetotelluric data, two approaches have been pursued; one is the linear CG inversion in which each step of the Gauss-Newton iteration is incompletely solved using a truncated CG technique, and the other is referred to as the nonlinear CG inversion in which CG is directly applied to the minimization of objective functional for a nonlinear inverse problem. In each procedure we only need to compute the effect of the sensitivity matrix or its transpose multiplying an arbitrary vector, significantly reducing the computational requirements needed to do large-scale inversion.

• Journal of the Korean earth science society
• /
• v.30 no.2
• /
• pp.153-164
• /
• 2009

We examined the regional 1-D deep resistivity structure of the Korean Peninsula using MT data acquired at seven sites located in the Kyongsang Basin and Kyonggi Massif. At the sites located in the Kyongsang Basin, surrounding sea distorts observed MT response and hence this distortion, so called "sea effect", is corrected using an iterative tensor stripping method. The 1-D layered inversion results for the seven MT sites reveal 4 layered structure, which is composed of 1) near surface layer, 2) upper crust, 3) lower crust and upper mantle, and 4) asthenosphere from the surface downward. Conrad interface, which is a boundary between upper and lower crust, is distinctly identified beneath all the MT sites. Conrad interface depth is estimated to about be 17km in the Kyongsang Basin and about 12km in the Kyonggi Massif, while the upper crust of the Kyongsang Basin is about 5 times more resistive than that of the Kyonggi Massif. Finally, asthenosphere is inferred to exist below a depth of approximately 100km with a resistivity of 200-300 ohm-m.

• Economic and Environmental Geology
• /
• v.50 no.6
• /
• pp.487-495
• /
• 2017

A granite rock body, called 'Bukbawi', located on a mountaineering trail at Mt. Palgong Provincial Park is popular among the public because it resembles a percussion instrument. If someone hits the specific surface area of this rock body, people can hear drum-like sound. Such phenomenon may be geologically associated with exfoliation process of the granite body or miarolitic cavity developed after gasses escaped during formation of granite. To understand better the inner structure causing drum-like sound, we carried out a non-destructive ground-penetrating radar survey. In this study, as our primary target is very close to the surface, we utilized 1 GHz antennas to produce high-resolution near-surface images. In order to construct 3-D internal images, the measurements were conducted along a pre-defined grid. The processed radargrams revealed that the locations associated with 'drum' sound coincide with strong reflections. In addition, both reflection patterns of fracture and cavity were observed. To further quantify the observed reflections, we simulated GPR scans from a synthetic fracture in a granite body, filled with different materials. The simulated results suggest that both exfoliation process and miarolitic cavity may have contributed to the 'drum' phenomena. Furthermore, the radargrams showed a well-developed cavity signature where two major reflection planes were crossed. Thus, our study is an example of non-destructive geophysical studies that can promote Earth Science in the broader community by examining geological structures attracting the public.

• Journal of the Korean earth science society
• /
• v.28 no.2
• /
• pp.187-201
• /
• 2007

The first magnetotelluric (MT) transect across the Korean Peninsula was obtained traversing from the East Sea shoreline to the Yellow Sea shoreline. The MT survey profile was designed perpendicular to the strike of the principal geologic structure of the Korean Peninsula $(N30^{\circ}E)$, so-called 'China direction'. MT data were achieved at 50 sites with spacings of $3{\sim}8km$ along the 240 km survey line. The impedance responses are divided into four subsets reflecting typical geological units: the Kyonggi Massif, the Okchon Belt, the western part of the Kyongsang Basin, and the eastern part of the Kyongsang Basin. In the western part of the Kyongsang Basin, the thickness of the sedimentary layer is estimated to be about 3 km to 8 km and its resistivity is a few hundred ohm-m. A highly conductive layer with a resistivity of 1 to 30 ohm-m was detected beneath the sedimentary layer. The MT data at the Okchon Belt show peculiar responses with phase exceeding $90^{\circ}$. This feature may be explained by an electrically anisotropic structure which is composed of a narrow anisotropic block and an anisotropic layer. The Kyonggi Massif and the eastern part of Kyongsang Basin play a role of window to the deep geoelectrical structure because of the very high resistivity of upper crust. The second layers with highest resistivities in 1-D conductivity models occupy the upper crust with thicknesses of 13 km in the Kyonggi Massif and 18 km in the eastern Kyongsang Basin, respectively.

• Korean Journal of Remote Sensing
• /
• v.36 no.5_4
• /
• pp.1267-1276
• /
• 2020

To quantitatively predict the impacts of large-scale volcanic eruptions of Mt. Baekdu on air quality and damage around the Korean Peninsula, a three-dimensional chemistry-transport modeling system (Weather Research & Forecasting - Sparse Matrix Operation Kernel Emission - Comunity Multi-scale Air Quality) was adopted. A worst-case meteorology scenario was selected to estimate the direct impact on Korea. This study applied the typical worst-case scenarios that are likely to cause significant damage to Korea among worst-case volcanic eruptions of Mt. Baekdu in the past decade (2005~2014) and assumed a massive VEI 4 volcanic eruption on May 16, 2012, to analyze the concentration of PM_2.5 caused by the volcanic eruption. The effects of air quality in each region-cities, counties, boroughs-were estimated, and vulnerable areas were derived by conducting an exposure assessment reflecting vulnerable groups. Moreover, the effects of cities, counties, and boroughs were analyzed with a high-resolution scale (9 km × 9 km) to derive vulnerable areas within the regions. As a result of analyzing the typical worst-case volcanic eruptions of Mt. Baekdu, a discrepancy was shown in areas between high PM_2.5 concentration, high population density, and where vulnerable groups are concentrated. From the result, PM_2.5 peak concentration was about 24,547 ㎍/㎥, which is estimated to be a more serious situation than the eruption of Mt. St. Helensin 1980, which is known for 540 million tons of volcanic ash. Paju, Gimpo, Goyang, Ganghwa, Sancheong, Hadong showed to have a high PM_2.5 concentration. Paju appeared to be the most vulnerable area from the exposure assessment. While areas estimated with a high concentration of air pollutants are important, it is also necessary to develop plans and measures considering densely populated areas or areas with high concentrations of susceptible population or vulnerable groups. Also, establishing measures for each vulnerable area by selecting high concentration areas within cities, counties, and boroughs rather than establishing uniform measures for all regions is needed. This study will provide the foundation for developing the standards for disaster declaration and preemptive response systems for volcanic eruptions.