In this study, to investigate an optimal configuration method for the modeling system, we performed an optimization experiment by controlling the types of compilers and libraries, and the number of CPU cores because it was important to provide reliable model data very quickly for the national air quality forecast. We were made up the optimization experiment of twelve according to compilers (PGI and Intel), MPIs (mvapich-2.0, mvapich-2.2, and mpich-3.2) and NetCDF (NetCDF-3.6.3 and NetCDF-4.1.3) and performed wall clock time measurement for the WRF and CMAQ models based on the built computing resources. In the result of the experiment according to the compiler and library type, the performance of the WRF (30 min 30 s) and CMAQ (47 min 22 s) was best when the combination of Intel complier, mavapich-2.0, and NetCDF-3.6.3 was applied. Additionally, in a result of optimization by the number of CPU cores, the WRF model was best performed with 140 cores (five calculation servers), and the CMAQ model with 120 cores ( five calculation servers). While the WRF model demonstrated obvious differences depending on the number of CPU cores rather than the types of compilers and libraries, CMAQ model demonstrated the biggest differences on the combination of compilers and libraries.

Weather is the most influential factor for crop cultivation. Weather information for cultivated areas is necessary for growth and production forecasting of agricultural crops. However, there are limitations in the meteorological observations in cultivated areas because weather equipment is not installed. This study tested methods of predicting the daily mean temperature in onion fields using geostatistical models. Three models were considered: inverse distance weight method, generalized additive model, and Bayesian spatial linear model. Data were collected from the AWS (automatic weather system), ASOS (automated synoptic observing system), and an agricultural weather station between 2013 and 2016. To evaluate the prediction performance, data from AWS and ASOS were used as the modeling data, and data from the agricultural weather station were used as the validation data. It was found that the Bayesian spatial linear regression performed better than other models. Consequently, high-resolution maps of the daily mean temperature of Jeonnam were generated using all observed weather information.

Agricultural meteorological information is an important resource that affects farmersʼ income, food security, and agricultural conditions. Thus, such data are used in various fields that are responsible for planning, enforcing, and evaluating agricultural policies. The meteorological information obtained from automatic weather observation systems operated by rural development agencies contains missing values owing to temporary mechanical or communication deficiencies. It is known that missing values lead to reduction in the reliability and validity of the model. In this study, the hierarchical Bayesian spatio–temporal model suggests replacements for missing values because the meteorological information includes spatio–temporal correlation. The prior distribution is very important in the Bayesian approach. However, we found a problem where the spatial decay parameter was not converged through the trace plot. A suitable spatial decay parameter, estimated on the bias of root–mean–square error (RMSE), which was determined to be the difference between the predicted and observed values. The latitude, longitude, and altitude were considered as covariates. The estimated spatial decay parameters were 0.041 and 0.039, for the spatio-temporal model with latitude and longitude and for latitude, longitude, and altitude, respectively. The posterior distributions were stable after the spatial decay parameter was fixed. root mean square error (RMSE), mean absolute error (MAE), mean absolute percentage error (MAPE), and bias were calculated for model validation. Finally, the missing values were generated using the independent Gaussian process model.

본 연구는 기상학적 및 수문학적 가뭄지수를 이용하여 가뭄사상의 첨두 심도 발생시점과 가뭄발생 기간에 대한 연구를 수행하였다. 연구를 수행하기 위해 사용한 가뭄지수로 기상학적 가뭄지수는 표준강수지수(Standardized Precipitation Index, SPI)를 사용하였으며, 수문학적 가뭄지수로는 하천가뭄지수(Streamflow Drought Index, SDI)와 표준하천유량지수(Standardized Streamflow Index, SSI)를 이용하였다. 연구 대상지역은 농촌과 도시가 공존하는 청미천 유역을 선택하였으며, 평가기간은 1985년 1월부터 2017년 6월까지 32.5년을 평가하였다 하천유량은 SWAT 모형을 이용하여 산정하였다. 수집한 데이터를 이용하여 가뭄지수를 산정한 후에 시계열을 토대로 가뭄의 특성을 분석하였다. 그 결과 수문학적 가뭄은 기상학적 가뭄이 발생한 후에 발생하는 것을 확인할 수 있다. SDI가 SSI보다 SPI와의 첨두 발생시점, 가뭄 시작일의 차이와 평균 가뭄기간이 더 크게 나지만, 가뭄지수의 심도를 비교해보면 일반적으로 SSI가 SDI 보다 심각한 심도를 나타내고 있다. 그러므로 가뭄의 특성을 확인하기 위해서는 기상학적, 수문학적 가뭄지수 등 다양한 가뭄지수를 검토해야 한다.

