• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• Korean Journal of Remote Sensing
• /
• v.34 no.2_1
• /
• pp.191-201
• /
• 2018

In this present study, we, for the first time, retrieved total column of ozone ($O_3$) and tropospheric ozone vertical profile using the Optimal Estimation (OE) method based on the MAX-DOAS measurement at the Yonsei University in Seoul, Korea. The optical density fitting is carried out using the OE method to calculate ozone columns. The optical density between the MAX-DOAS data obtained by dividing the measured intensities for each viewing elevated angle by those at the zenith angle. The retrieved total columns of the ozone are 375.4 and 412.6 DU in the morning (08:13) and afternoon (17:55) on 23 May, 2017, respectively. In addition, under 10 km altitude, the $O_3$ vertical profile was retrieved with about 5% of retrieval uncertainty. However, above 10 km altitude, the $O_3$ vertical profile retrieval uncertainty was increased (>10%). The spectral fitting errors are 16.8% and 19.1% in the morning and afternoon, respectively. The method suggested in this present study can be useful to measure the total ozone column using the ground-based hyper-spectral UV sensors.

• Geophysics and Geophysical Exploration
• /
• v.8 no.1
• /
• pp.7-17
• /
• 2005

The determination of seismic velocities in refractors for near-surface seismic refraction investigations is an ill-posed problem. Small variations in the computed time parameters can result in quite large lateral variations in the derived velocities, which are often artefacts of the inversion algorithms. Such artefacts are usually not recognized or corrected with forward modelling. Therefore, if detailed refractor models are sought with model based inversion, then detailed starting models are required. The usual source of artefacts in seismic velocities is irregular refractors. Under most circumstances, the variable migration of the generalized reciprocal method (GRM) is able to accommodate irregular interfaces and generate detailed starting models of the refractor. However, where the very-near-surface environment of the Earth is also irregular, the efficacy of the GRM is reduced, and weathering corrections can be necessary. Standard methods for correcting for surface irregularities are usually not practical where the very-near-surface irregularities are of limited lateral extent. In such circumstances, the GRM smoothing statics method (SSM) is a simple and robust approach, which can facilitate more-accurate estimates of refractor velocities. The GRM SSM generates a smoothing 'statics' correction by subtracting an average of the time-depths computed with a range of XY values from the time-depths computed with a zero XY value (where the XY value is the separation between the receivers used to compute the time-depth). The time-depths to the deeper target refractors do not vary greatly with varying XY values, and therefore an average is much the same as the optimum value. However, the time-depths for the very-near-surface irregularities migrate laterally with increasing XY values and they are substantially reduced with the averaging process. As a result, the time-depth profile averaged over a range of XY values is effectively corrected for the near-surface irregularities. In addition, the time-depths computed with a Bero XY value are the sum of both the near-surface effects and the time-depths to the target refractor. Therefore, their subtraction generates an approximate 'statics' correction, which in turn, is subtracted from the traveltimes The GRM SSM is essentially a smoothing procedure, rather than a deterministic weathering correction approach, and it is most effective with near-surface irregularities of quite limited lateral extent. Model and case studies demonstrate that the GRM SSM substantially improves the reliability in determining detailed seismic velocities in irregular refractors.

• Journal of the Korean earth science society
• /
• v.28 no.7
• /
• pp.936-946
• /
• 2007

Analytical geopotential field in balance with the sectoral mode (the first symmetric mode with respect to the equator) of the Rossby-Haurwitz wave on the inclined rotation axis was derived in presence of superrotation background flow. The balanced field was obtained by inverting the divergence equation with the time derivative being zero. The inversion consists of two steps, i.e., the evaluation of nonlinear forcing terms and the finding of analytical solutions based on the Poisson's equation. In the second step, the forcing terms in the from of Legendre function were readily inverted due to the fact that Legendre function is the eigenfunction of the spherical Laplacian operator, while other terms were solved either by introducing a trial function or by integrating the Legendre equation. The balanced field was found to be expressed with six zonal wavenumber components, and shown to be of asymmetric structure about the equator. In association with asymmetricity, the advantageous point of the balanced field as a validation method for the numerical model was addressed. In special cases where the strength of the background flow is a half of or exactly the same as the rotation rate of the Earth it was revealed that one of the zonal wavenumber components vanishes. The analytical balanced field was compared with the geopotential field which was obtained using a spherical harmonics spectral model. It was found that the normalized difference lied in the order of machine rounding, indicating the reliability of the analytical results. The stability of the sectoral mode of Rossby-Haurwitz wave and the associated balanced field was discussed, comparing with the flrst antisymmetric mode.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.6 no.4
• /
• pp.327-343
• /
• 2004

