• Korean Journal of Remote Sensing
• /
• v.37 no.5_1
• /
• pp.1135-1147
• /
• 2021

In order to geometrically correct high-resolution satellite imagery, the sensor modeling process that restores the geometric relationship between the satellite sensor and the ground surface at the image acquisition time is required. In general, high-resolution satellites provide RPC (Rational Polynomial Coefficient) information, but the vendor-provided RPC includes geometric distortion caused by the position and orientation of the satellite sensor. GCP (Ground Control Point) is generally used to correct the RPC errors. The representative method of acquiring GCP is field survey to obtain accurate ground coordinates. However, it is difficult to find the GCP in the satellite image due to the quality of the image, land cover change, relief displacement, etc. By using image maps acquired from various sensors as reference data, it is possible to automate the collection of GCP through the image matching algorithm. In this study, the RPC of KOMPSAT-3A satellite image was corrected through the extracted matching point using the UAV (Unmanned Aerial Vehichle) imagery. We propose a pre-porocessing method for the extraction of matching points between the UAV imagery and KOMPSAT-3A satellite image. To this end, the characteristics of matching points extracted by independently applying the SURF (Speeded-Up Robust Features) and the phase correlation, which are representative feature-based matching method and area-based matching method, respectively, were compared. The RPC adjustment parameters were calculated using the matching points extracted through each algorithm. In order to verify the performance and usability of the proposed method, it was compared with the GCP-based RPC correction result. The GCP-based method showed an improvement of correction accuracy by 2.14 pixels for the sample and 5.43 pixelsfor the line compared to the vendor-provided RPC. In the proposed method using SURF and phase correlation methods, the accuracy of sample was improved by 0.83 pixels and 1.49 pixels, and that of line wasimproved by 4.81 pixels and 5.19 pixels, respectively, compared to the vendor-provided RPC. Through the experimental results, the proposed method using the UAV imagery presented the possibility as an alternative to the GCP-based method for the RPC correction.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.23 no.4
• /
• pp.198-221
• /
• 2021

After large-scale reforestation in the 1960s and 1970s, forests in Korea have gradually been aging. Net ecosystem CO₂ exchange of old-growth forests is theoretically near zero; however, it can be a CO₂ sink or source depending on the intervention of disturbance or management. In this study, we report the CO₂ budget dynamics of the Gwangneung deciduous old-growth forest (GDK) in Korea and examined the following two questions: (1) is the preserved GDK indeed CO₂ neutral as theoretically known? and (2) can we explain the dynamics of CO₂ budget by the common mechanisms reported in the literature? To answer, we analyzed the 15-year long CO₂ flux data measured by eddy covariance technique along with other biometeorological data at the KoFlux GDK site from 2006 to 2020. The results showed that (1) GDK switched back-and-forth between sink and source of CO₂ but averaged to be a week CO₂ source (and turning to a moderate CO₂ source for the recent five years) and (2) the interannual variability of solar radiation, growing season length, and leaf area index showed a positive correlation with that of gross primary production (GPP) (R²=0.32~0.45); whereas the interannual variability of both air and surface temperature was not significantly correlated with that of ecosystem respiration (RE). Furthermore, the machine learning-based model trained using the dataset of early monitoring period (first 10 years) failed to reproduce the observed interannual variations of GPP and RE for the recent five years. Biomass data analysis suggests that carbon emissions from coarse woody debris may have contributed partly to the conversion to a moderate CO₂ source. To properly understand and interpret the long-term CO₂ budget dynamics of GDK, new framework of analysis and modeling based on complex systems science is needed. Also, it is important to maintain the flux monitoring and data quality along with the monitoring of coarse woody debris and disturbances.

• Tunnel and Underground Space
• /
• v.32 no.1
• /
• pp.30-58
• /
• 2022

The deep geological repository for high-level radioactive waste disposal is a multi barrier system comprised of engineered barriers and a natural barrier. The long-term integrity of the deep geological repository is affected by the coupled interactions between the individual barrier components. Erosion and piping phenomena in the compacted bentonite buffer due to buffer-rock interactions results in the removal of bentonite particles via groundwater flow and can negatively impact the integrity and performance of the buffer. Rapid groundwater inflow at the early stages of disposal can lead to piping in the bentonite buffer due to the buildup of pore water pressure. The physiochemical processes between the bentonite buffer and groundwater lead to bentonite swelling and gelation, resulting in bentonite erosion from the buffer surface. Hence, the evaluation of erosion and piping occurrence and its effects on the integrity of the bentonite buffer is crucial in determining the long-term integrity of the deep geological repository. Previous studies on bentonite erosion and piping failed to consider the complex coupled thermo-hydro-mechanical-chemical behavior of bentonite-groundwater interactions and lacked a comprehensive model that can consider the complex phenomena observed from the experimental tests. In this technical note, previous studies on the mechanisms, lab-scale experiments and numerical modeling of bentonite buffer erosion and piping are introduced, and the future expected challenges in the investigation of bentonite buffer erosion and piping are summarized.