Background : This study researches on the microclimate, photosynthesis and growth characteristics for the development of shading materials proper for the wide and inclined ginseng cultivation facility which can respond to climate change and save labor. Methods and Results : The wide shading facilities were installed on the area of 1,000 ㎡ in 2014 and 4 facilities were installed on the test ginseng cultivation area. On Mar. 29, 2017, 2 blue shading nets [with the sun blocking rate of 85% (200 g/㎡) and 90% (220 g/㎡)] were installed each for 4 facilities. On June 26, 2017, the aluminum screen and black shading net (with the sun blocking rate of 40% each) were installed during the period of high temperature (30℃ or higher) for each facility. The maximum light intensity under the shading facility was high with blue shading net 85% and PE black shading net 40%, or blue shading net 85% and aluminum screen 40%, higher than other treatments. They were higher by 4.8 - 7.3%, 5.3 - 7.8% each in July and August. Among the coating materials for reducing the high temperature, the aluminum screen coating had less water leakage in early July, late July, mid-August and late August when the precipitation was more than 100 ㎜. The death of the aerial part of ginseng occurred less until October. The growth of the aerial part of 4-year ginseng was better in blue shading net 85% and PE black shading net 40%, blue shading net 85% and aluminum screen 40% or blue shading net 90% and PE black shading net 40% than in blue shading net 90% and aluminum screen 40% The photosynthesis rate was the highest in June with 3.67 μmol CO2/㎡/s under the blue shading net 90% and aluminum screen 40% and with 3.55 μmol CO2/㎡/s under blue shading net 85% and PE black shading net 40% for K-1. As for the land races ginseng, it was the highest with 3.55 μmol CO2/㎡/s under the blue shading net 85% and PE black shading net 40%. For the growth of the underground part of the 4-year ginseng, the blue shading net 85% + PE black shading net 40% or the blue shading net 85% + aluminum screen 40% was the best with respect to the growth of the ground part of the ginseng such as the length of root, the length of main root, diameter of root and weight of root than other treatment. Conclusion : Best coating materials for the wide shading facilities are the blue shading net 85% and aluminum screen 40%.

Background : The continuous cropping of Cnidium officinale is a serious problem for the cultivation practices, which is an unelucidated subject. This study is concerned mainly with rhizosphere microbiome and meteological factors on the cause of physiological damage in the continuous cropping of Cnidium officinale. Methods and Results : Microbial population and community dynamics was evaluated with metagenomic DNA by IonTorrent PGM. Results of HPLC profiling revealed that metabolic components of symbiotic interaction with Cnidium officinale was not detected in cultivated soils. Proteobacteria groups such as nitrogen fixing bacteria, Pseudomonas and Burkholderia of rhizosphere soil in continuous cropped fields mainly decreased compared to the first cropped soil. Principle component analysis of bacterial community showed a significantly differentiated vector between first cropping field and continuous cropped fields. Although growth characteristics including height, leaf length, leaf diameter amd stem diameter etc., was not different with continuous cultivation year until mid-July, physiological damage was dramatically started from late July. Yield of Rhizoma in continuous cropped fields significantly decreased compared to first cropped field. Evapotranspiration of Cnidium officinale with lysimeter for summer season was evaluated. It showed high relationship between solar radiation and evapotranspiration with R2 = 0.7778 and 41% of solar radiation converted into evapotranspiration during 16 days. This result imply that evapotranspiration is mainly controlled by radiation energy in clear days. Water and heat cycle through evapotranspiration is suppose to be one of the important factors related with physiological disorder of Cnidium officinale. Conclusion : This result imply that physiological damage resulted from continuous cropping is involved in decrease of Proteobacteria at rhizosphere soils under stressed conditions.