During tunnel construction, ground failures often occur due to existence of weak zones, such as faults, joints, and cavities, ahead of tunnel face. It is hard to detect effectively weak zones, which can lead underground structure to fail after excavation and before supporting, by using conventional characterization methods. In this study, an enhanced analytical method of predicting weak zones ahead of tunnel face is developed to overcome some problems in the conventional geophysical exploration methods. The analytical method is based on Coulomb's and Gauss' laws with considering the characteristics of electric fields subjected to rock mass. Using the developed method, closed form solutions are obtained to detect a spherical shaped zone and an oriented fault ahead of tunnel face respectively. The analytical results suggest that the presence of weak zones and their sizes, location, and states can be accurately predicted by combining a proper inversion process with resistance measured from several electrodes on the tunnel face. It appears that the skin depth or resistivity in rock mass is affected by the diameter of tunnel face, natural electric potential and noises induced by experimental measurement and spatial distribution of uncertain properties. The developed analytical solution is verified through experimental tests. About 1800 concrete blocks of 5cm by 5cm by 5cm in size are prepared and used to model a joint rock mass around tunnel face. Weak zones are simulated ahead of tunnel face with a material which has relatively higher conductivity than concrete blocks. Experimental results on the model test show a good agreement with analytical results.

• Korean Journal of Remote Sensing
• /
• v.33 no.6_1
• /
• pp.961-970
• /
• 2017

Despite of large demand for high spatial resolution products of aerosol properties from satellite remote sensing, it has been very difficult due to the weak signal by a single pixel and higher noise from clouds. In this study, aerosol retrieval algorithm with the high spatial resolution ($500m{\times}500m$) was developed using Geostationary Ocean Color Imager (GOCI) data during the Korea-US Air Quality (KORUS-AQ) period in May-June, 2016.Currently, conventional GOCI Yonsei aerosol retrieval(YAER) algorithm provides $6km{\times}6km$ spatial resolution product. The algorithm was tested for its best possible resolution of 500 m product based on GOCI YAER version 2 algorithm. With the new additional cloud masking, aerosol optical depth (AOD) is retrieved using the inversion method, aerosol model, and lookup table as in the GOCI YAER algorithm. In some cases, 500 m AOD shows consistent horizontal distribution and magnitude of AOD compared to the 6 km AOD. However, the 500 m AOD has more retrieved pixels than 6 km AOD because of its higher spatial resolution. As a result, the 500 m AOD exists around small clouds and shows finer features of AOD. To validate the accuracy of 500 m AOD, we used dataset from ground-based Aerosol Robotic Network (AERONET) sunphotometer over Korea. Even with the spatial resolution of 500 m, 500 m AOD shows the correlation coefficient of 0.76 against AERONET, and the ratio within Expected Error (EE) of 51.1%, which are comparable to the results of 6 km AOD.

• Journal of the Korean earth science society
• /
• v.29 no.1
• /
• pp.1-12
• /
• 2008

