• Journal of the Korean Geophysical Society
• /
• v.8 no.1
• /
• pp.15-21
• /
• 2005

The apparent source spectrum of the Fukuoka earthquake is estimated at the seismic basement by removing from the observed spectra at Korean seismic stations the path and site responses that were previously revealed through inversion process applied to large spectral D/B accumulated until 2004. The approximate source spectrum is also estimated by using data recorded near the epicenter from various Japanese seismic networks and compared with the Korean source spectrum. The comparison result shows that there is good agreement among source spectra estimated based on the data from seismic networks of Korea at large distances (190 km

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.85-92
• /
• 2008

Magnetotelluric (MT) methods are widely applied as an effective exploration technique to geothermal surveys. Two-dimensional (2-D) analysis is frequently used to investigate a complicated subsurface structure in a geothermal region. A 2-D finite-element method (FEM) is usually applied to the MT analysis, but we must pay attention to the accuracy of so-called auxiliary fields. Rodi (1976) proposed an algorithm of improving the accuracy of auxiliary fields, and named it as the MOM method. Because it introduces zeros into the diagonal elements of coefficient matrix of the FEM total equation, a pivoting procedure applied to the symmetrical band matrix makes the numerical solution far less efficient. The MOM method was devised mainly for the inversion analysis, in which partial derivatives of both electric and magnetic fields with respect to model parameters are required. In the case of forward modeling, however, we do not have to resort to the MOM method; there is no need of modifying the coefficient matrix, and the auxiliary fields can be elicited from the regular FEM solution. The computational efficiency of the MOM method, however, can be greatly improved through a sophisticated rearrangement of the total equation.

• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.64-73
• /
• 2015

Recently, hydraulic fracturing is used in various fields and microseismic monitoring is one of the best methods for judging where hydraulic fractures exist and how they are developing. When locating microseismic events using single vertical well data, distances from the vertical array and depths from the surface are generally decided using time differences between compressional (P) wave and shear (S) wave arrivals and azimuths are calculated using P wave hodogram analysis. However, in field data, it is sometimes hard to acquire P wave data which has smaller amplitude than S wave because microseismic data often have very low signal to noise (S/N) ratio. To overcome this problem, in this study, we developed a grid search algorithm which can find event location using all combinations of arrival times recorded at receivers. In addition, we introduced and analyzed the method which calculates azimuths using S wave. The tests of synthetic data show the inversion method using all combinations of arrival times and receivers can locate events without considering the origin time even using only single phase. In addition, the method can locate events with higher accuracy and has lower sensitivity on first arrival picking errors than conventional method. The method which calculates azimuths using S wave can provide reliable results when the dip between event and receiver is relatively small. However, this method shows the limitation when dip is greater than about $20^{\circ}$ in our model test.

• Geophysics and Geophysical Exploration
• /
• v.15 no.1
• /
• pp.33-38
• /
• 2012

A resistivity method has been applied to wide range of engineering and environmental problems with the help of automatic and precise data acquisition. Thus, more accurate modeling and inversion of time-lapse monitoring data are required since resistivity monitoring has been introduced to quantitatively find out subsurface changes With respect to time. Here, we used the finite element method (FEM) for 3D resistivity modeling since the method is easy to realize complex topography and arbitrary shaped anomalous bodies. In the FEM, the linear elements, also referred to as first order elements, have certain advantages of simple formulation and narrow bandwidth of system equation. However, the linear elements show the poor accuracy and slow convergence of the solution with respect to the number of elements or nodes. To achieve the higher accuracy of finite element solution, high order elements are generally used. In this study, we developed a 3D resistivity modeling program using high order Serendipity elements. Comparing the Serendipity element solutions for a cube model with the linear element solutions, we assured that the Serendipity element solutions are more accurate than the linear element solutions in the 3D resistivity modeling.

• Economic and Environmental Geology
• /
• v.35 no.3
• /
• pp.273-284
• /
• 2002

The Gravity-Geologic Method (GGM) was implemented for bathymetric determinations in the Drake Passage, Antarctica, using global marine Free-air Gravity Anomalies (FAGA) data sets by Sandwell and Smith (1997) and local echo sounding measurements. Of the 6548 bathymetric sounding measurements, two thirds of these points were used as control depths, while the remaining values were used as checkpoints. A density contrast of 9.0 gm/㎤ was selected based on the checkpoints predictions with changes in the density contrast assumed between the seawater and ocean bottom topographic mass. Control depths from the echo soundings were used to determine regional gravity components that were removed from FAGA to estimate the gravity effects of the bathymetry. These gravity effects were converted to bathymetry by inversion. In particular, a selective merging technique was developed to effectively combine the echo sounding depths with the GGM bathymetiy to enhance high frequency components along the shipborne sounding tracklines. For the rugged bathymetry of the research area, the GGM bathymetry shows correlation coefficients (CC) of 0.91, 0.92, and 0.85 with local shipborne sounding by KORDI, GEODAS, and a global ETOPO5 model, respectively. The enhanced GGM by selective merging shows imploved CCs of 0.948 and 0.954 with GEODAS and Smith & Sandwell (1997)'s predictions with RMS differences of 449.8 and 441.3 meters. The global marine FAGA data sets and other bathymetric models ensure that the GGM can be used in conjunction with shipborne bathymetry from echo sounding to extend the coverage into the unmapped regions, which should generate better results than simply gridding the sparse data or relying upon lower resolution global data sets such as ETOPO5.

• Economic and Environmental Geology
• /
• v.31 no.2
• /
• pp.127-139
• /
• 1998

A new technique of simultaneous inversion for 3-D seismic velocity structure by using direct, reflected, and refracted waves is applied to the southeast part of the Korean Peninsula including Pohang Basin, Kyongsang Basin and Ryongnam Massif. Pg, Sg, PmP, SmS, Pn, and Sn arrival times of 44 events with 554 seismic rays are inverted for locations and crustal structure. $6{\times}6$ with $0.5^{\circ}$ and 8 layers (4 km each layer) model was inverted. 3-D seismic crustal velocity tomography including eight sections from surface to Moho, ten profiles along latitude and longitude are analyzed. The results are as follows: 1) the average velocity and thickness of sediment are 5.04 km/s and 3-4 km, and the velocity of basement is 6.11 km/s. The shape of velocity in shallower layer is agreement with Bouguer gravity anomaly (Cho et al., 1997). 2) the velocities fluctuate strongly in the upper crust. The velocity distribution of the lower crust under Conrad appears basically horizontal. 3) the average depth of Moho is 30.4 km, and velocity is 8.01 km/s. 4) from the velocity and depth of the sediment, the thickness, velocity and form of the upper crust, and the depth and form of Moho, we can find the obvious differences among Ryongnam Massif, Kyongsang Basin and Pohang Basin. 5) the deep faults (a Ulsan series faults) near Kyongju and Pohang areas can be found to be normal and/or thrust faults with detachment extended to the bottom of the upper crust.

