• Journal of the Korean earth science society
• /
• v.36 no.1
• /
• pp.109-124
• /
• 2015

A novel approach, hybrid surface rainfall (KNU-HSR) technique developed by Kyungpook Natinal University, was utilized for improving the radar rainfall estimation. The KNU-HSR technique estimates radar rainfall at a 2D hybrid surface consistings of the lowest radar bins that is immune to ground clutter contaminations and significant beam blockage. Two HSR techniques, static and dynamic HSRs, were compared and evaluated in this study. Static HSR technique utilizes beam blockage map and ground clutter map to yield the hybrid surface whereas dynamic HSR technique additionally applies quality index map that are derived from the fuzzy logic algorithm for a quality control in real time. The performances of two HSRs were evaluated by correlation coefficient (CORR), total ratio (RATIO), mean bias (BIAS), normalized standard deviation (NSD), and mean relative error (MRE) for ten rain cases. Dynamic HSR (CORR=0.88, BIAS= $-0.24mm\;hr^{-1}$, NSD=0.41, MRE=37.6%) shows better performances than static HSR without correction of reflectivity calibration bias (CORR=0.87, BIAS= $-2.94mm\;hr^{-1}$, NSD=0.76, MRE=58.4%) for all skill scores. Dynamic HSR technique overestimates surface rainfall at near range whereas it underestimates rainfall at far ranges due to the effects of beam broadening and increasing the radar beam height. In terms of NSD and MRE, dynamic HSR shows the best results regardless of the distance from radar. Static HSR significantly overestimates a surface rainfall at weaker rainfall intensity. However, RATIO of dynamic HSR remains almost 1.0 for all ranges of rainfall intensity. After correcting system bias of reflectivity, NSD and MRE of dynamic HSR are improved by about 20 and 15%, respectively.

• Spatial Information Research
• /
• v.22 no.1
• /
• pp.55-63
• /
• 2014

Near surface air temperature data which are one of the essential factors in hydrology, meteorology and climatology, have drawn a substantial amount of attention from various academic domains and societies. Meteorological observations, however, have high spatio-temporal constraints with the limits in the number and distribution over the earth surface. To overcome such limits, many studies have sought to estimate the near surface air temperature from satellite image data at a regional or continental scale with simple regression methods. Alternatively, we applied various Kriging methods such as ordinary Kriging, universal Kriging, Cokriging, Regression Kriging in search of an optimal estimation method based on near surface air temperature data observed from automatic weather stations (AWS) in South Korea throughout 2010 (365 days) and MODIS land surface temperature (LST) data (MOD11A1, 365 images). Due to high spatial heterogeneity, auxiliary data have been also analyzed such as land cover, DEM (digital elevation model) to consider factors that can affect near surface air temperature. Prior to the main estimation, we calculated root mean square error (RMSE) of temperature differences from the 365-days LST and AWS data by season and landcover. The results show that the coefficient of variation (CV) of RMSE by season is 0.86, but the equivalent value of CV by landcover is 0.00746. Seasonal differences between LST and AWS data were greater than that those by landcover. Seasonal RMSE was the lowest in winter (3.72). The results from a linear regression analysis for examining the relationship among AWS, LST, and auxiliary data show that the coefficient of determination was the highest in winter (0.818) but the lowest in summer (0.078), thereby indicating a significant level of seasonal variation. Based on these results, we utilized a variety of Kriging techniques to estimate the surface temperature. The results of cross-validation in each Kriging model show that the measure of model accuracy was 1.71, 1.71, 1.848, and 1.630 for universal Kriging, ordinary Kriging, cokriging, and regression Kriging, respectively. The estimates from regression Kriging thus proved to be the most accurate among the Kriging methods compared.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.2
• /
• pp.87-95
• /
• 2016

This study examined the characteristics of the hourly demand of city gas in Korea and proposed multiple regression models to obtain precise estimates of the hourly demand of city gas. Forecasting the hourly demand of city gas with accuracy is essential in terms of safety and cost. If underestimated, the pipeline pressure needs to be increased sharply to meet the demand, when safety matters. In the opposite case, unnecessary inventory and operation costs are incurred. Data analysis showed that the hourly demand of city gas has a very high autocorrelation and that the 24-hour demand pattern of a day follows the previous 24-hour demand pattern of the same day. That is, there is a weekly cycle pattern. In addition, some conditions that temperature affects the hourly demand level were found. That is, the absolute value of the correlation coefficient between the hourly demand and temperature is about 0.853 on average, while the absolute value of the correlation coefficient on a specific day improves to 0.861 at worst and 0.965 at best. Based on this analysis, this paper proposes a multiple regression model incorporating the hourly demand ahead of 24 hours and the hourly demand ahead of 168 hours, and another multiple regression model with temperature as an additional independent variable. To show the performance of the proposed models, computational experiments were carried out using real data of the domestic city gas demand from 2009 to 2013. The test results showed that the first regression model exhibits a forecasting accuracy of MAPE (Mean Absolute Percentage Error) around 4.5% over the past five years from 2009 to 2013, while the second regression model exhibits 5.13% of MAPE for the same period.

