• Title/Summary/Keyword: 평균하부량

Search Result 125, Processing Time 0.035 seconds

Structure and Variation of Tidal Flat Temperature in Gomso Bay, West Coast of Korea (서해안 곰소만 갯벌 온도의 구조 및 변화)

  • Lee, Sang-Ho;Cho, Yang-Ki;You, Kwang-Woo;Kim, Young-Gon;Choi, Hyun-Yong
    • The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
    • /
    • v.10 no.1
    • /
    • pp.100-112
    • /
    • 2005
  • Soil temperature was measured from the surface to 40 cm depth at three stations with different heights in tidal flat of Gomso Bay, west coast of Korea, for one month in every season 2004 to examine the thermal structure and the variation. Mean temperature in surface layer was higher in summer and lower in winter than in lower layer, reflecting the seasonal variation of vertically propagating structure of temperature by heating and cooling from the tidal flat surface. Standard deviation of temperature decreased from the surface to lower layer. Periodic variations of solar radiation energy and tide mainly caused short term variation of soil temperature, which was also intermittently influenced by precipitation and wind. Time series analysis showed the power spectral energy peaks at the periods of 24, 12 and 8 hours, and the strongest peak appeared at 24 hour period. These peaks can be interpreted as temperature waves forced by variations of solar radiation, diurnal tide and interaction of both variations, respectively. EOF analysis showed that the first and the second modes resolved 96% of variation of vertical temperature structure. The first mode was interpreted as the heating antl cooling from tidal flat surface and the second mode as the effect of phase lag produced by temperature wave propagation in the soil. The phase of heat transfer by 24 hour period wave, analyzed by cross spectrum, showed that mean phase difference of the temperature wave increased almost linearly with the soil depth. The time lags by the phase difference from surface to 10, 20 and 40cm were 3.2,6.5 and 9.8 hours, respectively. Vertical thermal diffusivity of temperature wave of 24 hour period was estimated using one dimensional thermal diffusion model. Average diffusivity over the soil depths and seasons resulted in $0.70{\times}10^{-6}m^2/s$ at the middle station and $0.57{\times}10^{-6}m^2/s$ at the lowest station. The depth-averaged diffusivity was large in spring and small in summer and the seasonal mean diffusivity vertically increased from 2 cm to 10 cm and decreased from 10 cm to 40 cm. Thermal propagation speeds were estimated by $8.75{\times}10^{-4}cm/s,\;3.8{\times}10{-4}cm/s,\;and\;1.7{\times}10^{-4}cm/s$ from 2 cm to 10 cm, 20 cm and 40 cm, respectively, indicating the speed reduction with depth increasing from the surface.

The Geochemical Characteristics and Environmental Factors on the Marine Shellfish Farm in Namhae-po Tidal Flat of Taean (태안 남해포 갯벌 패류양식해역의 환경특성)

  • Choi, Yoon Seok;Park, Kwang Jae;Yoon, Sang Pil;Chung, Sang Ok;An, Kyoung Ho;Song, Jae Hee
    • The Korean Journal of Malacology
    • /
    • v.29 no.1
    • /
    • pp.51-63
    • /
    • 2013
  • To assess the effect of environmental factors on the sustainability of cultured production shellfish, we investigated the habitat characteristics of tidal flat (Namhae-po in Taean). We measured the physiochemical parameters (temperature, salanity, pH, dissolved oxygen and nutrients) and the geochemical characteristics (chemical oxygen demand, ignition loss, C/N ratio and C/S ratio). Surface sediments were collected from several site of tidal flat to examine the geochemical characteristics of both the benthic environment and heavy metal pollution. The grain size for research area of tidal flat were similar at the ratio of silt and clay in comparison with the other site of it. The C/N ratio was more than 5.0, reflecting the range arising from the mix of marine organism and organic matter. The C/S ratio (about 2.8) showed that survey area had anoxic or sub-anoxic bottom conditions. The enrichment factor (Ef) and index of accumulation rate (Igeo) of the metals showed that those research areas can be classified as heavily polluted, heavily to moderately polluted, or more or less unpolluted, respectively. Adult surf clam (Mactra veneriformis) density was highest at St. 2 (middle part of the Namhae-po), on the other hand, surf clam spat density was highest at St. 3 (lower part of the Namhae-po). Heavy rain, terrigenous suspended clay with fresh water from neighboring agricultural land, and severe high air temperature during summer could be thought as detrimental causes of spat and adult mortality in Namhae-po tidal flat. We suggested that the growth of shellfish in the tidal flat was effected by the various environmental conditions, so an improvement in the cultured method was needed.

