• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.26 no.1
• /
• pp.11-36
• /
• 2021

The solar annual (Sa) and semiannual (Ssa) tides account for much of the non-uniform annual and seasonal variability observed in sea levels. These non-equilibrium tides depend on atmospheric variations, forced by changes in the Sun's distance and declination, as well as on hydrographic conditions. Here we employ tidal harmonic analyses to calculate Sa and Ssa harmonic constants for 21 Korean coastal tidal stations (TS), operated by the Korea Hydrographic and Oceanographic Agency. We used 19 year-long (1999 to 2017) 1 hr-interval sea level records from each site, and used two conventional harmonic analysis (HA) programs (Task2K and UTide). The stability of Sa harmonic constants was estimated with respect to starting date and record length of the data, and we examined the spatial distribution of the calculated Sa and Ssa harmonic constants. HA was performed on Incheon TS (ITS) records using 369-day subsets; the first start date was January 1, 1999, the subsequent data subset starting 24 hours later, and so on up until the final start date was December 27, 2017. Variations in the Sa constants produced by the two HA packages had similar magnitudes and start date sensitivity. Results from the two HA packages had a large difference in phase lag (about 78°) but relatively small amplitude (<1 cm) difference. The phase lag difference occurred in large part since Task2K excludes the perihelion astronomical variable. Sensitivity of the ITS Sa constants to data record length (i.e., 1, 2, 3, 5, 9, and 19 years) was also tested to determine the data length needed to yield stable Sa results. HA results revealed that 5 to 9 year sea level records could estimate Sa harmonic constants with relatively small error, while the best results are produced using 19 year-long records. As noted earlier, Sa amplitudes vary with regional hydrographic and atmospheric conditions. Sa amplitudes at the twenty one TS ranged from 15.0 to 18.6 cm, 10.7 to 17.5 cm, and 10.5 to 13.0 cm, along the west coast, south coast including Jejudo, and east coast including Ulleungdo, respectively. Except at Ulleungdo, it was found that the Ssa constituent contributes to produce asymmetric seasonal sea level variation and it delays (hastens) the highest (lowest) sea levels. Comparisons between monthly mean, air-pressure adjusted, and steric sea level variations revealed that year-to-year and asymmetric seasonal variations in sea levels were largely produced by steric sea level variation and inverted barometer effect.

• Korean Journal of Remote Sensing
• /
• v.38 no.6_1
• /
• pp.1285-1300
• /
• 2022

Turbidity, the measure of the cloudiness of water, is used as an important index for water quality management. The turbidity can vary greatly in small river systems, which affects water quality in national rivers. Therefore, the generation of high-resolution spatial information on turbidity is very important. In this study, a turbidity retrieval model using the Korea Multi-Purpose Satellite-3 and -3A (KOMPSAT-3/3A) images was developed for high-resolution turbidity mapping of Han River system based on eXtreme Gradient Boosting (XGBoost) algorithm. To this end, the top of atmosphere (TOA) spectral reflectance was calculated from a total of 24 KOMPSAT-3/3A images and 150 Landsat-8 images. The Landsat-8 TOA spectral reflectance was cross-calibrated to the KOMPSAT-3/3A bands. The turbidity measured by the National Water Quality Monitoring Network was used as a reference dataset, and as input variables, the TOA spectral reflectance at the locations of in situ turbidity measurement, the spectral indices (the normalized difference vegetation index, normalized difference water index, and normalized difference turbidity index), and the Moderate Resolution Imaging Spectroradiometer (MODIS)-derived atmospheric products(the atmospheric optical thickness, water vapor, and ozone) were used. Furthermore, by analyzing the KOMPSAT-3/3A TOA spectral reflectance of different turbidities, a new spectral index, new normalized difference turbidity index (nNDTI), was proposed, and it was added as an input variable to the turbidity retrieval model. The XGBoost model showed excellent performance for the retrieval of turbidity with a root mean square error (RMSE) of 2.70 NTU and a normalized RMSE (NRMSE) of 14.70% compared to in situ turbidity, in which the nNDTI proposed in this study was used as the most important variable. The developed turbidity retrieval model was applied to the KOMPSAT-3/3A images to map high-resolution river turbidity, and it was possible to analyze the spatiotemporal variations of turbidity. Through this study, we could confirm that the KOMPSAT-3/3A images are very useful for retrieving high-resolution and accurate spatial information on the river turbidity.

