• Journal of the Korean Geotechnical Society
• /
• v.36 no.2
• /
• pp.5-15
• /
• 2020

The liquefaction phenomenon that occurred during the Pohang earthquake (M_L=5.4) brought new awareness to the people about the risk of liquefaction caused by the earthquake. Liquefaction hazard maps with 2 km grid made in 2014 used more than 100,000 borehole data for the whole country, and regions without soil investigation data were produced using interpolation. In the mapping of macro liquefaction hazard for the whole country, the site amplification effect and the ground water level 0 m were considered. Recently, the Ministry of Public Administration and Security (2018) published a new site classification method and amplification coefficient of the common standard for seismic design. Therefore, it is necessary to rewrite the liquefaction hazard map reflecting the revised amplification coefficient. In this study, the results of site classification according to the average shear wave velocity in soils before and after revision were compared in the whole country. Also, liquefaction assessment results were compared in Gangseo-gu, Busan. At this time, two ground accelerations corresponding to the 500 and 1,000 years of return period and two ground water table, 5 m for the average condition and 0 m the extreme condition were applied. In the drawing of liquefaction hazard map, a 500 m grid was applied to secure a resolution higher than the previous 2 km grid. As a result, the ground conditions that were classified as S_C and S_D grounds based on the existing site classification standard were reclassified as S₂, S₃, and S₄ through the revised site classification standard. Also, the result of the Liquefaction assessments with a return period of 500 years and 1,000 years resulted in a relatively overestimation of the LPI applied with the ground amplification factor before revision. And the results of this study have a great influence on the liquefaction assessment, which is the basis of the creation of the regional liquefaction hazard map using the amplification factor.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.24 no.2
• /
• pp.63-77
• /
• 2022

Our aim was to reduce estimation errors of a wind velocity model used as an early warning system for weather risk management in the agricultural sector. The Rural Development Administration (RDA) agricultural weather observation network's wind velocity data and its corresponding estimated data from January to December 2020 were used to calculate linear regression equations (Y = aX + b). In each linear regression, the wind estimation error at 87 points and eight time slots per day (00:00, 03:00, 06:00, 09.00, 12.00, 15.00, 18.00, and 21:00) is the dependent variable (Y), while the estimated wind velocity is the independent variable (X). When the correlation coefficient exceeded 0.5, the regression equation was used as the wind velocity correction equation. In contrast, when the correlation coefficient was less than 0.5, the mean error (ME) at the corresponding points and time slots was substituted as the correction value instead of the regression equation. To enable the use of wind velocity model at a national scale, a distribution map with a grid resolution of 250 m was created. This objective was achieved b y performing a spatial interpolation with an inverse distance weighted (IDW) technique using the regression coefficients (a and b), the correlation coefficient (R), and the ME values for the 87 points and eight time slots. Interpolated grid values for 13 weather observation points in rural areas were then extracted. The wind velocity estimation errors for 13 points from January to December 2019 were corrected and compared with the system's values. After correction, the mean ME of the wind velocities reduced from 0.68 m/s to 0.45 m/s, while the mean RMSE reduced from 1.30 m/s to 1.05 m/s. In conclusion, the system's wind velocities were overestimated across all time slots; however, after the correction model was applied, the overestimation reduced in all time slots, except for 15:00. The ME and RMSE improved b y 33% and 19.2%, respectively. In our system, the warning for wind damage risk to crops is driven by the daily maximum wind speed derived from the daily mean wind speed obtained eight times per day. This approach is expected to reduce false alarms within the context of strong wind risk, by reducing the overestimation of wind velocities.

• Journal of the Korean Association of Geographic Information Studies
• /
• v.18 no.2
• /
• pp.174-189
• /
• 2015