수문 ․ 기상레이더는 강우량을 바로 추정하지 못하고 여러 단계의 정량적 강우량 추정과정을 거치게 되므로 많은 불확실성 발생요소가 존재한다. 불확실성 관련한 기존 연구들은 정량적 레이더기반 강우량 추정과정에서 보정방법을 이용하여 각 단계별 불확실성을 줄이는 연구들을 수행하였다. 하지만 기존 연구들은 전체 과정에 대한 포괄적인 불확실성을 나타내지 못하고 각 단계별 불확실성의 상대적인 비율도 제시하지 못하는 단점이 있다. 본 연구에서는 정량적 레이더강우량 추정과정의 각 단계별 불확실성을 정량화하고 불확실성 전파를 나타낼 수 있는 적합한 방법을 제시하였다. 첫 번째로 초기와 최종 불확실성, 각 단계별 불확실성의 변동과 상대적인 비율을 나타낼 수 있는 새로운 개념을 제안하였다. 두 번째로 레이더기반 추정과정의 불확실성 정량화와 전파과정을 분석하기 위해 Maximum Entropy Method (MEM)와 Uncertainty Delta Method (UMD)를 적용하였다. 세 번째로 레이더기반 강우량 추정과정의 불확실성 정량화를 위해 2개 품질관리 알고리즘, 2개 강우량 추정방법, 2개 후처리 강우량 보정방법을 2012년 여름철 18개 사례에 대하여 사용하였다. 적용결과, MEM에서 최종 불확실성(후처리 강우량 보정 불확실성: ME = 3.81)이 초기 불확실성(품질관리 불확실성: ME = 4.28)보다 작게 나타났으며, UMD에서도 최종 불확실성(UMD = 4.75)이 초기 불확실성(UMD = 5.33)보다 작게 나타나 불확실성이 감소하는 것으로 나타났다. 하지만 레이더강우량 추정단계의 불확실성은 증가하는 것으로 나타났다. 또한 레이더강우량 추정과정에서 각 단계별로 적합한 방법을 선정하는 것이 각 단계별로 불확실성이 감소시킬 수 있음을 확인하였다. 따라서 본 연구는 새로운 방법이 명확히 불확실성을 정량화할 수 있으며 정확한 정량적 레이더 강우추정에 기여할 것으로 판단한다.

가뭄은 여러 가지 요인에 의해 복합적으로 발생하는 현상으로 자연적 원인과 인위적 원인으로 구분할 수 있다. 우리나라는 기후학적 특성상 여름철에 태평양 고기압의 발달 시기가 빠르거나 평년보다 강하면 장마 기간이 짧아져 자연적인 원인에 의해 가뭄이 발생한다. 가뭄은 발생과정과 피해 영향에 따라 기상학적, 농업적, 수문학적, 사회경제적 가뭄으로 구분할 수 있으며 직 ․ 간접적으로 다른 가뭄에 영향을 미친다. 가뭄의 종류가 기상학적 가뭄에서 농업적 가뭄 혹은 기상학적 가뭄에서 수문학적 가뭄으로 변화되는 현상을 가뭄 전이라 한다. 본 연구에서는 국내의 가뭄 전이 발생여부와 발생 패턴을 검토하기 위해 수문기상 정보를 이용하여 기상학적 가뭄에서 농업적 가뭄으로의 전이 관계를 분석하였다. 가뭄 전이 발생 현황 및 특성 분석을 위해 가뭄 발생의 유형을 5가지로 구분하였으며, 유형에 따라 가뭄 전이가 발생하지 않거나 최대 3개월까지 지체되는 것을 확인하였다. 향후 더 많은 가뭄 지표들과의 분석을 통해 가뭄 전이 관계를 일반화한다면 기상학적 가뭄 발생 시 농업적 가뭄 예측을 위한 인자로 활용할 수 있을 것으로 기대된다.

Both the propagation velocity and the direction of atmospheric waves are important factors for analyzing and forecasting meteo-tsunami. In this study, a total of 14 events of meteo-tsunami over 11 years (2006-2016) are selected through analyzing sea-level data observed from tidal stations along the west coast of the Korean peninsula. The propagation velocity and direction are calculated by tracing the atmospheric disturbance of each meteo-tsunami event predicted by the WRF model. Then, the Froude number is calculated using the propagation velocity of atmospheric waves and oceanic long waves from bathymetry data. To derive the critical condition for the occurrence of meteo-tsunami, supervised learning using a logistic regression algorithm is conducted. It is concluded that the threshold distance of meteo-tsunami occurrence, from a propagation direction, can be calculated by the amplitude of air-pressure tendency and the resonance factor, which are found using the Froude number. According to the critical condition, the distance increases logarithmically with the ratio of the amplitude of air-pressure tendency and the square of the resonance factor, and meteo-tsunami do not occur when the ratio is less than 5.11 hPa/10 min.