• Journal of Korean Society of Forest Science
• /
• v.111 no.4
• /
• pp.622-629
• /
• 2022

The study aim was to derive a stem taper equation for Phyllostachys pubescens, a type of bamboo in South Korea, and to develop a stem volume table. To derive the stem taper equation, three stem taper models (Max & Burkhart, Kozak, and Lee) were used. Since bamboo stalks are hollow because of its woody characteristics, the outer and inner diameters of the tree were calculated, and connecting them enabled estimating the tree curves. The results of the three equations for estimating the outer and inner diameters led to selection of the Kozak model for determining the optimal stem taper because it had the highest fitness index and lowest error and bias. We used the Kozak model to estimate the diameter of Phyllostachys pubescens by stem height, which proved optimal, and drew the stem curve. After checking the residual degree in the stem taper equation, all residuals were distributed around "0", which proved the suitability of the equation. To calculate the stem volume of Phyllostachys pubescens, a rotating cube was created by rotating the stem curve with the outer diameter at 360°, and the volume was calculated by applying Smalian's method. The volume of Phyllostachys pubescens was calculated by deducting the inner diameter calculated volume from the outer diameter calculated volume. The volume of Phyllostachys pubescens was only 20~30% of the volume of Larix kaempferi, which is a general species. However, considering the current trees/ha of Phyllostachys pubescens and the amount of bamboo shoots generated every year, the individual tree volume was predicted to be small, but the volume/ha was not very different or perhaps more. The significance of this study is the stem taper equation and stem volume table for Phyllostachys pubescens developed for the first time in South Korea. The results are expected to be used as basic data for bamboo trading that is in increasing public and industrial demand and carbon absorption estimation.

• The Journal of Engineering Geology
• /
• v.32 no.4
• /
• pp.611-625
• /
• 2022

One of the applications of geomagnetic paleo-secular variation (PSV) is the age dating of archeological remains (i.e., the archeomagnetic dating technique). This application requires the local model of PSV that reflects non-dipole fields with regional differences. Until now, the tentative Korean paleosecular variation (t-KPSV) calculated based on JPSV (SW Japanese PSV) has been applied as a reference curve for individual archeomagnetic directions in Korea. However, it is less reliable due to regional differences in the non-dipole magnetic field. Here, we present PSV curves for AD 1 to 600, corresponding to the Korean Three Kingdoms (including the Proto Three Kingdoms) Period, using the results of archeomagnetic studies in the Korean Peninsula and published research data. Then we compare our PSV with the global geomagnetic prediction model and t-KPSV. A total of 49 reliable archeomagnetic directional data from 16 regions were compiled for our PSV. In detail, each data showed statistical consistency (N > 6, 𝛼₉₅ < 7.8°, and k > 57.8) and had radiocarbon or archeological ages in the range of AD 1 to 600 years with less than ±200 years error range. The compiled PSV for the initial six centuries (KPSV0.6k) showed declination and inclination in the range of 341.7° to 20.1° and 43.5° to 60.3°, respectively. Compared to the t-KPSV, our curve revealed different variation patterns both in declination and inclination. On the other hand, KPSV0.6k and global geomagnetic prediction models (ARCH3K.1, CALS3K.4, and SED3K.1) revealed consistent variation trends during the first six centennials. In particular, the ARCH3K.1 showed the best fitting with our KPSV0.6k. These results indicate that contribution of the non-dipole field to Korea and Japan is quite different, despite their geographical proximity. Moreover, the compilation of archeomagnetic data from the Korea territory is essential to build a reliable PSV curve for an age dating tool. Lastly, we double-check the reliability of our KPSV0.6k by showing a good fitting of newly acquired age-controlled archeomagnetic data on our curve.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.27 no.4
• /
• pp.194-210
• /
• 2022