The objective of this study was to reconstruct spatial information using the results of the investigation and evaluation of the health of the living organisms, habitat, and water quality at the investigation points for the aquatic ecosystem health of the Nakdong River basin, to support the rational decision making of the aquatic ecosystem preservation and restoration policies of the Nakdong River basin using spatial analysis techniques, and to present efficient management methods. To analyze the aquatic ecosystem health of the Nakdong River basin, punctiform data were constructed based on the position information of each point with the aquatic ecosystem health investigation and evaluation results of 250 investigation sections. To apply the spatial analysis technique, the data need to be reconstructed into areal data. For this purpose, spatial influence and trends were analyzed using the Kriging interpolation(ArcGIS 10.1, Geostatistical Analysis), and were reconstructed into areal data. To analyze the spatial distribution characteristics of the Nakdong River basin health based on these analytical results, hotspot(Getis-Ord Gi, $G^*_i$), LISA(Local Indicator of Spatial Association), and standard deviational ellipse analyses were used. The hotspot analysis results showed that the hotspot basins of the biotic indices(TDI, BMI, FAI) were the Andong Dam upstream, Wangpicheon, and the Imha Dam basin, and that the health grades of their biotic indices were good. The coldspot basins were Nakdong River Namhae, the Nakdong River mouth, and the Suyeong River basin. The LISA analysis results showed that the exceptional areas were Gahwacheon, the Hapcheon Dam, and the Yeong River upstream basin. These areas had high bio-health indices, but their surrounding basins were low and required management for aquatic ecosystem health. The hotspot basins of the physicochemical factor(BOD) were the Nakdong River downstream basin, Suyeong River, Hoeya River, and the Nakdong River Namhae basin, whereas the coldspot basins were the upstream basins of the Nakdong River tributaries, including Andong Dam, Imha Dam, and Yeong River. The hotspots of the habitat and riverside environment factor(HRI) were different from the hotspots and coldspots of each factor in the LISA analysis results. In general, the habitat and riverside environment of the Nakdong River mainstream and tributaries, including the Nakdong river upstream, Andong Dam, Imha Dam, and the Hapcheon Dam basin, had good health. The coldspot basins of the habitat and riverside environment also showed low health indices of the biotic indices and physicochemical factors, thus requiring management of the habitat and riverside environment. As a result of the time-series analysis with a standard deviation ellipsoid, the areas with good aquatic ecosystem health of the organisms, habitat, and riverside environment showed a tendency to move northward, and the BOD results showed different directions and concentrations by the year of investigation. These aquatic ecosystem health analysis results can provide not only the health management information for each investigation spot but also information for managing the aquatic ecosystem in the catchment unit for the working research staff as well as for the water environment researchers in the future, based on spatial information.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.26-32
• /
• 1999

Among the various seismic data processing sequences, the velocity analysis is the most time consuming and man-hour intensive processing steps. For the production seismic data processing, a good velocity analysis tool as well as the high performance computer is required. The tool must give fast and accurate velocity analysis. There are two different approches in the velocity analysis, batch and interactive. In the batch processing, a velocity plot is made at every analysis point. Generally, the plot consisted of a semblance contour, super gather, and a stack pannel. The interpreter chooses the velocity function by analyzing the velocity plot. The technique is highly dependent on the interpreters skill and requires human efforts. As the high speed graphic workstations are becoming more popular, various interactive velocity analysis programs are developed. Although, the programs enabled faster picking of the velocity nodes using mouse, the main improvement of these programs is simply the replacement of the paper plot by the graphic screen. The velocity spectrum is highly sensitive to the presence of the noise, especially the coherent noise often found in the shallow region of the marine seismic data. For the accurate velocity analysis, these noise must be removed before the spectrum is computed. Also, the velocity analysis must be carried out by carefully choosing the location of the analysis point and accuarate computation of the spectrum. The analyzed velocity function must be verified by the mute and stack, and the sequence must be repeated most time. Therefore an iterative, interactive, and unified velocity analysis tool is highly required. An interactive velocity analysis program, xva(X-Window based Velocity Analysis) was invented. The program handles all processes required in the velocity analysis such as composing the super gather, computing the velocity spectrum, NMO correction, mute, and stack. Most of the parameter changes give the final stack via a few mouse clicks thereby enabling the iterative and interactive processing. A simple trace indexing scheme is introduced and a program to nike the index of the Geobit seismic disk file was invented. The index is used to reference the original input, i.e., CDP sort, directly A transformation techinique of the mute function between the T-X domain and NMOC domain is introduced and adopted to the program. The result of the transform is simliar to the remove-NMO technique in suppressing the shallow noise such as direct wave and refracted wave. However, it has two improvements, i.e., no interpolation error and very high speed computing time. By the introduction of the technique, the mute times can be easily designed from the NMOC domain and applied to the super gather in the T-X domain, thereby producing more accurate velocity spectrum interactively. The xva program consists of 28 files, 12,029 lines, 34,990 words and 304,073 characters. The program references Geobit utility libraries and can be installed under Geobit preinstalled environment. The program runs on X-Window/Motif environment. The program menu is designed according to the Motif style guide. A brief usage of the program has been discussed. The program allows fast and accurate seismic velocity analysis, which is necessary computing the AVO (Amplitude Versus Offset) based DHI (Direct Hydrocarn Indicator), and making the high quality seismic sections.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.48 no.5
• /
• pp.25-33
• /
• 2011

