• Geophysics and Geophysical Exploration
• /
• v.15 no.1
• /
• pp.16-22
• /
• 2012

In this study, we setup a ground-based scatterometer using an antenna of which the center frequency is 9.5 GHz (X-band), and measured radar backscatterings from air/ice and ice/water interfaces to extract ice thickness. Both of air/ice and ice/water interfaces make strong radar backscatterings and so we can clearly identify two peaks in measured data by scatterometer. By using the distance of two peaks and refractive index of ice, we confirmed that it is possible to measure ice thickness. Field survey was performed at the downstream of Jiam River flowing into Chuncheon Lake. We measured radar backscattering from river ice along a survey path and extracted ice thickness. The ice thickness map of the downstream of Jiam River was produced by using kriging which is one of well known interpolation methods. The ice thickness was about 50 cm along the mainstream while ice was thin as 30 ~ 40 cm at a fast-flowing meander. Ice thickness was particularly thinner at some locations than that of surrounding areas even in the mainstream region of constant flow. This was because of impurities in ice or artificially formed refrozen holes after fishing. We expect that this study helps to expand utilization field of X-band SAR and airborne scatterometer system.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.414-414
• /
• 2017

토양수분은 지면환경에서 일어나는 수문 및 에너지 순환을 이해하는 데 있어 중요한 기상인자이다. 토양수분 현장관측은 땅속에 매설된 센서에 의해 상당히 정확하게 이루어지만, 관측점 수가 충분치 않아 공간적 연속성을 확보하지 못하는 어려움이 존재한다. 이에 광역적 및 연속적 관측이 가능한 마이크로파 위성센서가 토양수분 정보 획득을 위한 보조수단으로서 그 중요성이 부각되고 있다. 마이크로파 위성센서는 구름 등 기상조건의 제약을 받지 않으며, 1978년 이래 현재까지 여러 위성에 의해 25 km 및 10 km 해상도의 전지구 토양수분자료가 생산되어 왔다. 마이크로파 센서를 이용한 토양수분자료는 동일지점에 대하여 하루 2회 정도 산출되므로 적절한 시간분해능을 가지지만, 공간해상도가 최고 10 km로서 지역규모의 수문분석에 적용하기에는 충분치 않다. 이러한 토양수분자료의 공간해상도 문제 해결을 위하여 다양한 지면환경요소를 활용한 통계적 다운스케일링이 대안으로 제시되었다. 최근의 선행연구들은 대부분 방정식을 이용한 결합모형을 통해 통계적 다운스케일링을 수행하였는데, 회귀식과 같은 선형결합뿐 아니라 신경망이나 기계학습 등의 비선형결합에서도, 불가피하게 발생할 수밖에 없는 잔차(residual)로 인하여 다운스케일링 전후의 공간분포 패턴이 달라져버리는 문제를 안고 있었다. 회귀분석에 잔차의 공간내삽을 결합시킨 회귀크리깅(regression kriging)은 잔차보정을 통해 이러한 문제를 해결함으로써 다운스케일링 전후의 공간분포 일관성을 보장하는 기법이다. 이 연구에서는 회귀크리깅을 이용하여 일자별 AMSR2(Advanced Microwave Scanning Radiometer 2) 토양수분 자료를 10 km에서 1 km 해상도로 다운스케일링하고, 다운스케일링 전후의 자료패턴 일관성을 평가한다. 지면온도(LST), 지면온도상승률(RR), 식생온도건조지수(TVDI)는 일자별로 DB를 구축하였고, 식생지수(NDVI), 수분지수(NDWI), 지면알베도(SA)는 8일 간격으로 DB를 구축하였다. 이러한 8일 간격의 자료를 일자별로 변환하기 위하여 큐빅스플라인(cubic spline)을 이용하여 시계열내삽을 수행하였다. 또한 상이한 공간해상도의 자료는 최근린법을 이용하여 다운스케일링 목표해상도인 1 km에 맞도록 변환하였다. 우선 저해상도 스케일에서 추정치를 산출하기 위해서는 저해상도 픽셀별로 이에 해당하는 복수의 고해상도 픽셀을 평균화하여 대응시켜야 하며, 이를 통해 6개의 설명변수(LST, RR, TVDI, NDVI, NDWI, SA)와 AMSR2 토양수분을 반응변수로 하는 다중회귀식을 도출하였다. 이식을 고해상도 스케일의 설명변수들에 적용하면 고해상도 토양수분 추정치가 산출되는데, 이때 추정치와 원자료의 차이에 해당하는 잔차에 대한 보정이 필요하다. 저해상도 스케일로 존재하는 잔차를 크리깅 공간내삽을 통해 고해상도로 변환한 후 이를 고해상도 추정치에 부가해주는 방식으로 잔차보정이 이루어짐으로써, 다운스케일링 전후의 자료패턴 일관성이 유지되는(r>0.95) 공간상세화된 토양수분 자료를 생산할 수 있다.

• Journal of Soil and Groundwater Environment
• /
• v.7 no.3
• /
• pp.3-17
• /
• 2002