The monthly inventory of dissolved inorganic carbon (C_T) and its fluxes were simulated using a box-model for the southeastern Yellow Sea, bordering the northern East China Sea. The monthly C_T data was constructed by combining the observed data representing four seasons with the data adopted from the recent publications. A 2-box-model of the surface and deep layers was used, assuming that the annual C_T inventory was at the steady state and its fluctuations due to the advection in the surface box were negligible. Results of the simulation point out that the monthly C_T inventory variation between the surface and deep box was driven primarily by the mixing flux due to the variation of the mixed layer depth, on the scale of -40~35 mol C m^-2 month^-1. The air to sea CO₂ flux was about 2 mol C m^-2 yr^-1 and was lower than 1/100 of the mixing flux. The biological pump flux estimated magnitude, in the range of 4-5 mol C m^-2 yr^-1, is about half the in situ measurement value reported. The C_T inventory of the water column was maximum in April, when mixing by cooling ceases, and decreases slightly throughout the stratified period. Therefore, the total C_T inventory is larger in the stratified period than that of the mixing period. In order to maintain a steady state, 18 mol C m^-2 yr^-1 (= 216 g C m^-2 yr^-1), the difference between the maximum and minimum monthly CT inventory, should be transported out to the East China Sea. Extrapolating this flux over the entire southern Yellow Sea boundary yields 4 × 10⁹ g C yr^-1. Conceptually this flux is equivalent to the proposed continental shelf pump. Since this flux must go through the vast shelf area of the East China Sea before it joins the open Pacific waters the actual contribution as a continental shelf pump would be significantly lower than reported value. Although errors accompanied the simple box model simulation imposed by the paucity of data and assumptions are considerably large, nevertheless it was possible to constrain the relative contribution among the major fluxes and their range that caused the C_T inventory variations, and was able to suggest recommendations for the future studies.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.25 no.1
• /
• pp.28-36
• /
• 2023

Solar radiation that is measured at relatively small number of weather stations is one of key inputs to crop models for estimation of crop productivity. Solar radiation products derived from GK-2A and Himawari 8 satellite data have become available, which would allow for preparation of input data to crop models, especially for assessment of crop productivity under an agrivoltaic system where crop and power can be produced at the same time. The objective of this study was to compare the degree of agreement between the solar radiation products obtained from those satellite data. The sub hourly products for solar radiation were collected to prepare their daily summary for the period from May to October in 2020 during which both satellite products for solar radiation were available. Root mean square error (RMSE) and its normalized error (NRMSE) were determined for daily sum of solar radiation. The cumulative values of solar radiation for the study period were also compared to represent the impact of the errors for those products on crop growth simulations. It was found that the data product from the Himawari 8 satellite tended to have smaller values of RMSE and NRMSE than that from the GK-2A satellite. The Himawari 8 satellite product had smaller errors at a large number of weather stations when the cumulative solar radiation was compared with the measurements. This suggests that the use of Himawari 8 satellite products would cause less uncertainty than that of GK2-A products for estimation of crop yield. This merits further studies to apply the Himawari 8 satellites to estimation of solar power generation as well as crop yield under an agrivoltaic system.

• Journal of Korean Society of Forest Science
• /
• v.112 no.2
• /
• pp.127-135
• /
• 2023

