• Journal of Korean Society of Forest Science
• /
• v.102 no.4
• /
• pp.467-476
• /
• 2013

Site index information given from forest site map only exist in the sampled locations. In this study, site index for unsampled locations were estimated using kriging interpolation method which can interpolate values between point samples to generate a continuous surface. Site index of Pinus densiplora in Danyang area were calculated using Chapman-Richards model by plot unit. Then site index for unsampled locations were interpolated by theoretical variogram models and ordinary kriging. Also in order to assess parameter selection, cross-validation was performed by calculating mean error (ME), average standard error (ASE) and root mean square error (RMSE). In result, gaussian model was excluded because of the biggest relative nugget (37.40%). Then spherical model (16.80%) and exponential model (8.77%) were selected. Site index estimates of Pinus densiplora throughout the entire area in Danyang showed 4.39~19.53 based on exponential model, and 4.54~19.23 based on spherical model. By cross-validation, RMSE had almost no difference. But ME and ASE from spherical model were slightly lower than exponential model. Therefore site index prediction map from spherical model were finally selected. Average site index from site prediction map was 10.78. It can be expected that regional variance can be considered by site index prediction map in order to estimate forest biomass which has big spatial variance and eventually it is helpful to improve an accuracy of forest carbon estimation.

• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.1
• /
• pp.23-32
• /
• 1999

A new numerical method using solution-adaptive grids (SAG) is developed to solve the Richards' equation (RE) for unsaturated flow in porous media. Using a grid generation technique, the SAG method automatically redistributes a fixed number of grid points during the flow process, so that more grid points are clustered in regions of large solution gradients. The method uses the coordinate transformation technique to employ a new transformed RE, which is solved with the standard finite difference method. The movement of grid points is incorporated into the transformed RE, and therefore all computation is performed on fixed grid points of the transformed domain without using any interpolation techniques. Thus, numerical difficulties arising from the movement of the wetting front during the infiltration process have been substantially overcome by the new method. Numerical experiments for an one-dimensional infiltration problem are presented to compare the SAG method to the modified Picard method using a fixed grid. Results show that accuracy of a SAG solution using 41 nodes is comparable with the solution of the fixed grid method using 201 nodes, while it requires only 50% of the CPU time. The global mass balance and the convergence of SAG solutions are strongly affected by the time step size (Δt) and the weighting parameter (${\gamma}$) used for generating solution-adaptive grids. Thus, the method requires automated readjustment of Δt and ${\gamma}$ to yield mass-conservative and convergent solutions, although it may increase computational costs. The method can be effective especially for simulating unsaturated flow and other transport problems involving the propagation of a sharp-front.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.2 s.72
• /
• pp.171-185
• /
• 2006

A sophisticated reliability analysis method is proposed to evaluate the reliability of real nonlinear complicated dynamic structural systems excited by short duration dynamic loadings like earthquake motions by intelligently integrating the response surface method, the finite element method, the first-order reliability method, and the iterative linear interpolation scheme. The method explicitly considers all major sources of nonlinearity and uncertainty in the load and resistance-related random variables. The unique feature of the technique is that the seismic loading is applied in the time domain, providing an alternative to the classical random vibration approach. The four-parameter Richard model is used to represent the flexibility of connections of real steel frames. Uncertainties in the Richard parameters are also incorporated in the algorithm. The laterally flexible steel frame is then reinforced with reinforced concrete shear walls. The stiffness degradation of shear walls after cracking is also considered. The applicability of the method to estimate the reliability of real structures is demonstrated by considering three examples; a laterally flexible steel frame with fully restrained connections, the same steel frame with partially restrained connections with different rigidities, and a steel frame reinforced with concrete shear walls.

• Journal of the Korean Association of Geographic Information Studies
• /
• v.15 no.3
• /
• pp.36-51
• /
• 2012

The demand for a climatological dataset with a regular spaced grid is increasing in diverse fields such as ecological and hydrological modeling as well as regional climate impact studies. PRISM(Precipitation-Elevation Regressions on Independent Slopes Model) is a useful method to estimate high-altitude precipitation. However, it is not well discussed over the optimization of PRISM parameters and DEM(Digital Elevation Model) resolution in South Korea. This study developed the PRISM and then optimized parameters of the model and DEM resolution for producing a gridded annual average precipitation data of South Korea with 1km spatial resolution during the period 2000-2005. SCE-UA (Shuffled Complex Evolution-University of Arizona) method employed for the optimization. In addition, sensitivity analysis investigates the change in the model output with respect to the parameter and the DEM spatial resolution variations. The study result shows that maximum radius within which station search will be conducted is 67km. Minimum radius within which all stations are included is 31km. Minimum number of stations required for cell precipitation and elevation regression calculation is four. Optimizing DEM resolution is $1{\times}1km$. This study also shows that the PRISM output very sensitive to DEM spatial resolution variations. This study contributes to improving the accuracy of PRISM technique as it applies to South Korea.