Global climate change caused by industrialization has caused abnormal weather conditions such as urban temperatures and tropical nights, urban heat waves, heat waves, and heavy rains. Therefore, the study tried to analyze climate conditions and weather conditions in the streets and analyze climate factors and meteorological factors that lead to inconvenience to citizens. In the case of trees, the overall temperature, surface temperature, solar irradiance, and net radiation were measured low, and the temperature was lower in the Pedestrian road than in roads. The dry bulb temperature, the black bulb temperature, and the wet bulb temperature for the thermal evaluation showed the same tendency. In the case of thermal evaluation, there was a similar tendency to temperature in WBGT, MRT, and UTCI, and varied differences between types. Although the correlation between the meteorological environment and the thermal environment showed a statistically significant significance, the difference between the measured items was not significant. The study found that the trees were generally pleasant to weather and thermal climate in the form of trees, and the differences were mostly documented.

In this study, we investigated the impact of different initial data on atmospheric modeling results using the Weather Research and Forecast (WRF) model. Four WRF simulations were conducted with different initialization in March 2015, which showed the highest monthly mean PM10 concentration in the recent ten years (2006-2015). The results of WRF simulations using NCEP-FNL and ERA-Interim were compared with observed surface temperature and wind speed data, and the difference of grid nudging effect on WRF simulation between the two data were also analyzed. The FNL simulation showed better accuracy in the simulated temperature and wind speed than the Interim simulation, and the difference was clear in the coastal area. The grid nudging effect on the Interim simulation was larger than that of the FNL simulation. Despite of the higher spatial resolution of ERA-Interim data compared to NCEP-FNL data, the Interim simulation showed slightly worse accuracy than those of the FNL simulation. It was due to uncertainties associated with the Sea Surface Temperature (SST) field in the ERA-Interim data. The results from the Interim simulation with different SST data showed significantly improved accuracy than the standard Interim simulation. It means that the SST field in the ERA-Interim data need to be optimized for the better WRF simulation. In conclusion, although the WRF simulation with ERA-Interim data does not show reasonable accuracy compared to those with NCEP-FNL data, it would be able to be Improved by optimizing the SST variable.

In this study, a weighted ensemble method of numerical weather prediction by ensemble models is applied for PyeongChang area. The post-processing method takes into account combination and calibration of forecasts from different numerical models, assigning greater weight to ensemble models that exhibit the better performance. Three different numerical models, including European Center Medium-Range Weather Forecast, Ensemble Prediction System for Global, and Limited Area Ensemble Prediction System, were used to perform the post-processing method. We compared the model outputs from the weighed combination of ensembles with those from the Ensemble Model Output Statistics (EMOS) model for each raw ensemble model. The results showed that the weighted ensemble method can significantly improve the post-processing performance, compared to the raw ensemble method of the numerical models.

This study analyzes relationship between land cover change and local climate data of 46 weather stations in South Korea from the 1980s to 2000s. The used area shows proportional relationships to mean daily min. temperature and mean temperature, and a reverse relationship is found between the used area and relative humidity. However, the forest indicates reverse relationships to mean daily min. temperature and mean temperature. The agricultural land causes to increase in relative humidity and decrease in mean wind speed, while the water increases mean daily min. temperature, daily min. temperature and mean temperature as well as mean wind speed. The urbanization type (used area + barren) shows high correlations with temperature and humidity. The suburbanization type (agricultural land + forest + grass + wetland) has high correlations with temperature and wind speed. High correlations are also found between the waterfront type (wetland + water), and temperature and wind speed. It can be concluded that change in land cover around weather station obviously influences on climate data of the weather station and it is also expected that reliability and homogeneity of climate data from a weather station can be enhanced by this study.

최근 국지성 집중호우 및 급격한 기상변화로 인해 돌발홍수와 같은 기상재해의 발생빈도가 증가함에 따라 고해상도의 기상레이더 강수자료를 사용한 수공학 분야의 연구가 활발하게 진행되고 있다. 기상레이더 강수자료를 활용하는 목적은 기상레이더 강수자료가 제공하는 공간분포를 최대한 활용하는데 있다. 본 연구에서는 고해상도 기상레이더 강수자료의 공간적 특성을 유지하면서 지상 강수자료의 양적특성을 극대화할 수 있는 조건부 합성기법을 보간법에 따라 분석 하였다. 기상레이더 강수자료와 지상 강수자료를 조건부 합성하기 위하여 Kriging, 역거리 가중법 및 Spline 보간법을 적용하였다. 조건부 합성결과는 지상 강수패턴을 현실성 있게 재현하였으며 추가적으로 보간법에 적용되지 않은 강수자료와 모형검증을 수행한 결과 조건부 합성을 통하여 생산된 공간적 강수정보의 수문학적 활용이 가능할 것으로 판단된다.