Data analysis of groundwater monitoring wells and geostatistical methods are used to identify the local characteristics of sea water intrusion and the range of sea water intrusion at the southeastern coastal area of Busan, Korea. Rainfall and groundwater level of two monitoring wells show a linear correlation because of the direct groundwater recharge by the precipitation. However, rainfall and electric conductivity have the inverse relationship because of the increase of groundwater. Electric conductivity rapidly increased at 24m depth and exceeded 20,000$\mu\textrm{s}$/cm near 26m depth in the monitoring wells. The variations of groundwater level and electric conductivity show that the interface between sea water and fresh water tends to move upward when groundwater level goes down. In the cross correlation analysis, groundwater level versus rainfall represents the largest cross correlation coefficient in 0 time lag but the cross correlation coefficient of electric conductivity versus rainfall is the largest when the time lag is 9 days. This suggests that the fluctuations of groundwater level respond to rainfall in a short time, but the interface between sea water and fresh water respond very slow to rainfall. Horizontal extents of sea water intrusion are estimated to 14 m from the east of Line 1, and 25 m from the southeast end of Line 2 in the inversion of dipole-dipole profiling data of two survey lines. The data of VES by the Schulumberger array in May and July show lognormal distributions. In the kriged apparent resistivity and earth resistivity distributions, the resistivities of July are increased comparing to those of May. This reflects that the concentration of sea water in aquifer is reduced due to the increased groundwater recharge from the rainfall in June and July. In analyzing the vertical and horizontal apparent resistivities and earth resistivity distributions, the geostatistical methods are very useful to identify the variations of earth resistivity distributions at the coastal area.

• Journal of Korean Society of Forest Science
• /
• v.110 no.3
• /
• pp.322-340
• /
• 2021

Due to climate change and its consequential spring temperature rise, flowering time of Robinia pseudoacacia has advanced and a simultaneous blooming phenomenon occurred in different regions in South Korea. These changes in flowering time became a major crisis in the domestic beekeeping industry and the demand for accurate prediction of flowering time for R. pseudoacacia is increasing. In this study, we developed and compared performance of four different models predicting flowering time of R. pseudoacacia for the entire country: a Single Model for the country (SM), Modified Single Model (MSM) using correction factors derived from SM, Group Model (GM) estimating parameters for each region, and Local Model (LM) estimating parameters for each site. To achieve this goal, the bloom date data observed at 26 points across the country for the past 12 years (2006-2017) and daily temperature data were used. As a result, bloom dates for the north central region, where spring temperature increase was more than two-fold higher than southern regions, have advanced and the differences compared with the southwest region decreased by 0.7098 days per year (p-value=0.0417). Model comparisons showed MSM and LM performed better than the other models, as shown by 24% and 15% lower RMSE than SM, respectively. Furthermore, validation with 16 additional sites for 4 years revealed co-krigging of LM showed better performance than expansion of MSM for the entire nation (RMSE: p-value=0.0118, Bias: p-value=0.0471). This study improved predictions of bloom dates for R. pseudoacacia and proposed methods for reliable expansion to the entire nation.

• Korean Journal of Soil Science and Fertilizer
• /
• v.22 no.3
• /
• pp.163-172
• /
• 1989

In order to study spatial variability of soil physical properties accross the land-partitioned boundary on Hwadong silt clay loam soil (Fine clayey, mixed, mesic family of Aquic Hapludalfs), all measured data were analyzed by means of kriging, fractile diagram, smooth frequency distribution, and autocorrelation. Sampling for soil particle size distribution analysis was made at 225 intersections of $15{\times}15$ grid with 10m interval. Field capacity, bulk density and saturated hydraulic conductivity were measured in situ at 594 intersections of $33{\times}18$ grid with 2.5m interval in only $6,000m^2$ reselected from $22,500m^2$ of sampling area for particle size distribution analysis. Sampled or measured soil depths were 0 to 10cm 25 to 35cm and 50 to 60cm at each intersections. The results are summarized as follows: 1. The coefficient of variance (CV) of various physical properties ranges from 4.8 to 128.8%. Saturated hydraulic conductivity is classified into the high variation group with CV greater than 100%, while the low variation group with CV smaller than 10% consists of bulk density. Other properties belong to the medium variation group with CV between 10 and 100%. 2. The appropriate number of soil samples for the determination of various physical properties with error smaller than 10% are calculated as one for bulk density, six for field moisture capacity, 16 for silt, 19 for clay, 69 for sand and 686 for saturated hydraulic conductivity. 3. Smooth frequency distribution and fractile diagram show that saturated hydraulic conductivity is in lognormal distribution while other physical properties are in normal distribution. 4. Serial correlation analysis reveals that the soil physical properties have spatial dependence between two nearest neighbouring grid points. Autocorrelation analysis of physical properties measured between the serial grid points in the direction of south to north following section boundary shows that the zone of influence showing stationarity ranges from 7.5 to 40m. In the direction of east to west across section boundary, the autocorrelogram of many physical properties shows peaks with the periodic interval of 30m, which are similar to the partitioned land width. This reveals that the land-partitioned boundary causes soil variability.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.38 no.2
• /
• pp.357-367
• /
• 2018