This work proposes a skinny-type 10b 50MS/s 0.13um CMOS three-step pipeline ADC for CIS applications. Analog circuits for CIS applications commonly employ a high supply voltage to acquire a sufficiently acceptable dynamic range, while digital circuits use a low supply voltage to minimize power consumption. The proposed ADC converts analog signals in a wide-swing range to low voltage-based digital data using both of the two supply voltages. An op-amp sharing technique employed in residue amplifiers properly controls currents depending on the amplification mode of each pipeline stage, optimizes the performance of op-amps, and improves the power efficiency. In three FLASH ADCs, the number of input stages are reduced in half by the interpolation technique while each comparator consists of only a latch with low kick-back noise based on pull-down switches to separate the input nodes and output nodes. Reference circuits achieve a required settling time only with on-chip low-power drivers and digital correction logic has two kinds of level shifter depending on signal-voltage levels to be processed. The prototype ADC in a 0.13um CMOS to support 0.35um thick-gate-oxide transistors demonstrates the measured DNL and INL within 0.42LSB and 1.19LSB, respectively. The ADC shows a maximum SNDR of 55.4dB and a maximum SFDR of 68.7dB at 50MS/s, respectively. The ADC with an active die area of 0.53$mm^2$ consumes 15.6mW at 50MS/s with an analog voltage of 2.0V and two digital voltages of 2.8V ($=D_H$) and 1.2V ($=D_L$).

• Journal of the Korean Association of Geographic Information Studies
• /
• v.14 no.3
• /
• pp.78-95
• /
• 2011

This study has developed various DEMs with different spatial resolutions using many different interpolation methods with the aid of a 1:5,000 digital map. In addition, this study has evaluated the vertical accuracy of various DEMs constructed by check point data obtained from the network RTK GPS survey. The obtained results suggest that a DEM developed from the TIN-based Terrain method performs well in evaluating height restriction of buildings in a flying safety area considering general RMSE values, land-type RMSE values and profile evaluation results, etc. And, it has been found that three meters is the right spatial resolution for a DEM in evaluating height restriction of buildings in a flying safety area. Meanwhile, elevation values obtained by the DEM are not point estimation values but interval estimation values. This can be used to check whether the height of buildings in the vicinity of an airfield violates height limitation values of the area. To check whether the height of buildings measured in interval estimation values violates height limitation values of the area, this study has adopted three steps: 1) high probability of violation, 2) low probability of violation, 3) inconclusiveness about the violation. The obtained results will provide an important basis for developing a GIS related to the evaluation of height restriction of buildings in the vicinity of an airfield. Furthermore, although results are limited to the study area, the vertical accuracy values of the DEM constructed from a two-dimensional digital map may provide useful information to researchers who try to use DEMs.

• Journal of the Korean Association of Geographic Information Studies
• /
• v.19 no.3
• /
• pp.61-74
• /
• 2016

The purpose of this paper is to analyze the relationship between the location of the epicenter of a medium-sized earthquake(magnitude 4.8) that occurred on January 20, 2007 in the Odaesan area with lineament features using a shaded relief map(1/25,000 scale) and satellite images from LANDSAT-8 and KOMPSAT-2. Previous studies have analyzed lineament features in tectonic settings primarily by examining two-dimensional satellite images and shaded relief maps. These methods, however, limit the application of the visual interpretation of relief features long considered as the major component of lineament extraction. To overcome some existing limitations of two-dimensional images, this study examined three-dimensional images, produced from a Digital Elevation Model and drainage network map, for lineament extraction. This approach reduces mapping errors introduced by visual interpretation. In addition, spline interpolation was conducted to produce density maps of lineament frequency, intersection, and length required to estimate the density of lineament at the epicenter of the earthquake. An algorithm was developed to compute the Value of the Relative Density(VRD) representing the relative density of lineament from the map. The VRD is the lineament density of each map grid divided by the maximum density value from the map. As such, it is a quantified value that indicates the concentration level of the lineament density across the area impacted by the earthquake. Using this algorithm, the VRD calculated at the earthquake epicenter using the lineament's frequency, intersection, and length density maps ranged from approximately 0.60(min) to 0.90(max). However, because there were differences in mapped images such as those for solar altitude and azimuth, the mean of VRD was used rather than those categorized by the images. The results show that the average frequency of VRD was approximately 0.85, which was 21% higher than the intersection and length of VRD, demonstrating the close relationship that exists between lineament and the epicenter. Therefore, it is concluded that the density map analysis described in this study, based on lineament extraction, is valid and can be used as a primary data analysis tool for earthquake research in the future.

• Journal of Korean Society of Forest Science
• /
• v.106 no.3
• /
• pp.320-329
• /
• 2017