• Economic and Environmental Geology
• /
• v.31 no.3
• /
• pp.235-247
• /
• 1998

A new technique of simultaneons inversion for 3-D seismic velocity structure by using direct, reflected, and refracted waves is applied to the center of the Korean Peninsula including Pyongnam Basin, Kyonggi Massif, Okchon Fold Zone, Taebaeksan Fold Zone, Ryongnam Massif and Kyongsang Basin. Pg, Sg, PmP, SmS, Pn, and Sn arrival times of 32 events with 404 seismic rays are inverted for locations and crustal structure. 5 ($1^{\circ}$ along the latitude)${\times}6$ ($0.5^{\circ}$ along the longitude) ${\times}8$ block (4 km each layer) model was inverted. 3-D seismic crustal velocity tomography including eight sections from the surface to the Moho, eight profiles along latitude and longitude and the Moho depth distribution was determined. The results are as follows: (1) the average velocity and thickness of sediment are 5.15 km/sec and 3-4 km, and the velocity of basement is 6.12 km/sec. (2) the velocities fluctuate strongly in the upper crust, and the velocity distribution of the lower crust under Conrad appears basically horizontal. (3) the average depth of Moho is 29.8 km and velocity is 7.97 km/sec. (4) from the sedimentary depth and velocity, basement thickness and velocity, form of the upper crust, the Moho depth and form of the remarkable crustal velocity differences among Pyongnam Basin, Kyonggi Massif, Okchon Zone, Ryongnam Massif and Kyongsang Basin can be found. (5) The different crustal features of ocean and continent crust are obvious. (6) Some deep index of the Chugaryong Rift Zone can be located from the cross section profiles. (7) We note that there are big anisotropy bodies near north of Seoul and Hongsung in the upper crust, implying that they may be related to the Chugaryong Rift Zone and deep fault systems.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.3C
• /
• pp.143-155
• /
• 2006