• Journal of Korea Water Resources Association
• /
• v.55 no.10
• /
• pp.761-774
• /
• 2022

Accurate hydrologic prediction is essential to analyze the effects of drought, flood, and climate change on flow rates, water quality, and ecosystems. Disentangling the uncertainty of the hydrological model is one of the important issues in hydrology and water resources research. Hydrologic data assimilation (DA), a technique that updates the status or parameters of a hydrological model to produce the most likely estimates of the initial conditions of the model, is one of the ways to minimize uncertainty in hydrological simulations and improve predictive accuracy. In this study, the two ensemble-based sequential DA techniques, ensemble Kalman filter, and particle filter are comparatively analyzed for the daily discharge simulation at the Yongdam catchment using airGRdatassim. The results showed that the values of Kling-Gupta efficiency (KGE) were improved from 0.799 in the open loop simulation to 0.826 in the ensemble Kalman filter and to 0.933 in the particle filter. In addition, we analyzed the effects of hyper-parameters related to the data assimilation methods such as precipitation and potential evaporation forcing error parameters and selection of perturbed and updated states. For the case of forcing error conditions, the particle filter was superior to the ensemble in terms of the KGE index. The size of the optimal forcing noise was relatively smaller in the particle filter compared to the ensemble Kalman filter. In addition, with more state variables included in the updating step, performance of data assimilation improved, implicating that adequate selection of updating states can be considered as a hyper-parameter. The simulation experiments in this study implied that DA hyper-parameters needed to be carefully optimized to exploit the potential of DA methods.

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

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

• Journal of Radiation Protection and Research
• /
• v.23 no.3
• /
• pp.139-147
• /
• 1998

Purpose: Tissue inhomogeneity such as lung affects tumor dose as well as transmission dose in new concept of on-line dosimetry which estimates tumor dose from transmission dose using the new algorithm. This study was carried out to confirm accuracy of correction by tissue density in tumor dose estimation utilizing transmission dose. Methods: Cork phantom (CP, density $0.202\;gm/cm^3$) having similar density with lung parenchyme and polystyrene phantom (PP, density $1.040\;gm/cm^3$) having similar density with soft tissue were used. Dose measurement was carried out under condition simulating human chest. On simulating AP-PA irradiation, PPs with 3 cm thickness were placed above and below CP, which had thickness of 5, 10, and 20 cm. On simulating lateral irradiation, 6 cm thickness of PP was placed between two 10 cm thickness CPs additional 3 cm thick PP was placed to both lateral sides. 4, 6, and 10 MV x-ray were used. Field size was in the range of $3{\times}3$ cm through $20{\times}20$ cm, and phantom-chamber distance (PCD) was 10 to 50 cm. Above result was compared with another sets of data with equivalent thickness of PP which was corrected by density. Result: When transmission dose of PP was compared with equivalent thickness of CP which was corrected with density, the average error was 0.18 (${\pm}0.27$) % for 4 MV, 0.10 (${\pm}0.43$) % for 6 MV, and 0.33 (${\pm}0.30$) % for 10 MV with CP having thickness of 5 cm. When CP was 10 cm thick, the error was 0.23 (${\pm}0.73$) %, 0.05 (${\pm}0.57$) %, and 0.04 (${\pm}0.40$) %, while for 20 cm, error was 0.55 (${\pm}0.36$) %, 0.34 (${\pm}0.27$) %, and 0.34 (${\pm}0.18$) % for corresponding energy. With lateral irradiation model, difference was 1.15 (${\pm}1.86$) %, 0.90 (${\pm}1.43$) %, and 0.86 (${\pm}1.01$) % for corresponding energy. Relatively large difference was found in case of PCD having value of 10 cm. Omitting PCD with 10 cm, the difference was reduced to 0.47 (${\pm}$1.17) %, 0.42 (${\pm}$0.96) %, and 0.55 (${\pm}$0.77) % for corresponding energy. Conclusion When tissue inhomogeneity such as lung is in tract of x-ray beam, tumor dose could be calculated from transmission dose after correction utilizing tissue density.