This study aimed to develop biomass allometric equations and estimate carbon emission factors, such as the wood density, biomass-expansion factor, and root-to-shoot ratio, for Platanus occidentalis and Metasequoia glyptostroboides planted in urban areas. Twenty M. glyptostroboides and 25 P. occidentalis trees were harvested, and the dry weights and stem volumes of stems, branches, leaves, and roots (>5 mm) were measured. The wood densities of M. glyptostroboides and P. occidentalis were 0.293 ± 0.008 g cm^-3 and 0.509 ± 0.018 g cm^-3, and the biomass-expansion factors were 1.738 ± 0.031 and 1.561 ± 0.035. The root-to-shoot ratios were 0.446 ± 0.009 and 0.402 ± 0.012. The uncertainty tests (coefficient of variation, %) gave 2.8% and 3.5% values for wood density, 1.8% and 2.3% for biomass-expansion factor, and 2.1% and 2.9% for root-to-shoot ratio, respectively. Among the developed allometric equations, Model I using the diameter at breast height (DBH) was suitable. The allometric equations of M. glyptostroboides and P. occidentalis above ground were y = 1.679 (DBH)^1.315 and y = 0.505 (DBH)^1.896, and the allometric equations of the root and total were y = 0.746 (DBH)^1.315, y = 0.301 (DBH)^1.751, y = 2.422 (DBH)^1.316, and y = 0.787 (DBH)^1.858. If the carbon-emission factors of this study and biomass allometric equations of the three developed models are used to estimate the carbon storage and biomass of urban forests, errors caused by not considering the use of fixed factors and the environmental differences can be reduced.

• Journal of Korean Society of Forest Science
• /
• v.112 no.4
• /
• pp.417-425
• /
• 2023

The aim of this study was to derive stem taper equations for Quercus acuta, one of main evergreen broad-leaved tree species found in warm temperate regions, and to prepare a stem volume table using those stem taper equations. A total of 688 individual trees were used in the analysis, which were collected from Jeonnam-do, Gyeongnam-do, and Jeju-do. The stem taper models applied to derive the stem curve pattern were the Max and Burkhart, Kozak, and Lee models. Among the three stem taper models, the best explanation of the stem curve shape of Q. acuta was found to be given by the Kozak model, which showed a fitness index of 0.9583, bias of 0.0352, percentage of estimated standard error of 1.1439, and mean absolute deviation of 0.6751. Thus, the stem taper of Q. acuta was estimated using the Kozak model. Moreover,thestemvolumecalculationwasperforme d by applying the Smalian formula to the diameter and height of each stem interval. In addition, an analysis of variance (ANOVA) was conducted to compare the two existing Q. acuta stem volume tables (2007 and 2010) and the newly created stem volume table (2023). This analysis revealed that the stem volume table constructed in the Wando region in 2007 included about twice as much as the stem volume tables constructed in 2010 and 2023. The stem volume table (2023) developed in this study is not only based on the regional collection range and number of utilized trees but also on a sound scientific basis. Therefore, it can be used at the national level as an official stem volume table for Q. acuta.

• 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.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.25 no.3
• /
• pp.129-141
• /
• 2023

Crop models have been used to predict yield under diverse environmental and cultivation conditions, which can be used to support decisions on the management of forage crop. Cultivar parameters are one of required inputs to crop models in order to represent genetic properties for a given forage cultivar. The objectives of this study were to compare calibration and ensemble approaches in order to minimize the uncertainty of crop yield estimates using the SIMPLE crop model. Cultivar parameters were calibrated using Log-likelihood (LL) and Generic Composite Similarity Measure (GCSM) as an objective function for Metropolis-Hastings (MH) algorithm. In total, 20 sets of cultivar parameters were generated for each method. Two types of ensemble approach. First type of ensemble approach was the average of model outputs (E_em), using individual parameters. The second ensemble approach was model output (E_pm) of cultivar parameter obtained by averaging given 20 sets of parameters. Comparison was done for each cultivar and for each error calculation methods. 'Jowoo' and 'Yeongwoo', which are forage rice cultivars used in Korea, were subject to the parameter calibration. Yield data were obtained from experiment fields at Suwon, Jeonju, Naju and I ksan. Data for 2013, 2014 and 2016 were used for parameter calibration. For validation, yield data reported from 2016 to 2018 at Suwon was used. Initial calibration indicated that genetic coefficients obtained by LL were distributed in a narrower range than coefficients obtained by GCSM. A two-sample t-test was performed to compare between different methods of ensemble approaches and no significant difference was found between them. Uncertainty of GCSM can be neutralized by adjusting the acceptance probability. The other ensemble method (E_pm) indicates that the uncertainty can be reduced with less computation using ensemble approach.