Recent earthquakes in Gyeongju and Pohang have raised interest in liquefaction in South Korea. Liquefaction, which is a phenomenon that excessive pore pressure is generated and the shear strength of soil is decreased by repeated loads such as earthquakes, causes severe problems such as ground subsidence and overturning of structures. Therefore, it is necessary to identify and prepare for the possibility of liquefaction in advance. In general, the possibility of liquefaction is quantitatively assessed using the Liquefaction Potential Index (LPI), but it takes a lot of time and effort for performing site response analysis which is essential for the liquefaction evaluation. In this study, a simple method to evaluate the liquefaction potential without executing the site response analysis in a downtown area with a lot of borehole data was proposed. In this simple method, the correlation between the average shear wave velocity of the target location ground and the LPI divided by thickness of liquefiable layer was established. And the applicable correlation equation for various rock outcrop accelerations were derived. Using the 104 boreholes information in Seoul, the correlation equation between LPI and the shear wave velocity (ground water level: 0m, 1m, 2m, 3m) is obtained and the possibility of liquefaction occurrence in Seoul and Gyeongju is evaluated. The applicability of the proposed simple method was verified by comparing the LPI values calculated from the correlation equation and the LPI values derived using the existing site response analysis. Finally, the distribution map of LPI calculated from the correlation was drawn using Kriging, a geostatistical technique.

• Economic and Environmental Geology
• /
• v.38 no.4 s.173
• /
• pp.435-450
• /
• 2005

The purpose of the study is to analyze the hydrogeochemical characteristics by multivariate statistical method, to interpret the hydrogeochemical processes for the new variables calculated from principal components analysis (PCA), and to infer the groundwater flow and circulation mechanism by applying the geostatistical methods for each element and principal component. Chloride and nitrate are the most influencing components for groundwater quality, and the contents of $NO_3$ increased by the input of agricultural activities show the largest variation. The results of PCA, a multivariate statistical method, show that the first three principal components explain $73.9\%$ of the total variance. PC1 indicates the increase of dissolved ions, PC2 is related with the dissolution of carbonate minerals and nitrate contamination, and PC3 shows the effect of cation exchange process and silicate mineral dissolution. From the results of experimental semivariogram, the components of groundwater are divided into two groups: one group includes electrical conductivity (EC), Cl, Na, and $NO_3$, and the other includes $HCO_3,\;SiO_2,$ Ca, and Sr. The results for spatial distribution of groundwater components showed that EC, Cl, and Na increased with approaching the coastal line and nitrate has close relationship with the presence of agricultural land. These components are also correlated with the topographic features reflecting the groundwater recharge effect. The kriging analysis by using principal components shows that PC 1 has the different spatial distribution of Cl, Na, and EC, possibly due to the influence of pH, Ca, Sr, and $HCO_3$ for PC1. It was considered that the linear anomaly zone of PC2 in western area was caused by the dissolution of carbonate mineral. Consequently, the application of multivariate and geostatistical methods for groundwater in the study area is very useful for determining the quantitative analysis of water quality data and the characteristics of spatial distribution.

• Korean Journal of Soil Science and Fertilizer
• /
• v.22 no.4
• /
• pp.257-264
• /
• 1989

In order to study spatial variability of soil chemical properties across the land-partitioned boundary on Hwadong silt clay loam soil (Fine clayey, mixed, mesic family of Aquic Hapludalfs) in the experimental fie ld of the wheat and Barley Research Institute in Suwon, all measured data were analyzed by means of kriging, fractile diagram, smooth frequency distribution, and autocorrelation. Sampling for soil chemical property analysis was made at 225 intersections of 15x 15 grid with 10m interval from three soil depths (0-10cm, 25-35cm, 50-60cm) in the seven patitioned fields. 1. The coefficient of variance (CV) of various chemical properties varied from 5.4 to 72.7%. Soil pH was classified into the low variation group with CV smaller than 10%, while the other chemical properties belonged to the medium variation group with C.V. between 10 and 100% 2. The approximate number of soil samples for the determination of various chemical properties with error smaller than 10% were two for pH, ten for CEC, 15 for exchangeable Ca, 32 for total nitrogen content, 39 for exchangeable Mg, 40 for exchangeable K, 61 for exchangeable Na, 82 for organic matter content, 212 for available phosphate,. 3. Smooth frequency distribution and fractile diagram showed that available phosphate was in log-normal distribution while others were in normal distribution. 4. Serial correlation analysis revaled that the soil chemical properties had spatial dependence between two nearest neighbouring grid points. Autocorrelation analysis of chemcial properties measured between the serial grid points in the direction of south to north following land-partitioned boundary showed that the zone of influence showing stationarity ranged from 20 to 50m. In the direction of east to west accross land-partitioned boundary, the autocorrelogram of many chemical properies showed peaks with the periodic interval of 30m, which were similar to the partitioned land width. This reveals that the land-partitioned boundary causes soil variability.