A STUDY ON THE IONOSPHERE AND THERMOSPHERE INTERACTION BASED ON NCAR-TIEGCM: DEPENDENCE OF THE INTERPLANETARY MAGNETIC FIELD (IMF) ON THE MOMENTUM FORCING IN THE HIGH-LATITUDE LOWER THERMOSPHERE (NCAR-TIEGCM을 이용한 이온권과 열권의 상호작용 연구: 행성간 자기장(IMF)에 따른 고위도 하부 열권의 운동량 강제에 대한 연구)

  • Kwak, Young-Sil;Richmond, Arthur D.;Ahn, Byung-Ho;Won, Young-In
    • Journal of Astronomy and Space Sciences
    • /
    • v.22 no.2
    • /
    • pp.147-174
    • /
    • 2005
  • To understand the physical processes that control the high-latitude lower thermospheric dynamics, we quantify the forces that are mainly responsible for maintaining the high-latitude lower thermospheric wind system with the aid of the National Center for Atmospheric Research Thermosphere-Ionosphere Electrodynamics General Circulation Model (NCAR-TIEGCM). Momentum forcing is statistically analyzed in magnetic coordinates, and its behavior with respect to the magnitude and orientation of the interplanetary magnetic field (IMF) is further examined. By subtracting the values with zero IMF from those with non-zero IMF, we obtained the difference winds and forces in the high-latitude 1ower thermosphere(<180 km). They show a simple structure over the polar cap and auroral regions for positive($B_y$ > 0.8|$\overline{B}_z$ |) or negative($B_y$ < -0.8|$\overline{B}_z$|) IMF-$\overline{B}_y$ conditions, with maximum values appearing around -80$^{\circ}$ magnetic latitude. Difference winds and difference forces for negative and positive $\overline{B}_y$ have an opposite sign and similar strength each other. For positive($B_z$ > 0.3125|$\overline{B}_y$|) or negative($B_z$ < -0.3125|$\overline{B}_y$|) IMF-$\overline{B}_z$ conditions the difference winds and difference forces are noted to subauroral latitudes. Difference winds and difference forces for negative $\overline{B}_z$ have an opposite sign to positive $\overline{B}_z$ condition. Those for negative $\overline{B}_z$ are stronger than those for positive indicating that negative $\overline{B}_z$ has a stronger effect on the winds and momentum forces than does positive $\overline{B}_z$ At higher altitudes(>125 km) the primary forces that determine the variations of tile neutral winds are the pressure gradient, Coriolis and rotational Pedersen ion drag forces; however, at various locations and times significant contributions can be made by the horizontal advection force. On the other hand, at lower altitudes(108-125 km) the pressure gradient, Coriolis and non-rotational Hall ion drag forces determine the variations of the neutral winds. At lower altitudes(<108 km) it tends to generate a geostrophic motion with the balance between the pressure gradient and Coriolis forces. The northward component of IMF By-dependent average momentum forces act more significantly on the neutral motion except for the ion drag. At lower altitudes(108-425 km) for negative IMF-$\overline{B}_y$ condition the ion drag force tends to generate a warm clockwise circulation with downward vertical motion associated with the adiabatic compress heating in the polar cap region. For positive IMF-$\overline{B}_y$ condition it tends to generate a cold anticlockwise circulation with upward vertical motion associated with the adiabatic expansion cooling in the polar cap region. For negative IMF-$\overline{B}_z$ the ion drag force tends to generate a cold anticlockwise circulation with upward vertical motion in the dawn sector. For positive IMF-$\overline{B}_z$ it tends to generate a warm clockwise circulation with downward vertical motion in the dawn sector.