• Journal of the Korean Geotechnical Society
• /
• v.34 no.5
• /
• pp.19-35
• /
• 2018

The necessity of predicting the spatial information of the site-specific seismic response, which is essential information for the comprehensive earthquake disaster countermeasures, is increasing for the mid-west urban areas where the earthquake-induced damages can be increased due to frequent occurrence of mid-scale earthquake such as 2016 Gyeongju Earthquake and 2017 Pohang Earthquake. Especially, researches on strategic securing of site survey datasets and understanding the site-specific site response characteristics were conducted for Gyeonggi-do, South Korea. In this study, a GIS-based framework for site-specific assessment of site response and induced vulnerable area in Gyeonggi-do, South Korea was proposed. Geo-Data based on GIS platform was constructed for regional estimation of geotechnical characteristics by collecting borehole and land coverage datasets. And the geo-spatial grid information was developed for deriving spatial distribution of geotechnical layer and site response parameters based on the optimization of the geostatistical interpolation method. Accordingly, base information for Improving earthquake preparedness measures was derived as seismic zonation map with administrative sub-units considering the quantitative site effect of Gyeonggi-do.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.21 no.3
• /
• pp.167-174
• /
• 2019

A process-oriented crop growth model can simulate the biophysical process of rice under diverse environmental and management conditions, which would make it more versatile than an empirical crop model. In the present study, we examined chronology and background of the development of the rice growth models in Korea, which would provide insights on the needs for improvement of the models. The rice crop growth models were introduced in Korea in the late 80s. Until 2000s, these crop models have been used to simulate the yield in a specific area in Korea. Since then, improvement of crop growth models has been made to take into account biological characteristics of rice growth and development in more detail. Still, the use of the crop growth models has been limited to the assessment of climate change impact on crop production. Efforts have been made to apply the crop growth model, e.g., the CERES-Rice model, to develop decision support system for crop management at a farm level. However, the decision support system based on a crop growth model was attractive to a small number of stakeholders most likely due to scarcity of on-site weather data and reliable parameter sets for cultivars grown in Korea. The wide use of the crop growth models would be facilitated by approaches to extend spatial availability of reliable weather data, which could be either measured on-site or estimates using spatial interpolation. New approaches for calibration of cultivar parameters for new cultivars would also help lower hurdles to crop growth models.

• Proceedings of the Korean Association of Geographic Inforamtion Studies Conference
• /
• 2003.04a
• /
• pp.363-368
• /
• 2003

We have successfully developed a more effective algorithm to extract the lineament in the area covered by wide alluvial deposits characterized by a relatively narrow range of brightness in the Landsat TM image, while the currently used algorithm is limited to the mountainous areas. In the new algorithm, flat areas mainly consisting of alluvial deposits were selected using the Local Enhancement from the Digital Elevation Model (DEM). The aspect values were obtained by 3${\times}$3 moving windowing of Zevenbergen & Thorno's Method, and then the slopes of the study area were determined using the aspect values. After the lineament factors in the alluvial deposits were revealed by comparing the threshold values, the first rank lineament under the alluvial deposits were extracted using the Hough transform In order to extract the final lineament, the lowest points under the alluvial deposits in a given topographic section perpendicular to the first rank lineament were determined through the spline interpolation, and then the final lineament were chosen through Hough transform using the lowest points. The algorithm developed in this study enables us to observe a clearer lineament in the areas covered by much larger alluvial deposits compared with the results extracted using the conventional existing algorithm. There exists, however, some differences between the first rank lineament, obtained using the aspect and the slope, and the final lineament. This study shows that the new algorithm more effectively extracts the lineament in the area covered with wide alluvlal deposits than in the areas of converging slope, areas with narrow alluvial deposits or valleys.

• Radiation Oncology Journal
• /
• v.15 no.1
• /
• pp.71-78
• /
• 1997