SPT-Uphole tomography method was introduced for the evaluation of near subsurface shear wave velocity (Vs) distribution map. In SPT-Uphole method, SPT (Standard Penetration Test) which is common in geotechnical site investigation was used as a source and several surface geophones in line were used as receivers. Vs distribution map which is the triangular shape around the boring point can be developed by tomography inversion. To obtain the exact travel time information of shear wave component, a procedure using the magnitude summation of vertical and horizontal components was used based on the evaluation of particle motion at the surface. It was verified that proposed method could give reliable Vs distribution map through the numerical study using the FEM (Finite Element Method) model. Finally, SPT-Uphole tomography method was performed at the weathered soil site where several boring data with SPT-N values are available, and the feasibility of proposed method was verified in the field.

• Economic and Environmental Geology
• /
• v.39 no.2 s.177
• /
• pp.129-149
• /
• 2006

Previously published stratigraphic, sedimentologic, paleontologic, paleomagnetic and geophysical data are reviewed to make an understanding on the tectonic evolution of the Cretaceous Gyeongsang (Kyongsang) basin, southeast Korea. A stratigraphic framework and a tectonic model on the formation and deformation of the Gyeongsang Basin are newly proposed on the basis of integration these data with magmatism and mineralization ages in the basin. A newly proposed stratigraphic framework indicates that strata in the basin can be subdivided into five distinct stratigraphic units that represent pre-rifting, syn-rifting, inversion I, II, and III stages. The Gyeongsang Basin was formed initially as a pre-rifting stage due to north-south extension in the Late Jurassic prior to a syn-riftins stage that resulted from east-west extension during the Early Cretaceous. In the Late Cretaceous, the basin was deformed by three-staged sequential deformation of north-south, northwest-southeast, and east-west compressions. The tectonic history of the basin has been largely controlled by the change of motion of the Izanagi Plate from north to northwest during the Cretaceous. In the early Cretaceous, the Izanagi Plate began to subduct northward beneath the Eurasian Plate and caused the left-lateral strike-slip fault systems in the southern part of the peninsula. The left-lateral wrenching of these fault systems was causally linked to development of pull-apart basins, such as the Gyeongsang Basin in the southeastern part of the peninsula. However, northwestward movement of the Izanagi Plate during the Late Cretaceous probably led to the extensive volcanism as well as sequential deformations in the basin. The stratigraphic and tectonic model, which is newly proposed as a result of this study, may be expected to enhancing the efficiency for exploration and exploitation of useful mineral resources in the basin as well as establishing geologic history in the Cretaceous Gyeongsang Basin. Together with the spatial and temporal correlation of the Cretaceous basins in adjacent areas, this stratigraphic and tectonic model provides a new geologic paradigm to delineate the sophisticated tectonic history of East Asia turing the Cretaceous.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.20 no.1
• /
• pp.73-87
• /
• 2018

A variety of micro meteorological variables such as air temperature, wind, solar radiation and latent heat at Gwangneung forests (conifer and broadleaved forests) and AWS (Automated Weather Station) of Pocheon urban area were used to quantify the air temperature reduction effect of forests, which is considered to be an eco-friendly solution for reducing the urban heat island intensity during summer. In June, July and August of 2016 and 2017, the average maximum air temperature differences between above and below canopy of forests, and between the forests and urban areas were $-1.9^{\circ}C$ and $-3.4^{\circ}C$ respectively, and they occurred at 17:00. However, there was no difference between conifer and broadleaved forests. The effect of air temperature reduction by the forests was positively correlated with accumulated evapotranspiration and solar radiation from 14:00 to 17:00 and showed a negative correlation with wind speed. We have developed a model to quantify the effect of air temperature reduction by forests using these variables. The nighttime air temperature reduction effect by forests was due to the generation of cold air from radiative cooling and the air temperature inversion phenomenon that occurs when the generated cold air moves down the side of mountain. The model was evaluated in Seoul by using 28 AWSs. The evaluation shows that the air temperature of each district in Seoul was negatively correlated with the area and size of the surrounding tall vegetation that drives vegetation evapotranspiration during the day. During the night, however, the size of the surrounding tall vegetation and the elevations of nearby mountains were the main influencing factors on the air temperature. Our research emphasizes the importance of the establishment and management of urban forests and the composition of wind roads from mountains for urban air temperature reduction.