• Geophysics and Geophysical Exploration
• /
• v.26 no.1
• /
• pp.18-30
• /
• 2023

Submarine mud volcanos are topographic features that resemble volcanoes, and are formed due to eruptions of fluidized or gasified sediment material. They have gained attention as a source of subsurface heat, sediment, or hydrocarbons supplied to the surface. In the continental slope of the Canadian Beaufort Sea, mud volcano exists at various water depths. The MV420, is an active mud volcano erupting at a water depth of 420 meters, and it has been the subject of extensive study. The Korea Polar Research Institute(KOPRI) collected high-resolution seismic data and heat flow data around the caldera of the mud volcano. By analyzing the multi-channel seismic data, we confirmed the reverse-polarity reflector assumed by a gas hydrate-related bottom simulating reflector(BSR). To further elucidate the relationship between the BSR and gas hydrates, as well as the thermal structure of the mud volcano, a numerical geothermal model was developed based on the steady-state heat equation. Using this model, we estimated the base of the gas hydrate stability zone and found that the BSR depth estimated by multi-channel seismic data and the bottom of the gas hydrate stability zone were in good agreement., This suggests the presence of gas hydrates, and it was determined that the depth of the gas hydrate was likely up to 50 m, depending on the distance from the mud conduit. Thus, this depth estimate slightly differs from previous studies.

• Journal of the Korean earth science society
• /
• v.43 no.6
• /
• pp.691-711
• /
• 2022

The state-of-the-art Earth system model as a virtual Earth is required for studies of current and future climate change or climate crises. This complex numerical model can account for almost all human activities and natural phenomena affecting the atmosphere of Earth. The Unified Model (UM) from the United Kingdom Meteorological Office (UK Met Office) is among the best Earth system models as a scientific tool for studying the atmosphere. However, owing to the expansive numerical integration cost and substantial output size required to maintain the UM, individual research groups have had to rely only on supercomputers. The limitations of computer resources, especially the computer environment being blocked from outside network connections, reduce the efficiency and effectiveness of conducting research using the model, as well as improving the component codes. Therefore, this study has presented detailed guidance for installing a new version of the UM on high-performance parallel computers (Linux clusters) owned by individual researchers, which would help researchers to easily work with the UM. The numerical integration performance of the UM on Linux clusters was also evaluated for two different model resolutions, namely N96L85 (1.875° ×1.25° with 85 vertical levels up to 85 km) and N48L70 (3.75° ×2.5° with 70 vertical levels up to 80 km). The one-month integration times using 256 cores for the AMIP and CMIP simulations of N96L85 resolution were 169 and 205 min, respectively. The one-month integration time for an N48L70 AMIP run using 252 cores was 33 min. Simulated results on 2-m surface temperature and precipitation intensity were compared with ERA5 re-analysis data. The spatial distributions of the simulated results were qualitatively compared to those of ERA5 in terms of spatial distribution, despite the quantitative differences caused by different resolutions and atmosphere-ocean coupling. In conclusion, this study has confirmed that UM can be successfully installed and used in high-performance Linux clusters.

• Economic and Environmental Geology
• /
• v.56 no.1
• /
• pp.103-114
• /
• 2023

The Christian Museum of Gongju Jeil Church was first built in 1931 and was largely damaged during the Korean War, but the walls and chimneys have been preserved. This building has a high architectural values in that the chapel was reconstructed in 1956, and maintains its original form through repair of damaged parts rather than new construction. The stained glass windows were as installed in 1979 and has a great significance in the Dalle de Verre method using lump glass. However, some of the stained glass damaged partially, such as various cracks and splits, and vertical and horizontal cracks in the joint fillers of supporting the colored glass. As the structural materials of the stained glass window, an iron frame and cement mortar filled with it were used, and corrosion of iron, cracking of mortar and granular decomposition appear partially due to weathering. In the joint fillers, the content of Ca and S is very high, indicating that gypsum were used as admixtures, and the gypsums grow in a rhombohedral and forms a bundle, which is investigated to have undergone recrystallization. As a result of modeling the ultrasonic velocity at the joint fillers, the left and right windows at the entrance show relatively weak in the range of 800 to 1,600m/s, and the lower right corner of the altar window and the upper left corner of the center window were also 1,000 to 1,800m/s, showing relatively low physical properties. And gypsums produced during the neutralization of lime mortar were detected in the joint fillers and contaminants on the surface. Such salts may cause damage to the joint material due to freezing and thawing, so appropriate preventive conservation is required. Also, since various damage types are complexly appearing in stained glass window and joint filler, customized conservation treatment should be reviewed through clinical tests.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.6-7
• /
• 2021