기후변화에 대한 우려와 함께 증가하고 있는 극한호우의 피해를 줄이기 위해서는 호우사상 발생 이전에 홍수위험을 미리 파악하여 피해를 대비 할 시간을 늘리는 것이 중요하다. 본 연구에서는 기상청 동네예보를 기반으로 하는 간단한 확률적 홍수위험 산정방법을 제시하였다. 예보강수를 조건부로 하는 6시간 강수량의 확률밀도함수를 이용해 다수의 임의 강수량을 생성한 후 추계학적 모형으로 1시간 단위로 분해하여 간단한 강우-유 출모형에 입력하는 방법을 사용하였다. 보청천 유역의 2017년 주요 강우사상에 제안된 방법을 적용한 결과, 7월 4일 최대홍수량이 나타났던 사상에 대해서는 예보강수를 이용한 모의는 홍수위험을 과소평가하였음을 확인하였고 반면 8월 15일 사상에 대한 동네예보는 강수량을 다소 과대추 정하였지만 홍수위험을 충분히 알릴 수 있는 정보로 평가되었다. 본 연구는 확정론적 모형과 확률론적 강수량을 결합하여 기상예보의 불확실성을 고려한 자료기반 홍수위험도 산정방법을 제시한다.

이 논문에서는 동아시아적 자연학의 한 갈래이자 계절의 변화를 천 문기상학적 이론으로 접근시킨 고대 월령론이 성립하고 전개되는 문제 를 살펴보았다. 이때 고대 동아시아적 월령론이 정초되는 사상과 이론 의 기초로서 천문학과 기상학 관점의 논의를 주목하였고, 특히 계절의 변화에 따라 배태되는 자연물의 관찰경험적 지식이 광범위한 물후론으 로 성립하는 측면을 드러내고자 노력하였다. 이를 위해 본고에서는 먼저 그 계절의 변화론을 천문관측이론으로 설명하려는 상서 ｢요전｣의 사중성 문제를 살펴보았고, 이를 확대하 여 매월의 별자리 관찰과 물후로 연결시킨 대대례기 ｢하소정｣의 월 령론을 고찰하였다. 이어서 진제국의 통일이념서 역할을 하였던 여씨춘추 ｢십이기｣에 서 비로소 본격적인 1년 12월령 체계가 수립되었고, 이는 전한시기 회남자 ｢시칙훈｣과 예기 ｢월령편｣으로 계승되었는 바, 이 세 종류 문헌에 수록된 구체적인 천문과 물후월령의 관계를 살펴보았다. 이러는 사이에 월령론은 자연의 계절변화를 설명하는 적합한 이론틀 로 인식되고 확산되었으며, 그 끝에 후한대 성립으로 추정되는 일주서  ｢시훈해｣는 1년 12월 24절기 체제를 더욱 공고하게 만든 1년 72후법 을 처음으로 규범화하고 체계적 설명틀로 수립하였다. 이제 자연의 계 절 변화가 매뉴얼화되어 달력이 없더라도 규범화된 각 절기의 물후 곧 절후법을 통해 쉽게 자연의 변화를 인지하고 생활 속에서 활용할 수 있게 된 것이다. 그 72후의 절후법에는 1년 태양의 변화에 따라 규칙적으로 발생하는 기상상태와 동식물의 반응양태가 다분히 경험학적으로 반영되었으므 로 이 규범체계의 확산력은 매우 광범위하였다 하겠다. 남북조시대에 들어서 6세기 중반에 반포 사용된 위서 ｢정광력｣은 드디어 72절후법을 역법의 주요 요소로 평가하여 수록하였으며, 다만 일주서와는 다른 계통의 72후 체계를 엮어 제시하였다. 이후로 ｢역 지｣ 계통에는 반드시 72후법을 수록하는 것이 항례가 되었다. 수나라 의 <황극력 72후법>이나 당나라 전기의 <인덕력 72후법>은 이를 전승한 주요한 흐름들이며, 당나라 중기 <대연력 72후법>은 기존의 72후 순서나 내용을 다소 개변하였고, 당나라 후기 <선명력 72후법> 은 이 대연력 계통의 72후법을 전승한 것이었다. 우리나라의 경우 고려가 국초부터 당나라 <선명력>을 수용하여 사 용하였으므로, <선명력 72후법>의 이론으로 자연변화를 인식하고 설 명하였을 것이다. 고려말에 도입된 <수시력 72후법>은 기존의 72후 내용에 약간의 변화를 주었으며, 이 체계가 조선전기에도 습용되어 조 선의 학인들이 자연을 인식하는 주요 규범틀로 작용하기에 이르렀다. 이와 같이 72후에 기초한 물후월령은 월령론의 천문기상학적 접근 이자 전근대시대 자연의 변화를 해석하는 주요 통로로 기능하였다고 정리할 수 있다. 이에 따라 일상생활을 영위하는 우리로서 그 물후월령 을 통해 자연과 동화되고 적응하는 삶을 구가하는 것이며, 이는 달리 말하면 우리의 일상생활을 우주론적으로 성화시키는 자연학적 기제가 여기에 놓인 것이라 이를 수가 있을 것이다. 이런 측면에서 본고의 논의 는 월령자연학으로 열리는 방향이라 하겠다. 조선 후기에 널리 성립하는 각종 세시기 문헌의 세시월령이나 각종 시가 작품으로 전하는 월령체 문학 또는 각종 구비전승으로 전개되는 민속월령 문화는 우리 문화사를 구성하는 주요 영역인 바, 그 문화의 기초가 되는 월령자연학의 역사적 흐름과 이에 대한 천문기상학적 원 리 해명을 시도한 것이 본고의 주요 관심사이자 연구 의의라 이를 수 있다.