• Journal of Sericultural and Entomological Science
• /
• v.8
• /
• pp.11-25
• /
• 1968

Various formulae for estimation of leaf production in mulberry trees were investigated and obtained. Four varieties of mulberry trees were used as the materials, and seven characters namely branch length. branch diameter, node number per branch, total branch weight, branch weight except leaves, leaf weight and leaf area, were studied. The formulae to estimate the leaf yield of mulberry trees are as follows: 1. Varietal differences were appeared in means, variances, standard devitations and standard errors of seven characters studied as shown in table 1. 2. Y$_1$=a$_1$X$_1$${\times}$P$_1$......(l) where Y$_1$ means yield per l0a by branch number and leaf weight determination. a$_1$.........leaf weight per branch. X$_1$.......branch number per plant. P$_1$........plant number per l0a. 3. Y$_2$=(a$_2$${\pm}$S. E.${\times}$X$_2$)+P$_1$.......(2) where Y$_2$ means leaf yield per l0a by branch length and leaf weight determination. a$_2$......leaf weight per meter of branch length. S. E. ......standard error. X$_2$....total branch length per plant. P$_1$........plant number per l0a as written above. 4. Y$_3$=(a$_3$${\pm}$S. E${\times}$X$_3$)${\times}$P$_1$.....(3) where Y$_3$ means of yield per l0a by branch diameter measurement. a$_3$.......leaf weight per 1cm of branch diameter. X$_3$......total branch diameter per plant. 5. Y$_4$=(a$_4$${\pm}$S. E.${\times}$X$_4$)P$_1$......(4) where Y$_4$ means leaf yield per 10a by node number determination. a$_4$.......leaf weight per node X$_4$.....total node number per plant. 6. Y$\sub$5/= {(a$\sub$5/${\pm}$S. E.${\times}$X$_2$)Kv}${\times}$P$_1$.......(5) where Y$\sub$5/ means leaf yield per l0a by branch length and leaf area measurement. a$\sub$5/......leaf area per 1 meter of branch length. K$\sub$v/......leaf weight per 100$\textrm{cm}^2$ of leaf area. 7. Y$\sub$6/={(X$_2$$\div$a$\sub$6/${\pm}$S. E.)}${\times}$K$\sub$v/${\times}$P$_1$......(6) where Y$\sub$6/ means leaf yield estimated by leaf area and branch length measurement. a$\sub$6/......branch length per l00$\textrm{cm}^2$ of leaf area. X$_2$, K$\sub$v/ and P$_1$ are written above. 8. Y$\sub$7/= {(a$\sub$7/${\pm}$S. E. ${\times}$X$_3$)}${\times}$K$\sub$v/${\times}$P$_1$.......(7) where Y$\sub$7/ means leaf yield estimates by branch diameter and leaf area measurement. a$\sub$7/......leaf area per lcm of branch diameter. X$_3$, K$\sub$v/ and P$_1$ are written above. 9. Y$\sub$8/= {(X$_3$$\div$a$\sub$8/${\pm}$S. E.)}${\times}$K$\sub$v/${\times}$P$_1$.......(8) where Y$\sub$8/ means leaf yield estimates by leaf area branch diameter. a$\sub$8/......branch diameter per l00$\textrm{cm}^2$ of leaf area. X$_3$, K$\sub$v/, P$_1$ are written above. 10. Y$\sub$9/= {(a$\sub$9/${\pm}$S. E.${\times}$X$_4$)${\times}$K$\sub$v/}${\times}$P$_1$......(9) where Y$\sub$7/ means leaf yield estimates by node number and leaf measurement. a$\sub$9/......leaf area per node of branch. X$_4$, K$\sub$v/, P$_1$ are written above. 11. Y$\sub$10/= {(X$_4$$\div$a$\sub$10/$\div$S. E.)${\times}$K$\sub$v/}${\times}$P$_1$.......(10) where Y$\sub$10/ means leaf yield estimates by leaf area and node number determination. a$\sub$10/.....node number per l00$\textrm{cm}^2$ of leaf area. X$_4$, K$\sub$v/, P$_1$ are written above. Among many estimation methods. estimation method by the branch is the better than the methods by the measurement of node number and branch diameter. Estimation method, by branch length and leaf area determination, by formulae (6), could be the best method to determine the leaf yield of mulberry trees without destroying the leaves and without weighting the leaves of mulberry trees.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.27 no.1
• /
• pp.1-12
• /
• 1991