Recent updated diagnostic methods for esophageal motility disorders (식도의 운동장애에 관한 최신지견)

  • Yoon, Seok-Hwan
    • Journal of radiological science and technology
    • /
    • v.27 no.4
    • /
    • pp.11-16
    • /
    • 2004
  • Classification of esophageal motility disorders not yet finalized and is still ongoing as the new disorders are reported, and the existing classification is changed or removed. In terms of radiology, the primary peristalsis does not exist, and the lower end of the esophagus show the smooth, tapered, beak-like appearance. The esophageal motility disorder, which mostly occurs in the smooth muscle area, show the symptoms of reduction or loss (hypomotility) or abnormal increase (hypermotility) of peristalsis of the esophagus. It is important to understand the anatomy and physiology of the esophagus for the appropriate radiological method and diagnosis. Furthermore, the symptom of the patient and the manometry finding must be closely referred for the radiological diagnosis. The lower esophageal sphincter can be normally functioning and open completely as the food moves lower. Sperandio M et al. argues that the name diffuse esophageal spasm must be changed to distal esophageal spasm (DES) as most of the spasm occurs in the distal esophagus, composed of the smooth muscle. According to Ott et al., usefulness of barium method for diagnosing the esophageal motility disorder is Achalasia 95%, DES 71% and NEMD 46%, with the overall sensitivity of 56%. However, excluding the nutcracker esophagus or nonspecific disorder which cannot be diagnosed with the radiological methods, the sensitivity increases to 89%. Using videofluoroscopy and 5 time swallows, the average sensitivity was over 90%. In conclusion, the barium method is a simple primary testing method for esophageal motility test. Using not only the image but also the videofluoroscopy with good knowledge of the anatomy and physiology, it is believed that the method will yield the accurate diagnosis.

  • PDF

Sutdy of Appropriate Media Selection and Early Life Cycle of Marine Free-Living Nematodes, Enoplolaimus sp. (Enoplida: Thoracostomopsidae) and Bathylaimus sp. (Enoplida: Tripyloididae) (해양 자유생활형 선충류 Enoplolaimus sp. (Enoplida: Thoracostomopsidae)와 Bathylaimus sp. (Enoplida: Tripyloididae)의 배양용 적합배지 선정 및 초기 생활사 연구)

  • SHIN, AYOUNG;KIM, DONGSUNG;KANG, TEAWOOK;OH, JE HYEOK;LEE, JIMIN;HONG, JAE-SANG
    • The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
    • /
    • v.23 no.3
    • /
    • pp.109-124
    • /
    • 2018
  • In order to find the optimum culture condition for marine free-living nematoda in the laboratory, various agar media were developed and experiments were carried on nematodes. Nematodes have collected from the bottom of the sandy sediments' surface layer (about 5cm sediment) of Taean Mallipo beach's intertidal zone. Especially, with regard to agar medium, Killian medium was transformed slightly, density of agar had made a difference and this agar medium(height 2.0 mm on diameter 60 mm Petri dish) was divided. It was mixed with 5 different species of microorganism as nematodes' live food and added each culture medium. Five individuals of Enoplolaimus sp. on each culture medium were grown in a culture medium which was set on $20^{\circ}C$ and light blocked. Moreover, as a result of the optimum culture condition, 1.0% density of agar showed the highest survival rates (the average time of survival is 246.8 hours). On the other hand, the 0.4% density of agar adding Killian medium, bacto peptone and beef extract showed the lowest survival rates, which indicates the average time of survival is 27.6 hours. About Bathylaimus sp., on Killian medium's 1.0% density of agar(no feeding amount) showed the highest survival rates, which connects that the hatching rate 94.7% after 99.5 hours and it laid 7 eggs averagely on the spawning amount and the hatching rate experiment.

Comparative Analysis of Shallow and Deep Groundwater Pumping Effects on Stream Depletion (천부와 심부지하수 양수에 따른 하천수 감소 영향의 비교분석)

  • Lee, Jeongwoo;Chung, Il-Moon
    • KSCE Journal of Civil and Environmental Engineering Research
    • /
    • v.40 no.4
    • /
    • pp.383-391
    • /
    • 2020
  • In this study, Hunt's analytical solution and Ward & Lough's analytical solution for two-layered leaky aquifer system were used to estimate stream depletions due to shallow and deep groundwater pumping, and their differences were compared. Depending on the combination of the separation distance between the stream and the well, the transmissivity and the storage coefficient of the aquifer, and the leakage coefficient between the upper and lower layers, the stream depletion, which is the amount of stream water reduction compared to the amount of groundwater pumping, for each of 45,000 cases was calculated for both shallow and deep groundwater pumping, and the differences were analyzed quantitatively. When the leakage coefficient was very small, with a value of 10-61/d, the difference in the average five-year stream depletion due to the pumping of shallow and deep groundwater showed a large deviation of up to 0.9 depending on the given hydraulic characteristics; this value exponentially decreased as the stream depletion factor (SDF) increased. This exponential relationship gradually weakened as the leakage coefficient increased due to interaction effects between layers, resulting in a small difference of up to 0.2 when the leakage coefficient reached 10-31/d. Under the condition of greater interlayer hydraulic connectivity, there was little influence of the depth of groundwater pumping on the stream water reduction.