Vegetation development in rivers is one of the important issues not only in academic fields such as geomorphology, ecology, hydraulics, etc., but also in river management practices. The problem of river vegetation is directly connected to the harmony of conflicting values of flood management and ecosystem conservation. In Korea, since the 2000s, the issue of river vegetation and land formation has been continuously raised under various conditions, such as the regulating rivers downstream of the dams, the small eutrophicated tributary rivers, and the floodplain sites for the four major river projects. In this background, this study proposes a method for classifying the distribution of vegetation in rivers based on remote sensing data, and presents the results of applying this to the Naeseong Stream. The Naeseong Stream is a representative example of the river landscape that has changed due to vegetation development from 2014 to the latest. The remote sensing data used in the study are images of Sentinel 1 and 2 satellites, which is operated by the European Aerospace Administration (ESA), and provided by Google Earth Engine. For the ground truth, manually classified dataset on the surface of the Naeseong Stream in 2016 were used, where the area is divided into eight types including water, sand and herbaceous and woody vegetation. The classification method used a random forest classification technique, one of the machine learning algorithms. 1,000 samples were extracted from 10 pre-selected polygon regions, each half of them were used as training and verification data. The accuracy based on the verification data was found to be 82~85%. The model established through training was also applied to images from 2016 to 2020, and the process of changes in vegetation zones according to the year was presented. The technical limitations and improvement measures of this paper were considered. By providing quantitative information of the vegetation distribution, this technique is expected to be useful in practical management of vegetation such as thinning and rejuvenation of river vegetation as well as technical fields such as flood level calculation and flow-vegetation coupled modeling in rivers.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.28 no.2
• /
• pp.183-193
• /
• 1995

Assuming that fishing nets are porous structures to suck water into their mouth and then filtrate water out of them, the flow resistance N of nets with wall area S under the velicity v was taken by $R=kSv^2$, and the coefficient k was derived as $$k=c\;Re^{-m}(\frac{S_n}{S_m})n(\frac{S_n}{S})$$ where $R_e$ is the Reynolds' number, $S_m$ the area of net mouth, $S_n$ the total area of net projected to the plane perpendicular to the water flow. Then, the propriety of the above equation and the values of c, m and n were investigated by the experimental results on plane nettings carried out hitherto. The value of c and m were fixed respectively by $240(kg\cdot sec^2/m^4)$ and 0.1 when the representative size on $R_e$ was taken by the ratio k of the volume of bars to the area of meshes, i. e., $$\lambda={\frac{\pi\;d^2}{21\;sin\;2\varphi}$$ where d is the diameter of bars, 21 the mesh size, and 2n the angle between two adjacent bars. The value of n was larger than 1.0 as 1.2 because the wakes occurring at the knots and bars increased the resistance by obstructing the filtration of water through the meshes. In case in which the influence of $R_e$ was negligible, the value of $cR_e\;^{-m}$ became a constant distinguished by the regions of the attack angle $ \theta$ of nettings to the water flow, i. e., 100$(kg\cdot sec^2/m^4)\;in\;45^{\circ}<\theta \leq90^{\circ}\;and\;100(S_m/S)^{0.6}\;(kg\cdot sec^2/m^4)\;in\;0^{\circ}<\theta \leq45^{\circ}$. Thus, the coefficient $k(kg\cdot sec^2/m^4)$ of plane nettings could be obtained by utilizing the above values with $S_m\;and\;S_n$ given respectively by $$S_m=S\;sin\theta$$ and $$S_n=\frac{d}{I}\;\cdot\;\frac{\sqrt{1-cos^2\varphi cos^2\theta}} {sin\varphi\;cos\varphi} \cdot S$$ But, on the occasion of $\theta=0^{\circ}$ k was decided by the roughness of netting surface and so expressed as $$k=9(\frac{d}{I\;cos\varphi})^{0.8}$$ In these results, however, the values of c and m were regarded to be not sufficiently exact because they were obtained from insufficient data and the actual nets had no use for k at $\theta=0^{\circ}$. Therefore, the exact expression of $k(kg\cdotsec^2/m^4)$, for actual nets could De made in the case of no influence of $R_e$ as follows; $$k=100(\frac{S_n}{S_m})^{1.2}\;(\frac{S_m}{S})\;.\;for\;45^{\circ}<\theta \leq90^{\circ}$$, $$k=100(\frac{S_n}{S_m})^{1.2}\;(\frac{S_m}{S})^{1.6}\;.\;for\;0^{\circ}<\theta \leq45^{\circ}$$