This study aims to evaluate whether the accuracy of LiDAR DEM is affected by the changes of the five input levels ('1','3','5','7' and '9') of median parameter ($F_{md}$), mean parameter ($F_{mn}$) of the Filtering Algorithm (FA) in the GroundFilter module and median parameter ($I_{md}$), mean parameter ($I_{mn}$) of the Interpolation Algorithm (IA) in the GridSurfaceCreate module of the FUSION in order to present the combination of parameter levels producing the most accurate LiDAR DEM. The accuracy is measured by the residuals calculated by difference between the field elevation values and their corresponding DEM elevation values. A multi-way ANOVA is used to statistically examine whether there are effects of parameter level changes on the means of the residuals. The Tukey HSD is conducted as a post-hoc test. The results of the multi- way ANOVA test show that the changes in the levels of $F_{md}$, $F_{mn}$, $I_{mn}$ have significant effects on the DEM accuracy with the significant interaction effect between $F_{md}$ and $F_{mn}$. Therefore, the level of $F_{md}$, $F_{mn}$, and the interaction between two variables are considered to be factors affecting the accuracy of LiDAR DEM as well as the level of $I_{mn}$. As the results of the Tukey HSD test on the combination levels of $F_{md}{\ast}F_{mn}$, the mean of residuals of the '$9{\ast}3$' combination provides the highest accuracy while the '$1{\ast}1$' combination provides the lowest one. Regarding $I_{mn}$ levels, the mean of residuals of the both '3' and '1' provides the highest accuracy. This study can contribute to improve the accuracy of the forest attributes as well as the topographic information extracted from the LiDAR data.

• The Journal of the Acoustical Society of Korea
• /
• v.13 no.6
• /
• pp.75-82
• /
• 1994

The ray paths and travel times of sound wave in the ocean depend on the physical properties of the propagating media. Ocean Acoustic Tomography(OAT), which is inversely estimate the travel time variations between fixed sources and receivers the physical properties of the corresponding media can he understood. To apply ocean survey technology by using the OAT, the tomographic procedure requires forward problem that variation of the travel times be identified with the variability of the medium. Also, received signals must be satisfied the necessary conditions of ray path stability, identification and resolution in order for OAT to work. The canonical ocean has been determined based on the historical data and its travel time and ray path are used as reference values. The sound speed of canonical ocean in the East Sea is about 1523 m/s at the surface and 1458 m/s at the sound channel axis(400m). Sound speeds in the East Sea are perturbed by warm eddy whose horizontal extension is more than 100 km with deeper than 200 m in depth scale. In this study, an acoustic source and receiver are placed at the depth above the sound channel axis, 350 m, and are separated by 200 km range. Ray paths are identified by the ray theory methed in a range dependent medium whose sound speeds are functions of a range and depth. The eigenray information obtained from interpolation between the rays bracketing the receiver are used to simulate the received signal by convolution of source signal with the eigenray informations. The source signal is taken as a 400 Hz rectangular pulse signal, bandwidth is 16 Hz and pulse length is 64 ms. According to the analysis of the received signal and identified ray path by using numerical model of underwater sound propagation, simulated signals satisfy the necessary conditions of OAT, applied in the East Sea.

• Economic and Environmental Geology
• /
• v.36 no.4
• /
• pp.305-312
• /
• 2003

This study was performed to delineate the subsurface geology, geologic structure, and distribution pattern of the Palgongsan granitic body, and to reveal the relationship between the Kyeongsang basin and Yongnam massif by gravity survey. The study area is located between the latitude of 35$^{\circ}$45'-36$^{\circ}$21'N and longitude of 128$^{\circ}$15'-129$^{\circ}$00'E. Total of 966 gravity data measured by Seoul National University, KlGAM(Korea Institute of Geology, Mining ＆ Materials), Pusan National University and Yonsei University were used. The Bouguer gravity anomaly in the study area ranges from -12.88 to 26.01 mgal with a mean value of 11.27 mgal. A very low anomaly zone is located in the Yongnam massif in west of the study area. The anomaly value increases going from west to east. A low anomaly distribution in Palgongsan granite and Yongnam massif is interpreted as the effect of their lower density than that of Kyeongsang Super Group. Power spectrum analysis is applied to evaluate the average depth of basement the Kyeongsang Basin and Conrad discontinuity from gravity anomaly. The average depths of density discontinuities are calculated 10.45 km and 4.9 km, and these are interpreted as Conrad discontinuity and depth of basement of the Kyeongsang Basin, respectively. The depth of Palgongsan granite is derived by means of 2-dimensional modeling and it decreases gradually toward the east. The gravity anomaly east of the study area decreases abruptly due to Shingryeong fault and Nogosan ring fault. Two deepest and sharp roots of Palgongsan granite are recognized by 2-dimensional modeling of each profiles. The depths of those roots are 5.3 km on a profile AA' and 7 km on a profile BB' which is the maximum depth of Palgongsan granite. Small granitic bodies are also seen to be intruded around the Palgongsan granite. The root of Palgongsan granite is shown by 3-dimensional analysis based on the interpolation of 2-dimensional modeling along each profiles to exist in the southwest vicinity of Palgongsan granite. The total volume of Palgongsan granite is approximately 31.211 $Km^3$.