Freezing rain is a phenomenon when precipitation falls as a liquid rain drop, but freezes when it comes into contact with surfaces or objects. In this study, we investigated the predictability of freezing rain and its characteristics, which are strongly related with the occurrence of black ice using synoptic scale meteorological observation data. Two different cases occurred at 2012 were analyzed and in the presented cases, freezing rain often occurs in the low-level low pressure with the warm front. The warm front due to the lower cyclone make suitable environment in which snow falling from the upper layer can change into supercooled water. The 0℃ temperature line to generate supercooling water is located at an altitude of 850 hPa in the vertical temperature distribution. And the ground temperature remained below zero, as is commonly known as a condition for black ice formation. It is confirmed that the formation rate of freezing rain is higher when the thickness after 1000-850 hPa is 1290-1310 m and the thickness of 850-700 hPa layer is larger than 1540 m in both cases. It can also be used to predict and estimate the generation of freezing rain by detecting and analyzing bright bands in radar observation.

The Weather Research and Forecasting (WRF) model and Vegetation Photosynthesis and Respiration Model (VPRM) were coupled to simulate atmospheric CO2 concentrations. The performance of the WRF-VPRM to simulate regional scale CO2 concentration was estimated over coastal basin areas. Either Hestia 2011(HST) or Vulcan 2002(VUL) anthropogenic CO2 emission data were used in two numerical experiments for the study regions. Simulated meteorological variables were validated with ground and background CO2 measurement data, and the results show that the model captured temporal variations of CO2 concentration on a daily basis. CO2 directional analysis revealed that the dominant CO2 emission sources are located S and SW. The simulated Net Ecosystem Exchange (NEE) agreed relatively well with measured CO2 fluxes at each vegetation class site, showing approximately 40% at max improvement at shrub areas.