A geophysical mapping was performed for Hwasan caldera which is located in Euisung Sub-basin of the southeastern part of the Korean Peninsula. In order to overcome the limitation of the previous studies, remote sensing technic was used and dense potential data were obtained and analyzed. First, we analyzed geological lineament for target area using geological map, digital elevation model (DEM) data and satellite imagery. The results were greatly consistent with the previous studies, and showed that N-S and NW-SE direction are the most dominant one in target area. Second, based on the lineament analysis, highly dense gravity data were acquired in Euisung Sub-basin and an integrated interpretation considering air-born magnetic data was made to investigate the regional structure of the target area. The results of power spectrum analysis for the acquired potential data revealed that the subsurface of Euisung Sub-basin have two density discontinuities at about 1 km and 3-5 km depth. A 1 km depth discontinuity is thought as the depth of pyroclastic sedimentary rocks or igneous rocks which were intruded at the ring vent of Hwasan caldera, while a 3-5 km depth discontinuity seems to be associated with the depth of the basin basement. In addition, three-dimensional gravity inversion for the total area of Euisung Sub-basin was carried out, and the inversion results indicated two followings; 1) Cretaceous Palgongsan granite and Bulguksa intrusion rocks, which are located in southeastern part and northeastern part of Euisung Sub-basin, show two major low density anomalies, 2) pyroclastic rocks around Hwasan caldera also have lower density when compared with those of neighborhood regions and are extended to 1.5 km depth. However, a poor vertical resolution of potential survey makes it difficult to accurately delineate the detailed structure caldera which has a vertically developed characteristic in general. To overcome this limitation, integrated analysis was carried out using the magnetotelluric data on the corresponding area with potential data and we could obtain more reasonable geologic structure.

• Korean Journal of Remote Sensing
• /
• v.39 no.5_3
• /
• pp.933-948
• /
• 2023

Atmospheric aerosols not only have adverse effects on human health but also exert direct and indirect impacts on the climate system. Consequently, it is imperative to comprehend the characteristics and spatiotemporal distribution of aerosols. Numerous research endeavors have been undertaken to monitor aerosols, predominantly through the retrieval of aerosol optical depth (AOD) via satellite-based observations. Nonetheless, this approach primarily relies on a look-up table-based inversion algorithm, characterized by computationally intensive operations and associated uncertainties. In this study, a novel high-resolution AOD direct retrieval algorithm, leveraging machine learning, was developed using top-of-atmosphere reflectance data derived from the Geostationary Ocean Color Imager-II (GOCI-II), in conjunction with their differences from the past 30-day minimum reflectance, and meteorological variables from numerical models. The Light Gradient Boosting Machine (LGBM) technique was harnessed, and the resultant estimates underwent rigorous validation encompassing random, temporal, and spatial N-fold cross-validation (CV) using ground-based observation data from Aerosol Robotic Network (AERONET) AOD. The three CV results consistently demonstrated robust performance, yielding R²=0.70-0.80, RMSE=0.08-0.09, and within the expected error (EE) of 75.2-85.1%. The Shapley Additive exPlanations(SHAP) analysis confirmed the substantial influence of reflectance-related variables on AOD estimation. A comprehensive examination of the spatiotemporal distribution of AOD in Seoul and Ulsan revealed that the developed LGBM model yielded results that are in close concordance with AERONET AOD over time, thereby confirming its suitability for AOD retrieval at high spatiotemporal resolution (i.e., hourly, 250 m). Furthermore, upon comparing data coverage, it was ascertained that the LGBM model enhanced data retrieval frequency by approximately 8.8% in comparison to the GOCI-II L2 AOD products, ameliorating issues associated with excessive masking over very illuminated surfaces that are often encountered in physics-based AOD retrieval processes.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.05a
• /
• pp.79-86
• /
• 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 (190km

• Geophysics and Geophysical Exploration
• /
• v.17 no.3
• /
• pp.129-136
• /
• 2014

Electromagnetic (EM) methods are generally divided into frequency-domain EM (FDEM) and time-domain EM (TDEM) methods, depending on the source waveform. The FDEM and TDEM fields are mathematically related by the Fourier transformation, and the TDEM field can thus be obtained as the Fourier transformation of FDEM data. For modeling in time-domain, we can use fast frequency-domain modeling codes and then convert the results to the time domain with a suitable numerical method. Thus, frequency-to-time transformations are of interest to EM methods, which is generally attained through fast Fourier transform. However, faster frequency-to-time transformation is required for the 3D inversion of TDEM data or for the processing of vast air-borne TDEM data. The diffusion expansion method (DEM) is one of smart frequency-to-time transformation methods. In DEM, the EM field is expanded into a sequence of diffusion functions with a known frequency dependence, but with unknown diffusion-times that must be chosen based on the data to be transformed. Especially, accuracy of DEM is sensitive to the diffusion-time. In this study, we developed a method to determine the optimum range of diffusion-time values, minimizing the RMS error of the frequency-domain data approximated by the diffusion expansion. We confirmed that this method produces accurate results over a wider time range for a homogeneous half-space and two-layered model.