Study on the Characteristics of Far Infrared Ray Drying for Rough Rice(I) (벼의 원적외선 건조특성에 관한 연구(I))

  • 김유호;조영길;조광환;이선호;김영민;한충수;이호필
    • Proceedings of the Korean Society for Agricultural Machinery Conference
    • /
    • 2002.02a
    • /
    • pp.355-361
    • /
    • 2002
  • 본 연구에서는 원적외선.열풍 복합건조특성을 구명하기 위하여 건조용량 150-500kg이고, 승강기, 상.하부스크류, 건조실, 템퍼링실, 송풍기 및 가열장치로 구성된 시뮬레이터를 제작하여 건조특성시험을 실시하였다. 그 결과를 요약하면 다음과 같다. 가. 열풍온도에 따른 곡온변화를 시험한 결과 열풍온도 45$^{\circ}C$일 때 곡온 32-33$^{\circ}C$를 유지하였으며, 48, 51$^{\circ}C$일 때는 곡온이 35$^{\circ}C$가 넘어서는 현상을 나타났다. 건조중 곡온이 35$^{\circ}C$를 넘어서게 되면 동할미 발생량이 많아지고 품질저하가 급격히 일어난다. 나. 템퍼링실의 온도편차가 2,5$^{\circ}C$ 정도로 고른 온도분포를 나타내었고, 버너 입구쪽과 템퍼링실 중앙지점에서 온도가 약간 높게 나타났으며, 원적외선방사체 표면온도분포는 열풍온도가 45$^{\circ}C$일 때 평균 17$0^{\circ}C$를 유지하였고, 48$^{\circ}C$, 51$^{\circ}C$일 때 각각 22$0^{\circ}C$, 23$0^{\circ}C$에서 유지하는 것으로 나타났다. 다. 원적외선방사체 길이방향으로 온도편차는 버너를 기준으로 해서 버너쪽에서 멀수록 온도가 높았고, 중간, 근거리 순으로 나타났다. 버너의 원거리쪽에서 온도가 높게 나타난 것은 원적외선방사체를 통과하는 열풍이 빠져나가도록 되어있는 열풍 유동관이 버너 원거리에 위치하고 있어 버너에 불꽃이 점화되면서 열풍이 방사체 끝쪽으로 일시 머물렀다가 배출되기 때문으로 판단된다. 라. 건조기의 송풍량을 요인으로 하여 건조속도와 건조에너지를 비교한 결과 송풍량이 30cmm일 때가 25cmm에서보다 약 33%의 건조속도가 증가되어 송풍량이 많을수록 건조속도가 빨라졌으나, 건조에너지는 1,391kcal/kg.water로 나타나 약 4.2%정도가 더 소요 되었다. 곡물순환속도를 요인으로 하여 비교 시험한 결과 곡물순환속도가 33kg/min일때가 26kg/min보다 약 25%의 건조속도가 증가되어 곡물의 순환속도가 빠를수록 건조속도가 빨라졌으며, 건조에너지도 1,334kcal/kg.water로 비슷하게 소요되었다. 마. 시험구와 대비구의 건감률은 시험구에서 1.08~1.36w.b./h로 나타나 대비구보다 약 9.9~18.3%가 높게 나타났고, 건조에너지는 10.2~14.6%가 절감되었다. 발아율은 열풍온도가 낮을수록 높게 나타났고 시험구가 대비구보다 발아율이 낮게 나타났으며, 동할률 증가량도 원적외선.열풍 복합건조방법이 높게 나타나 이것은 곡물 표면에 원적외선 방사에의한 복사열이 전달되어 열장해를 받았기 때문으로 판단되며, 금후 더 연구하여 적정 열풍온도 및 방사체 크기를 구명해야 할 것이다.

  • PDF