이 논문에서는 동아시아적 자연학의 한 갈래이자 계절의 변화를 천 문기상학적 이론으로 접근시킨 고대 월령론이 성립하고 전개되는 문제 를 살펴보았다. 이때 고대 동아시아적 월령론이 정초되는 사상과 이론 의 기초로서 천문학과 기상학 관점의 논의를 주목하였고, 특히 계절의 변화에 따라 배태되는 자연물의 관찰경험적 지식이 광범위한 물후론으 로 성립하는 측면을 드러내고자 노력하였다. 이를 위해 본고에서는 먼저 그 계절의 변화론을 천문관측이론으로 설명하려는 상서 ｢요전｣의 사중성 문제를 살펴보았고, 이를 확대하 여 매월의 별자리 관찰과 물후로 연결시킨 대대례기 ｢하소정｣의 월 령론을 고찰하였다. 이어서 진제국의 통일이념서 역할을 하였던 여씨춘추 ｢십이기｣에 서 비로소 본격적인 1년 12월령 체계가 수립되었고, 이는 전한시기회남자 ｢시칙훈｣과 예기 ｢월령편｣으로 계승되었는 바, 이 세 종류 문헌에 수록된 구체적인 천문과 물후월령의 관계를 살펴보았다. 이러는 사이에 월령론은 자연의 계절변화를 설명하는 적합한 이론틀 로 인식되고 확산되었으며, 그 끝에 후한대 성립으로 추정되는 일주서  ｢시훈해｣는 1년 12월 24절기 체제를 더욱 공고하게 만든 1년 72후법 을 처음으로 규범화하고 체계적 설명틀로 수립하였다. 이제 자연의 계 절 변화가 매뉴얼화되어 달력이 없더라도 규범화된 각 절기의 물후 곧 절후법을 통해 쉽게 자연의 변화를 인지하고 생활 속에서 활용할 수 있게 된 것이다. 그 72후의 절후법에는 1년 태양의 변화에 따라 규칙적으로 발생하는 기상상태와 동식물의 반응양태가 다분히 경험학적으로 반영되었으므 로 이 규범체계의 확산력은 매우 광범위하였다 하겠다. 남북조시대에 들어서 6세기 중반에 반포 사용된 위서 ｢정광력｣은 드디어 72절후법을 역법의 주요 요소로 평가하여 수록하였으며, 다만 일주서와는 다른 계통의 72후 체계를 엮어 제시하였다. 이후로 ｢역 지｣ 계통에는 반드시 72후법을 수록하는 것이 항례가 되었다. 수나라 의 <황극력 72후법>이나 당나라 전기의 <인덕력 72후법>은 이를 전승한 주요한 흐름들이며, 당나라 중기 <대연력 72후법>은 기존의 72후 순서나 내용을 다소 개변하였고, 당나라 후기 <선명력 72후법> 은 이 대연력 계통의 72후법을 전승한 것이었다. 우리나라의 경우 고려가 국초부터 당나라 <선명력>을 수용하여 사 용하였으므로, <선명력 72후법>의 이론으로 자연변화를 인식하고 설 명하였을 것이다. 고려말에 도입된 <수시력 72후법>은 기존의 72후 내용에 약간의 변화를 주었으며, 이 체계가 조선전기에도 습용되어 조 선의 학인들이 자연을 인식하는 주요 규범틀로 작용하기에 이르렀다. 이와 같이 72후에 기초한 물후월령은 월령론의 천문기상학적 접근 이자 전근대시대 자연의 변화를 해석하는 주요 통로로 기능하였다고 정리할 수 있다. 이에 따라 일상생활을 영위하는 우리로서 그 물후월령을 통해 자연과 동화되고 적응하는 삶을 구가하는 것이며, 이는 달리 말하면 우리의 일상생활을 우주론적으로 성화시키는 자연학적 기제가 여기에 놓인 것이라 이를 수가 있을 것이다. 이런 측면에서 본고의 논의 는 월령자연학으로 열리는 방향이라 하겠다. 조선 후기에 널리 성립하는 각종 세시기 문헌의 세시월령이나 각종 시가 작품으로 전하는 월령체 문학 또는 각종 구비전승으로 전개되는 민속월령 문화는 우리 문화사를 구성하는 주요 영역인 바, 그 문화의 기초가 되는 월령자연학의 역사적 흐름과 이에 대한 천문기상학적 원 리 해명을 시도한 것이 본고의 주요 관심사이자 연구 의의라 이를 수 있다.

Land cover around weather station can affect various climatic elements such as temperature, precipitation, wind and moisture. This study tries to reveal patterns and causes in land cover changes around weather stations in South Korea and to classify types of the land cover changes. The agricultural land, forest and used area have occupied most areas around un-relocated 47 weather stations with enough weather observation durations for the past 20-30 years. Area of the agricultural land shows steady decrease with remarkable areal decrease in the wetland, while steady areal increase is found in the used area. Due to urban expansions and new town developments around the weather stations, areas of the urbanization type (used area+barren) and suburbanization type (agricultural land+forest+grass+wetland) show steady increase and decrease, respectively. On the other hand, accurate changes in land cover are not reflected in the water-front type (wetland+water) due to errors and limitations in land cover classification using satellite images. It is expected that the results in this study can contribute to reliability and homogeneity enhancement of weather data.