An experiment has been carefully designed and performed to verify the theory for the echointergration technique of estimating the density of fish school by the use of steel spheres in a laboratory tank. The spheres used to simulate a fish school were randomly distributed throughout the insonified volume to produce the acoustic echoes similar to those scattered from real fish schools. The backscattered echoes were measured as a function of target density at tow frequencies of 50kHz and 200kHz. Data acquisition, processing and analysis were performed by means of the microcomputer-based sonar-echo processor including a FFT analyzer. Acoustic scattering characteristics of a 36cm mackerel was investigated by measuring fish echoes with frequencies ranging from 47.8kHz to 52.0kHz. The fluctuation of bottom echoes caused by the effects of fish-school attenuation and multiple scattering which occurred in dense aggregations of fishes was also examined by analyzing the echograms of sardine schools obtained by a 50kHz telesounder in the set-net's bagnet, and the echograms obtained by a scientific echo sounder of 50kHz in the East China Sea, respectively. The results obtained can be summarized as follows: 1. The measured and the calculated echo shapes on the steel sphere used to simulate a fish school were in close agreement. 2. The waveform and amplitude of echo signals by a mackerel without swimbladder fluctuated irregularly with the measuring frequency. 3. When a collection of 30 targets/m super(3) lied the shadow region behind another collection of 5 targets/m super(3), the mean losses in echo energy for the 30 targets/m super(3) were about -0.4dB at 50kHz and about -0.2dB at 200kHz, respectively. 4. In the echograms obtained in the East China Sea, the bottom echoes fluctuated remarkably when the dense aggregations of fish appeared between transducer and seabed. Especially, in the case of the echograms of sardine school obtained in a set-net's bagnet, the disappearance of bottom echoes and the lengthening of the echo trace by fish aggregations were observed. Then the mean density of the sardine school was estimated as 36 fish/m super(3). It suggests that when the distribution density of fishes in oceans is greater than this density, the effects of fish-school attenuation and multiple scattering must be taken into account as a possible source of error in fish abundance estimates. 5. The relationship between mean backscattering strength (, dB) and target density ($\rho$, No./m super(3)) were expressed by the equations: =-46.2+13.7 Log($\rho$) at 50kHz and =-43.9+13.4 Log($\rho$) at 200kHz. 6. The difference between the experimentally derived number and the actual number of targets gradually decreased with an increase in the target density and was within 20% when the density was 30 targets/m super(3). From these results, we concluded that when the number of targets in the insonified volume is large, the validity of the echo-integration technique of estimating the density of fish schools could be expected.

• Journal of the Korean Regional Science Association
• /
• v.34 no.2
• /
• pp.3-19
• /
• 2018

The objective of this study is to observe the efficiency of clinical performance on the supply and demand of medical resources in Korea. For the empirical analysis, we constructed the dataset on age standardized mortality rate, the number of physician, specialist, surgery, medical institution, ratio of general hospitals of 16 provinces in Korea from 2006 to 2013. The panel probability frontier model is employed as an analysis method and considered heteroscedasticity and autocorrelation of the error in panel data. In addition, the demographic and socioeconomic characteristics of the 16 provinces, unemployment rate, elderly population ratio, GRDP per capita, and ratio of hospitals in comparison to the general hospitals are used to find the effect on the technical efficiency of clinical performance on supply and demand of medical resources. The results are as follows. First, for the clinical performance, the supply side of human resources such as doctors and specialists and the demand side factors such as chronic illness clinic per unit population have a significant influence, respectively. Second, the technical efficiency of clinical performance on the supply and demand of medical resources of each input component was 59-70% in terms of clinical efficiency in each region. Third. estimates of technical efficiency of inputs that affect clinical performance showed a slight increase in all regions during the analysis period, but the increase trend decreased slightly. Fourth, the ratio of the elderly population and GRDP per capita have a positive influence on the technical efficiency of clinical performance on the supply and demand of medical resources. The difference of each efficiency by region is due to the regional differences of the input medical resources and the combination of them and the demographic and socioeconomic characteristics of the region. It is understood that the differences in technological efficiency due to the complexity of supply and demand of medical resources, demographic structure and economic difference affecting clinical performance by region are different.