The Comparison of the Ultra-Violet Radiation of Summer Outdoor Screened by the Landscaping Shade Facilities and Tree (조경용 차양시설과 수목에 의한 하절기 옥외공간의 자외선 차단율 비교)

  • Lee, Chun-Seok;Ryu, Nam-Hyong
    • Journal of the Korean Institute of Landscape Architecture
    • /
    • v.41 no.6
    • /
    • pp.20-28
    • /
    • 2013
  • The purpose of this study was to compare the ultra-violet(UV) radiation under the landscaping shade facilities and tree with natural solar UV of the outdoor space at summer middays. The UVA+B and UVB were recorded every minute from the $20^{th}$ of June to the $26^{th}$ of September 2012 at a height of 1.1m above in the four different shading conditions, with fours same measuring system consisting of two couple of analog UVA+B sensor(220~370nm, Genicom's GUVA-T21GH) and UVB sensor(220~320nm, Genicom's GUVA-T21GH) and data acquisition systems(Comfile Tech.'s Moacon). Four different shading conditions were under an wooden shelter($W4.2m{\times}L4.2m{\times}H2.5m$), a polyester membrane structure ($W4.9m{\times}L4.9m{\times}H2.6m$), a Salix koreensis($H11{\times}B30$), and a brick-paved plot without any shading material. Based on the 648 records of 17 sunny days, the time serial difference of natural solar UVA+B and UVB for midday periods were analysed and compared, and statistical analysis about the difference between the four shading conditions was done based on the 2,052 records of daytime period from 10 A.M. to 4 P.M.. The major findings were as follows; 1. The average UVA+B under the wooden shelter, the membrane and the tree were $39{\mu}W/cm^2$(3.4%), $74{\mu}W/cm^2$(6.4%), $87{\mu}W/cm^2$(7.6%) respectively, while the solar UVA+B was $1.148{\mu}W/cm^2$. Which means those facilities and tree screened at least 93% of solar UV+B. 2. The average UVB under the wooden shelter, the membrane and the tree were $12{\mu}W/cm^2$(5.8%), $26{\mu}W/cm^2$(13%), $17{\mu}W/cm^2$(8.2%) respectively, while the solar UVB was $207{\mu}W/cm^2$. The membrane showed the highest level and the wooden shelter lowest. 3. According to the results of time serial analysis, the difference between the three shaded conditions around noon was very small, but the differences of early morning and late afternoon were apparently big. Which seems caused by the matter of the formal and structural characteristics of the shading facilities and tree, not by the shading materials itself. In summary, the performance of the four landscaping shade facilities and tree were very good at screening the solar UV at outdoor of summer middays, but poor at screening the lateral UV during early morning and late afternoon. Therefore, it can be apparently said that the more delicate design of shading facilities and big tree or forest to block the additional lateral UV, the more effective in conditioning the outdoor space reducing the useless or even harmful radiation for human activities.

Estimation of Hydraulic Parameters from Slug, Single Well Pumping and Step-drawdown Tests (순간수위 변화시험, 단공양수시험 및 단계양수시험을 통한 수리상수 추정연구)

  • Jo, Yun-Ju;Lee, Jin-Yong;Jun, Seong-Chun;Cheon, Jeong-Yong;Kwon, Hyung-Pyo
    • The Journal of Engineering Geology
    • /
    • v.20 no.2
    • /
    • pp.203-212
    • /
    • 2010
  • The aim in this study is used to develop the remediation technologies for contaminated ground water. Slug, single well pumping and step-drawdown tests have been used to obtain hydraulic parameter estimates in the field. Slug tests yield hydraulic conductivity values using the Bouwer and Rice and C-B-P analysis methods. The mean and median hydraulic conductivity values of Bouwer and Rice method are $4.48{\times}10^{-3}$ and $1.16{\times}10^{-3}cm/sec$, respectively. On the other hand, C-B-P method gave mean and median hydraulic conductivity values of $2.37{\times}10^{-3}$ and $7.09{\times}10^{-4}cm/sec$, respectively. These analyses show a trend for the Bouwer and Rice method to yield lower hydraulic conductivity values in low permeability zones of granite in the study area. Sing well pumping test data were calculated through type curve in GW7, GW12 and MW9 wells. It could be interpreted that the differences of hydraulic conductivity and transmissivity values between GW7 and GW12, MW9 is related with fault clays and fractures in the bedrock among the wells. Step-drawdown tests were carried out in the KDPW1 and KDPW2 wells. The hydraulic parameter of KDPW1 and KDPW2 showed very litter difference between the values. The study of hydraulic parameter estimates can be used to purify in contaminated groundwater.

Simulation of Turbidity Flow in the Andon-Imha Linked Reservoir System (안동-임하호 연결 시스템의 탁수유동 모의)

  • Park, Hyung Seok;Chung, Se Woong
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2015.05a
    • /
    • pp.46-46
    • /
    • 2015
  • 강우가 지역별 계절별로 편중되어 있는 우리나라는 수자원의 안정적인 확보와 이용을 위해 다양한 형태의 댐을 건설하여 운영하고 있다. 그러나 대부분의 댐건설을 통해 형성된 저수지들은 탁수 장기화 및 녹조 발생 등의 환경, 생태적인 문제를 겪고 있으며, 그에 따른 사회적 우려로 인해 신규댐 건설을 통한 수자원확보는 더 이상 어려운 실정이다. 이러한 문제에 대응하기 위한 대안으로 기존 댐 저수지들(안동호-임하호)의 구조적 연계운영방안이 진행되고 있다. 본 연구의 목적은 2차원 CE-QUAL-W2모형을 활용하여 안동호와 임하호의 구조적 연결에 따른 탁수의 이동과 각 저수지 내에서의 유동 변화를 해석하는데 있다. 저수지 연계 시나리오는 EL. 138 m 위치에 길이 2 km, 직경 5.5 m 의 콘크리트관(마찰계수 0.05)이 안동호 좌안인 임동면 마리와 임하호 우안 망천리를 연결하는 것으로 가정하였다. 모델의 보정은 실측자료가 풍부한 2006년도 수문사상을 대상으로, 개별 저수지에 대해 수행하였고, 탁수 유동 시나리오 해석은 임하호에 심각한 탁수장기화 문제가 발생했던 2002년을 대상으로 댐 연계 탁수모의를 수행하였다. 안동호와 임하호의 댐 앞에서 모의값과 실측값을 오차를 분석한 결과 탁수예측오차는 AME 0.5~24 mg/L, RMSE 0.7~30.2mg/L의 범위로 비교적 실측값을 잘 반영한 것으로 나타났다. 임하댐의 경우 탁수층의 위치와 두께, 그리고 최고 탁도값을 적절히 재현 하였지만, 안동댐은 최고 탁도값 예측에서 다소 오차가 발생하는 것으로 나타났다. 안동호와 임하호 단독 운영시와 연계 운영시의 탁수변화 파악을 위해 초기 홍수사상이 발생한 8월 이후부터 저수지내의 TSS농도 분포를 비교하였다. 안동호의 경우 댐앞지점의 탁수분포는 수온성층구조에 영향을 받아, 단독 운영시(EL. 130 m)보다 연계운영시(EL. 140 m)에 탁수의 중심이 높은 위치에 형성되었다. 단독 운영시 10월 이후에 전도현상으로 인해 침강되지 않은 잔류 탁수층이 저수지 하부로 확산되었지만, 연계 운영시에는 재부상 되어 상층으로 확산되는 것으로 모의되었다. 또한 연계운영시 유량이동으로 인해 안동호의 탁수 댐앞 도달시간이 짧아지는 것으로 나타났다. 반면 임하호는 연계 운영시 안동댐으로 유출이 생기면서 중층에서 탁수량이 저감되는 것으로 모의되었다. 저수지 내 탁수량 분석을 위해 SS 15 mg/L 이상의 잔류 탁수량을 분석한 결과, 연계운영시 안동호의 평균 잔류탁수량 비율은 11.8% 증가, 임하호의 경우 11.7% 감소하였다. 또한, 탁수의 댐하류 방류일수도 SS 15 mg/L 기준 임하호 9일 저감, 안동호는 70일 증가하여 임하호의 탁수가 안동호의 탁수 장기화에 영향을 주는 것으로 나타났다.

  • PDF