Purpose : To collect beam data for dynamic wedge fields using conventional measurement tools without the multi-detector system, such as the linear diode detectors or ionization chambers. Materials and Methods : The accelerator CL 2100 C/D has two photon energies of 6MV and 15MV with dynamic wedge an91es of 15o, 30o, 45o and 60o. Wedge transmission factors, percentage depth doses(PDD's) and dose Profiles were measured. The measurements for wedge transmission factors are performed for field sizes ranging from $4\times4cm^2\;to\;20\times20cm^2$ in 1-2cm steps. Various rectangular field sizes are also measured for each photon energy of 6MV and 15MV, with the combination of each dynamic wedge angle of 15o 30o. 45o and 60o. These factors are compared to the calculated wedge factors using STT(Segmented Treatment Table) value. PDD's are measured with the film and the chamber in water Phantom for fixed square field. Converting parameters for film data to chamber data could be obtained from this procedure. The PDD's for dynamic wedged fields could be obtained from film dosimetry by using the converting parameters without using ionization chamber. Dose profiles are obtained from interpolation and STT weighted superposition of data through selected asymmetric static field measurement using ionization chamber. Results : The measured values of wedge transmission factors show good agreement to the calculated values The wedge factors of rectangular fields for constant V-field were equal to those of square fields The differences between open fields' PDDs and those from dynamic fields are insignificant. Dose profiles from superposition method showed acceptable range of accuracy(maximum 2% error) when we compare to those from film dosimetry. Conclusion : The results from this superposition method showed that commissionning of dynamic wedge could be done with conventional dosimetric tools such as Point detector system and film dosimetry winthin maximum 2% error range of accuracy.

• Korean Journal of Remote Sensing
• /
• v.37 no.6_1
• /
• pp.1611-1623
• /
• 2021

Thisstudy demonstratessoil moisture and evapotranspiration performance using Korea Land Data Assimilation System (KLDAS) under Korea Land Information System (KLIS). Spin-up was repeated 8 times in 2018. In addition, low-resolution and high-resolution meteorological data were generated using meteorological data observed by Korea Meteorological Administration (KMA), Rural Development Administration (RDA), Korea Rural Community Corporation (KRC), Korea Hydro & Nuclear Power Co.,Ltd. (KHNP), Korea Water Resources Corporation (K-water), and Ministry of Environment (ME), and applied to KLDAS. And, to confirm the degree of accuracy improvement of Korea Low spatial resolution (hereafter, K-Low; 0.125°) and Korea High spatial resolution (hereafter, K-High; 0.01°), soil moisture and evapotranspiration to which Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) and ASOS-Spatial (ASOS-S) used in the previous study were applied were evaluated together. As a result, optimization of the initial boundary condition requires 2 time (58 point), 3 time (6 point), and 6 time (3 point) spin-up for soil moisture. In the case of evapotranspiration, 1 time (58 point) and 2 time (58 point) spin-ups are required. In the case of soil moisture to which MERRA-2, ASOS-S, K-Low, and K-High were applied, the mean of R² were 0.615, 0.601, 0.594, and 0.664, respectively, and in the case of evapotranspiration, the mean of R² were 0.531, 0.495, 0.656, and 0.677, respectively, indicating the accuracy of K-High was rated as the highest. The accuracy of KLDAS can be improved by securing a large number of ground observation data through the results of this study and generating high-resolution grid-type meteorological data. However, if the meteorological condition at each point is not sufficiently taken into account when converting the point data into a grid, the accuracy is rather lowered. For a further study, it is expected that higher quality data can be produced by generating and applying grid-type meteorological data using the parameter setting of IDW or other interpolation techniques.

• The Journal of Engineering Geology
• /
• v.33 no.1
• /
• pp.105-119
• /
• 2023

Water supply is decreasing due to climate change, and coastal and island regions are highly dependent on groundwater, reducing the amount of available water. For sustainable water supply in coastal and island regions, it is necessary to accurately diagnose the current condition and efficiently distribute and manage water. For a precise analysis of the groundwater flow in the coastal island region, submarine fresh groundwater discharge was calculated for the Seongsan basin in the eastern part of Jeju Island. Two methods were used to estimate the thickness of the fresh groundwater. One method employed vertical interpolation of measured electrical conductivity in a multi depth monitoring well; the other used theoretical Ghyben-Herzberg ratio. The value using the Ghyben-Herzberg ratio makes it impossible to accurately estimate the changing salt-saltwater interface, and the value analyzed by electrical conductivity can represent the current state of the freshwater-saltwater interface. Observed parameter was distributed on a virtual grid. The average of submarine fresh groundwater discharge fluxes for the virtual grid was determined as the watershed's representative flux. The submarine fresh groundwater discharge and flux distribution by year were also calculated at the basin scale. The method using electrical conductivity estimated the submarine fresh groundwater discharge from 2018 to 2020 to be 6.27 × 10⁶ m³/year; the method using the Ghyben-Herzberg ratio estimated a discharge of 10.87 × 10⁶ m³/year. The results presented in this study can be used as basis data for policies that determine sustainable water supply by using precise water budget analysis in coastal and island areas.