The effects of climate change on the breeding ecology and phenological asynchrony of the Varied Tit (Sittiparus varius) in Korea

학위논문(박사)--서울대학교 대학원 :농업생명과학대학 산림과학부(산림환경학전공),2020. 2. 이우신.기후변화로 인한 영양단계간 생물계절성 불일치는 생태계 구성원간의 상호작용을 교란시킨다. 특히 산림성 조류와 곤충 먹이원의 생물계절성 변화정도의 차이는 산림성 조류의 번식시기와 곤충 먹이원의 최대 가용시기 간의 불일치로 이어져, 조류의 번식에 부정적인 영향을 미치고 더 나아가 개체군 감소와 산림생태계의 기능 저하를 초래할 수 있다. 그러나 기후변화가 국내 산림성 조류의 생물계절성에 미치는 영향을 파악하기 위한 광역적 규모의 장기 모니터링 체계가 없는 상황이며, 조류의 생물계절성과 번식생태에 대한 기후변화의 영향을 파악하고 예측하는 연구 역시 부족한 실정이다. 본 연구는 박새(Parus major)와 곤줄박이(Sittiparus varius)를 대상으로 광역적 규모의 생물계절성 장기 모니터링에 적용이 가능한 번식시기 자동 모니터링 개발, 장기적 기온 변화와 생물계절성 불일치가 곤줄박이의 번식생태에 미치는 영향 파악, 기후변화에 따른 곤줄박이 산란일 변화 예측을 목적으로 수행되었다. 이를 위해 강원도 원주시 백운산, 전라남도 구례군 지리산, 전라남도 광양시 백운산에서 고도가 다른 12개 조사구에 인공새집 576개를 설치하여 박새와 곤줄박이의 번식지표를 파악하고, 조사구별 기온을 측정하였다. 지리산에서는 2015년부터 2018년까지, 원주시 백운산과 광양시 백운산에서는 2016년부터 2018년까지 조사를 수행하였다. 첫째, 번식시기 자동 모니터링 개발을 위해 2016년부터 2018년까지 12개 조사구에서 소형 온도기록계를 통해 측정한 새집내부온도와 박새와 곤줄박이에 대한 인공새집 번식조사결과를 이용하였으며, 기계학습을 이용하여 번식상황 분류모형을 개발하고 산란일과 부화일을 추정하였다. 둘째, 장기적 기온 변화가 곤줄박이 번식생태에 미치는 영향을 파악하기 위하여, 구례군 지리산의 고도 300m, 900m, 1,300m에 위치한 3개 조사구에서 2015년부터 2018년까지 인공새집 번식을 조사하고, 곤줄박이의 산란일, 한배산란수, 번식실패 둥지 비율과 번식기 전 기온의 연도별 변화를 분석하였다. 셋째, 기후변화에 따른 곤줄박이 산란일 변화를 예측하기 위해, 2016년부터 2018년까지 12개 조사구에서 기록된 3월 기온(3월 10 – 21일)과 산란일(학습자료)을 이용하여 산란일 예측모형을 개발하였다. 또한 2006년부터 2015년까지 지리산 고도 300m, 900m 조사구의 3월 기온과 산란일(검증자료)을 이용하여 산란일 예측모형의 정확도를 검증하였으며, 산란일 예측모형을 통해 RCP (Representative Concentration Pathway) 4.5와 8.5 기후변화 시나리오에 따른 2100년까지의 산란일을 예측하였다. 이를 통해 산란일에 대한 기후변화 영향의 지역간 차이를 파악하고, 기후변화 완화(mitigation) 정도에 따른 산란일 변화 경향의 차이를 분석하였다. 넷째, 고도별 곤줄박이와 주요 먹이원인 나비목 유충 간의 생물계절성 일치도와 생물계절성 불일치가 곤줄박이 번식생태에 미치는 영향 파악을 위해 전라남도 광양 백운산의 4개 고도(200m, 400m, 800m, 1,000m)에 위치한 5개 조사구에서 현장조사를 실시하였다. 고도별 나비목 유충의 발생양상(최대 발생시기, 최대 발생량, 최대 발생기간)과 곤줄박이 유조의 최대 성장시기를 파악하고, 나비목 유충과 곤줄박이간 생물계절성 일치도를 파악하였다. 생물계절성 일치도에 따른 유조의 먹이원 비율과 이소 전 건강도 변화를 파악하기 위해, 곤충 및 기타 먹이원과 유조 혈액의 탄소 및 질소 안정성 동위원소비(δ13C, δ15N)를 분석하여 둥지별 먹이원 비율을 추정하였으며, 13일령 유조의 부척 길이와 체중을 측정하여 이소 전 건강도를 파악하였다. 연구 결과, 첫째, 박새와 곤줄박이에 대한 번식상황 분류모형의 정확도와 AUC (Area Under the Curve)가 모두 0.9이상으로 나타나, 모형의 수행능력이 뛰어난 것으로 확인되었다. 번식상황 추정결과를 이용하여 추정한 산란일과 부화일에 대한 평균 추정 오차(RMSE: Root Mean Square Error)는 박새 2.44 – 4.94일, 곤줄박이 1.31 – 3.03일로 나타나, 박새류에 대한 번식시기 추정의 정확도가 상당히 높고, 종간 번식시기 추정의 정확도에도 차이가 없는 것으로 나타났다. 따라서 박새류로 대표되는 산림성 조류의 생물계절성에 대한 기후변화의 영향을 파악하기 위한 광역적 규모의 장기 모니터링에 온도기록계를 이용한 번식시기 자동 모니터링을 적용한다면, 인적 오류(human error)를 줄이고 조사 방법을 표준화하며, 조사 비용과 시간을 절감할 수 있을 것으로 기대된다. 둘째, 전라남도 구례군 지리산 고도 300m, 900m, 1,300m 조사구의 번식기 전 기온의 연도별 변화를 확인한 결과, 저지대의 기온 변화는 미미한 반면, 고지대에서는 뚜렷한 기온 상승이 확인되었다. 그러나 곤줄박이 산란일은 세 조사구 모두에서 뚜렷한 연도별 변화가 나타나지 않았다. 고지대에서 확인된 기온 상승은 기후변화의 영향이 고지대에서 저지대에 비해 더욱 빠르게 심화될 수 있음을 의미한다. 그러나 고지대에서의 번식기 전 기온과 산란일의 장기적 변화 경향의 차이가 나타났음에도 불구하고, 한배산란수 감소 또는 번식실패 둥지 비율 증가 등 번식지표의 악화는 확인되지 않았다. 향후 곤줄박이 번식에 대한 기후변화의 영향을 파악하기 위해, 번식성공률 및 유조의 이소 전 건강도 등 다른 번식지표에 대한 추가적인 모니터링이 필요하다고 생각된다. 셋째, 곤줄박이 산란일 예측모형의 예측변수로는 3월 일평균기온의 평균과 일최고기온의 평균이 선정되었다. 검증자료와 학습자료에 대한 추정오차는 각각 5.21일, 5.60일로 큰 차이가 없었으며, 검증자료의 산란일과 검증자료의 3월 기온과 산란일 예측모형을 통해 추정된 산란일간 높은 상관관계가 나타났다. 기후변화 시나리오를 이용하여 곤줄박이의 산란일 변화를 예측한 결과, 이산화탄소 배출량이 현재 수준으로 유지될 경우(RCP 8.5), 2090년대의 산란일은 현재보다 10일 빨라질 것으로 나타났다. 반면 기후변화 완화를 통해 2100년 이산화탄소 농도를 540ppm까지 감축할 경우(RCP 4.5), 2090년대 산란일은 현재보다 4일 앞당겨질 것으로 예측되었다. 특히 RCP 8.5 시나리오에서는 고위도 및 고지대의 산란일이 다른 지역에 비해 더욱 빨라지는 것으로 나타나, 국내 고위도 및 고지대 곤줄박이 개체군에 대한 기후변화의 영향이 다른 지역에 비해 더욱 심화될 것으로 예측되었다. 넷째, 고도 상승에 따른 기온 저하로 인해 곤줄박이 유조의 최대 성장시기와 나비목 유충의 최대 발생시기 모두 고도가 높아질수록 지연되었다. 그러나 나비목 유충의 최대 발생시기가 빠를수록 곤줄박이 유조의 최대 성장시기는 유충 발생시기에 비해 상대적으로 늦어져, 기온 상승에 따라 나비목 유충과 곤줄박이간의 생물계절성 불일치가 증가하였다. 나비목 유충의 최대 발생량은 고도가 높아질수록 감소했으나, 최대 발생기간은 고도에 따른 유의한 변화가 나타나지 않았다. 안정성 동위원소비 분석을 통해 곤줄박이 유조의 먹이원 비율을 파악한 결과, 고도에 따라 유조 혈액의 δ13C가 증가하므로 나비목 유충 비율은 감소하고 메뚜기목과 거미목 비율이 높아질 것으로 추정되었다. 이는 나비목 유충의 출현이 늦고 발생량이 적은 고지대에서 번식하는 곤줄박이는 나비목 유충의 일부를 다른 먹이로 대체하여 유조에게 제공한 결과로 생각된다. 유조 혈액의 δ15N는 번식시기가 늦어질수록 감소하였으며, 주요 먹이원의 비율을 재구성한 결과, 번식시기가 늦을수록 메뚜기목의 비율이 유의하게 증가하였다. 이는 최대 발생시기 이후 가용도가 감소하는 나비목 유충 대신 가을까지 가용도가 증가하는 메뚜기목을 더 많이 급이한 결과로 생각된다. 유조의 이소 전 건강도는 생물계절성 불일치가 증가할수록 저하되는 것으로 나타났으며, 따라서 유조의 최대 성장시기에 나비목 유충이 충분하게 제공되는 것은 유조의 성장에 중요한 것으로 판단된다. 그러나 안정성 동위원소비 및 먹이원 비율에 따른 유조의 이소 전 건강도의 변화는 나타나지 않았다. 결론적으로 기온이 상승할수록 나비목 유충과 곤줄박이의 생물계절성 불일치가 증가하며, 이러한 생물계절성 불일치는 대체 먹이의 이용에도 불구하고 유조의 이소 전 건강도를 저하시켜 곤줄박이의 번식에 부정적인 영향을 미칠 수 있는 것으로 확인되었다. 본 연구 결과, 고지대에서 번식하는 곤줄박이 개체군의 산란일과 번식기 전 기온의 장기적 변화 경향의 차이가 나타났으나, 한배산란수 감소 또는 번식실패비율 증가 등의 뚜렷한 번식지표의 악화는 확인되지 않았다. 그러나 고도별 곤줄박이 유조의 최대 성장시기와 나비목 유충 최대 발생시기를 비교한 결과, 기온 상승에 따라 나비목 유충의 최대 발생시기가 앞당겨질수록 나비목 유충과 곤줄박이간의 생물계절성 불일치가 증가하며, 생물계절성 불일치의 증가는 유조의 이소 전 건강도에 부정적인 영향을 미치는 것으로 나타났다. 따라서 최근 10여년 동안 발생한 고지대에서의 급격한 기온 상승에 따라 곤줄박이와 나비목 유충간 생물계절성 불일치의 증가와 그로 인한 유조의 이소 전 건강도 저하가 발생했을 가능성이 있다고 판단된다. 기후변화 시나리오를 통해 예측된 곤줄박이의 산란일은 기후변화 완화정도에 관계없이 21세기 중반까지 지속적으로 앞당겨질 것으로 나타났으며, 또한 기후변화가 현재 속도로 진행될 경우에는 고위도 및 고지대에서의 기온 상승이 빠르게 진행될 것으로 예측되었다. 따라서 앞으로의 기온 상승으로 인해 나비목 유충과 곤줄박이간의 생물계절성 불일치 증가 및 그에 따른 유조의 이소 전 건강도 저하가 점차 심화될 가능성이 있다고 판단되며, 특히 고지대 및 고위도 개체군에 대한 영향이 더욱 클 것으로 예상된다. 곤줄박이와 같이 나비목 유충을 번식기 주요 먹이원으로 이용하는 산림성 조류와 나비목 유충간의 생물계절성 불일치 심화는 어린 새들의 건강도 저하, 나비목 유충 포식을 통한 잎 피해 저감 등 산림성 조류가 제공하고 유지하는 생태적 기능을 약화시켜, 산림생태계 전반에 부정적인 영향을 미칠 수 있다. 따라서 기후변화에 대한 산림성 조류와 산림생태계의 취약성을 평가하고 대응방안 마련을 위해, 나비목 유충과 산림성 조류간 생물계절성 불일치에 대한 지속적인 모니터링과 생물계절성 불일치로 인한 유조의 이소 전 건강도의 저하가 개체의 적응도와 개체군 동태에 미치는 영향 파악에 대한 연구가 추가적으로 수행되어야 할 것이다. 특히 장기자료를 통해 파악한 곤줄박이 산란일의 연도별 변화와 기후변화 시나리오를 통해 예측된 산란일 변화 모두 지역간 차이를 보였으며, 곤줄박이와 나비목 유충 간 생물계절성 일치정도가 지역에 따라 달랐다는 점에서, 국내 산림성 조류에 대한 기후변화의 영향을 파악하기 위해서는 광역적 규모의 장기 모니터링을 통해 생물계절성 불일치의 시공간적 변이를 파악할 필요가 있다고 생각된다.Climate change can alter and disrupt the phenological interaction between organisms from various trophic levels in ecosystems. In temperate forests, breeding populations of forest birds are particularly susceptible to climate change because their reproductive timing and the peak abundance of arthropod prey must be tightly synchronized. Phenological asynchrony between forest birds and arthropod prey due to climate change can therefore lead to a decline in the population of forest birds and negatively affect the functional associations the function of the forest ecosystem. In Korea, there is no broad-scale, long-term monitoring system to determine the effect of climate change on the breeding phenology of forest birds. Few studies have focused on quantitative analysis of the effect of climate change on the breeding ecology and phenological synchrony of forest birds. This study was aimed at developing automated monitoring system of the timing of breeding in the Great (Parus major) and Varied Tits (Sittiparus varius), investigating the effect of climate change on the breeding ecology of Varied Tits and phenological synchrony between Varied Tits and arthropod prey, and predicting shifts in the breeding phenology of Varied Tits. A nestbox survey was conducted to record breeding data of Great and Varied Tits in 576 nestboxes on Mt. Baegun in Wonju (four study plots), Mt. Jiri in Gurye (three study plots), and Mt. Baekwoon in Gwangyang (five study plots). The ambient temperature was also measured at all 12 study plots using a temperature/relative humidity data logger. The field survey was conducted between 2015 and 2018 on Mt. Jiri and between 2016 and 2018 on Mt. Baegun and Mt. Baekwoon. The breeding stages classification models were developed through machine learning using the breeding data and the nestbox temperature recorded by a small temperature data logger (TDL) from 2016 to 2018, thereby creating an automated monitoring system for the breeding timing of the two bird species. The first egg-laying dates and hatching dates were estimated through the classification models and, compared with the dates recorded from nestbox surveys. The annual trends of pre-breeding temperature, egg-laying dates, clutch size, and proportion of failed nests of Varied Tits at three study plots located at 300m, 900m, and 1,300m on Mt. Jiri were analyzed using the long-term nestbox monitoring data (2006 – 2013) and the breeding data from nestbox survey (2015 – 2018). The predictive model for the egg-laying dates was developed by using the egg-laying dates and ambient temperature in March at the 12 study plots from 2016 to 2018 as training data. Also, validation of the model was conducted by comparing observed and predicted egg-laying dates using data from 300m and 900m study plots on Mt. Jiri from 2008 to 2015 as the test data. The egg-laying dates were predicted until 2100 based on the Representative Concentration Pathway (RCP) 4.5 and 8.5 climate change scenarios to investigate how they were affected by climate change and the regional differences observed when climate change mitigation was applied. Further investigation was also conducted on the phenological synchrony between Varied Tits and Lepidoptera larvae along with elevation gradient, and its effects on the breeding ecology of Varied Tits. This was performed using records of the temporal pattern of availability of Lepidoptera larvae (peak date, peak biomass, and peak duration), and the peak date of growth of nestling Varied Tits at five study plots (200m, 400m, 800m, 1,000m) on Mt. Baekwoon (Gwangyang). The proportion of potential prey species in the diets of nestlings was estimated by using carbon and nitrogen stable isotope ratio (δ13C and δ15N) in the whole blood of nestlings and potential prey. The body condition index was then calculated using the weight and tarsus length of nestlings to investigate the effect of phenological asynchrony on the proportion of diets and body condition of nestlings. The accuracy and area under the curve (AUC) of breeding stage classification models were over 0.9, indicating that the performance of the classification model was excellent. In addition, the root mean square error (RMSE) for the estimation of egg-laying dates and hatching dates was 2.44 – 4.94 days for Great Tits and 1.31 – 3.03 days for Varied Tits, which suggested high accuracy in the estimation of breeding timing for both species. The automated monitoring of breeding timing using TDL may therefore be adapted to the broad-scale and long-term monitoring to decrease human error, reduce the cost for field surveys, and standardize the monitoring methods. The rise of pre-breeding temperature was strong in the higher-elevation plots at 900m and 1,300m on Mt. Jiri, but relatively stable in the low-elevation plot on Mt. Jiri over the last decade. Although the mean egg-laying dates shifted more in plots experiencing a greater increase in the pre-breeding temperature, the shift was subtle during the study period. Although there was a difference in the annual trend of pre-breeding temperature and egg-laying dates, there was no evidence of a decrease in clutch size or an increase in the proportion of failed nests. Overall mean of daily mean temperature and daily maximum temperature in March were selected as predictor variables in the predictive model of egg-laying dates of Varied Tits. The RMSE of prediction in the training and test data was 5.60 and 5.21 days, respectively. The observed egg-laying dates of the test data highly correlated with those estimated using the predictive model and temperature in March. Prediction based on the RCP climate change scenario showed that egg-laying dates would be advanced by more than 10 days compared to the present if no mitigation is undertaken (RCP 8.5). On the other hand, under the intermediate mitigation (RCP 4.5), the egg-laying dates were predicted to advance by about 4 days compared to the present. In addition, under the baseline scenario with no mitigation, the egg-laying dates will occur earlier in the higher elevation and latitude regions of South Korea, where the temperature will increase faster. The peak date of growth of the nestling Varied Tits and biomass of Lepidoptera larvae at the five study plots on Mt. Baekwoon in Gwangyang were delayed along the elevation gradient. However, with advancing peak dates of biomass of Lepidoptera larvae, the peak date of growth of the nestlings occurred later, indicating that the phenological asynchrony could increase along with rising temperature. Analysis of stable isotope in the blood of nestlings and potential prey showed that the values of δ13C of blood of nestlings increased with elevation gradient. Estimation of the proportion of potential prey through stable isotope mixing model indicated that such an increase in the values of δ13C were related to decreasing proportions of Lepidoptera larvae and increasing proportion of Araneae and Orthoptera. On the other hand, the values of δ15N in the blood of nestlings decreased over the breeding season. Estimating the proportion of potential prey showed a decrease in the proportion of Orthoptera in the nestlings of late breeding pairs. The body condition of nestlings decreased as phenological asynchrony increased but did not significantly change along with the change in the values of δ13C and δ15N of blood of nestlings and the estimated proportion of potential prey. Consequently, the results indicated that the phenological asynchrony could increase along with a rise in temperature, which could lead to the decrease in body condition of nestlings. This study showed that the advancement of egg-laying dates in Varied Tits was negligible at high elevation plots on Mt. Jiri, where the pre-breeding temperature rose significantly in the past decade, without significant decrease in breeding parameters (clutch size, proportion of failed nest). In comparison with the growth of nestlings along elevation gradient, the peak date of biomass of Lepidoptera larvae and growth of nestlings along with elevation gradients indicated that as peak date of biomass of Lepidoptera larvae advanced with the rise of temperature and, the phenological asynchrony increased, which negatively affected the body condition of nestlings. The inferred increase in phenological asynchrony might have led to the decrease in the body condition of nestlings at higher elevation regions on Mt. Jiri over the past decades. Predicted egg-laying dates of Varied Tits for the coming decades indicated that reproductive phenology would advance under both climate change scenarios (RCP 4.5 and RCP 8.5), and the egg-laying dates would occur earlier in the regions with higher elevation and latitude in South Korea where the temperature would increase faster. Therefore, warming temperature due to climate change may cause an increase in phenological asynchrony, and hence negatively affect the body condition of nestlings of Varied Tit. Such effects may be more severe in the breeding population at higher elevations and latitudes. Phenological asynchrony between forest songbirds and arthropod prey can negatively affect the functions of songbirds such as the top-down control of the density and interactions of leaf-eating Lepidoptera larvae. This can subsequently result in the disruption of structure, function, and services of the forest ecosystem. Therefore, the monitoring of the phenological synchrony between birds and prey should be conducted in order to assess vulnerability and develop a conservation and management strategy for forest birds and ecosystem. Studies on phenological-asynchrony-related decreases in body condition of nestlings and its effects on the population dynamics of birds are also warranted. Additionally, this study showed the regional differences in the annual trend of pre-breeding temperature in the past decade and the predicted shift of egg-laying dates during the 21st century. Broad-scale and long-term monitoring of phenology of forest songbirds should therefore be conducted to evaluate the effect of climate change on the forest songbirds in Korea.제 1장 서론 Ⅰ. 연구 배경 및 필요성 2 Ⅱ. 연구 목적 및 구성 9 Ⅲ. 연구사 1. 자동 기록장치를 이용한 조류 번식 모니터링 13 2. 산림성 조류 번식생태에 대한 기후변화의 영향 16 3. 기후변화에 따른 생물계절성 변화 예측 19 4. 곤충 먹이원과 산림성 조류간의 생물계절성 불일치 21 Ⅳ. 연구 방법 1. 연구 대상종 25 2. 연구 대상지 28 3. 세부 연구별 연구 방법 39 제 2장 온도기록계를 이용한 박새류 번식시기 자동 모니터링 Ⅰ. 서론 45 Ⅱ. 연구 방법 1. 연구 기간 및 대상지 49 2. 인공새집 번식조사 50 3. 인공새집 내부온도 측정 51 4. 박새류 번식시기 자동 모니터링 개발 순서 52 5. 인공새집 번식조사를 통한 번식시기 및 번식상황 파악 53 6. 새집내부온도를 이용한 예측변수 산출 55 7. 번식상황 분류모형 개발 및 검증 61 8. 번식시기 추정 및 정확도 검증 70 Ⅲ. 연구 결과 1. 번식상황 분류모형 개발 및 검증 71 2. 박새류 번식시기 추정 정확도 검증 82 3. 종간 번식시기 추정 정확도 비교 91 Ⅳ. 고찰 97 제 3장 고도별 기온 및 곤줄박이 번식지표의 장기적 변화 Ⅰ. 서론 102 Ⅱ. 연구 방법 1. 연구 기간 및 대상지 106 2. 곤줄박이 번식지표 107 3. 기온자료 111 4. 통계 분석 114 Ⅲ. 연구 결과 1. 고도별 기온 및 산란일 117 2. 번식기 전 기온과 곤줄박이 평균 산란일 120 3. 번식기 전 기온의 연도별 변화 124 4. 곤줄박이 번식지표의 연도별 변화 128 Ⅳ. 고찰 134 제 4 장 기후변화 시나리오에 따른 곤줄박이 산란일 변화 예측 Ⅰ. 서론 139 Ⅱ. 연구 방법 1. 연구 기간 및 대상지 143 2. 곤줄박이 산란일 144 3. 산란일 예측변수 산출 146 4. 곤줄박이 산란일 예측모형 개발 및 검증 148 5. 곤줄박이 산란일 변화 예측 151 Ⅲ. 연구 결과 1. 고도별 3월 기온 및 산란일 153 2. 곤줄박이 산란일 예측모형 개발 157 3. 곤줄박이 산란일 예측모형 검증 161 4. 기후변화에 따른 산란일 변화 예측 163 Ⅳ. 고찰 168 제 5장 곤줄박이와 먹이원의 번식기 생물계절성 불일치 Ⅰ. 서론 174 Ⅱ. 연구 방법 1. 연구 기간 및 대상지 178 2. 나비목 유충의 발생양상 179 3. 곤줄박이 유조의 최대 성장시기 및 생물계절성 일치도 182 4. 곤줄박이 유조의 안정성 동위원소 및 먹이원 비율 분석 183 5. 곤줄박이 유조의 이소 전 건강도 186 6. 통계 분석 187 Ⅲ. 연구 결과 1. 고도별 기온, 곤줄박이 유조 성장시기, 나비목 유충 발생양상 190 2. 고도에 따른 생물계절성 일치도 변화 195 3. 곤줄박이 유조 및 잠재 먹이원의 안정성 동위원소비 198 4. 곤줄박이 유조의 안정성 동위원소비 변화 201 5. 곤줄박이 유조의 먹이원 비율 변화 207 6. 곤줄박이 유조 먹이원 비율과 이소 전 건강도 216 7. 생물계절성 일치도에 따른 곤줄박이 유조의 이소 전 건강도 변화 221 Ⅳ. 고찰 225 제 6장 종합 고찰 및 결론 Ⅰ. 종합 고찰 233 Ⅱ. 결론 238 참 고 문 헌 241 Abstract 275Docto

기혼여성 시간제일자리의 상태의존성(state dependence) 변화

학위논문 (석사)-- 서울대학교 대학원 : 사회과학대학 경제학부, 2018. 2. 이철인.시간제일자리가 노동시장의 유연성과 효율성을 제고하는 역할을 수행 하려면 전일제-시간제-미취업 근로형태간의 원활한 이동이 필요하다. 본 연구에서는 우리나라 20-30대 기혼여성들의 시간제 등 근로형태간 이동성 지표로써 상태의존성(state dependence)을 추정하고, 최근의 변화를 분석하고자 하였다. 동태적다항로짓모형(dynamic multinomial logit model) 추정결과 시간제와 전일제 상호간 이동성이 오히려 낮아지면서 각 근로형태의 상태의존성은 강화되는 모습을 보였다. 이는 전일제일자리의 근로시간 유연성 제고 등으로 시간제로의 전환수요가 줄어드는 가운데 시간제 및 전일제 노동시장 간의 구분은 구조적으로 더 강화된 현상과 부합된다. 최근 시간제일자리는 전일제와 미취업 사이에서 근로형태간 이동을 활발하게 하는 다리(bridge)의 역할이 약화되고 전일제 외에 노동시장 참여옵션을 제공하는, 전일제와는 구분된 추가적 일자리형태로 기능하고 있는 것으로 판단된다.제 1 장 서론 1 제 2 장 선행연구 4 제 3 장 모형설정 7 제 1 절 근로형태 선택모형 7 제 2 절 미관측이질성(Unobserved Heterogeneity)을 고려한 동태적 다항로짓모형 8 제 4 장 분석자료 13 제 5 장 추정결과 18 제 1 절 합동다항로짓모형 추정결과 18 제 2 절 미관측이질성(Unobserved Heterogeneity)을 고려한 동태적 다항로짓모형 추정결과 20 제 3 절 여성 근로형태의 상태의존성 강화 배경과 시사점 26 제 6 장 결론 30 참고문헌 33 Abstract 37Maste

Design and Fabrication of Glass Fresnel Microlens for Optical Scanner

학위논문 (석사)-- 서울대학교 대학원 : 공과대학 전기·정보공학부, 2018. 8. 김용권.본 연구는 유리 재흐름 공정을 통해 실리콘 기판에 미리 결합된 광 스캐너용 유리 프레넬 마이크로렌즈를 제작하는 새로운 방법을 설명한다. 제작된 유리 프레넬 마이크로렌즈가 실리콘 기판에 미리 결합되어있기 때문에 추가적인 수작업 없이 반도체 미세공정을 통해 구동부와 결합이 가능하다. 프레넬 렌즈의 얇고 가벼운 구조적 특징 덕분에 렌즈의 크기가 크며 빠른 구동 속도를 가지는 광 스캐너를 설계할 수 있다. 또한, 기존 구면렌즈와 비교해 초점거리를 짧게 설계할 수 있으며 설계 자유도가 높다. 본 연구에서는 반도체 미세공정을 통한 유리 프레넬 마이크로렌즈를 학계 최초로 제작하였으며 마이크로스캐너와 결합하여 구동에 성공하였다. 제작된 유리 프레넬 마이크로렌즈는 3개의 유리 동심원으로 구성되어 있으며 1.4 ㎜의 지름을 가진다. 평균 표면거칠기는 36 ㎚로 측정되었다. 동일 곡률반경과 지름을 가지는 구면렌즈와 비교해 24 %의 부피와 42.4%의 두께를 가져 경량화에 성공하였다. 유리 동심원의 수의 증가를 통해 부피와 두께를 더욱 크게 감소시킬 수 있다. 마이크로스캐너는 4개의 folded-beam 스프링에 지지되어 있으며 중심에 프레넬 마이크로렌즈가 위치한다. 대기압에서 측정된 광 스캐너의 공진주파수는 1.292 kHz이며 Q값은 107.7로 측정되었다. 17.5 V 직류 바이어스 전압과 17.5 V 진폭의 교류 입력을 인가한 결과 ±50 ㎛ 이상의 큰 횡 방향 변위를 구동했다. 광학 스캔의 시연을 통해 ±1.20°의 레이저 조향 각도가 측정되었다. 제작된 프레넬 마이크로렌즈 스캐너는 지름이 1 ㎜ 이상의 유리 마이크로렌즈가 사용된 광 스캐너 중 가장 빠른 구동 속도를 달성하였으며 대면적 마이크로렌즈 경량화의 새로운 해결 방법이 될 것이다. 또한, 유리 열 재흐름 공정을 통한 프레넬 마이크로렌즈와 같은 비대칭 구조를 제작하는 방법은 보다 복잡한 구조물을 제작하는 데 활용될 수 있다.제 1 장 서론 1 1.1 연구의 배경 1 1.2 연구의 동기 3 1.3 연구의 목적 5 1.4 논문의 구성 5 제 2 장 설계 6 2.1 프레넬 마이크로렌즈 설계 6 2.2 마이크로 스캐너 설계 13 제 3 장 제작 19 3.1 전체 공정 과정 19 3.2 프레넬 마이크로렌즈공정 결과 26 3.3 마이크로스캐너 공정 결과 31 제 4 장 측정 34 4.1 PCB 제작 및 소자와 결합 34 4.2 동적 응답 측정 35 4.3 내열성 평가 40 4.4 광학 측정 42 제 5 장 결론 44 참고문헌 45 Abstract 47Maste

A Study on the Estimation of Non-Ideal Blasthole Pressure of Bulk Blend Explosives

학위논문 (박사)-- 서울대학교 대학원 : 공과대학 에너지시스템공학부, 2018. 8. 전석원.화약의 폭발현상과 암반 파쇄에 대한 과학적인 접근이 완전하지 않고 모델링이 발파분야에 도입되어 활용도를 넓혀가고 있으나, 폭발현상을 이상폭굉모델로부터 산정하기 때문에 밀폐도, 혼합도 등이 독립적인 변수로 작용하는 산업용 화약에는 부적합하다. 또한, 화약의 종류, 밀도 및 성분을 현장에서 조정, 제조하는 모바일 설비를 갖춘 벌크형 혼합폭약에 대해 암반파괴에 활용가능하고 쉽게 사용할 수 있는 비이상 폭발현상의 해석이 요구되고 있다. 화약의 폭발력에 대한 정량화 관점에서 발파공 압력이 발파설계의 입력 자료로 활용되지만, 폭발속도를 근거로 계산식에 의해 구하는 발파공 압력은 그 실제적인 의미와 응용방법에 대해 명확한 이론이 제시되고 있지 않은 실정이다. 폭발속도는 폭굉 진행 구역(DDZ, Detonation Driving Zone)내에서 방출되는 에너지에 의해 결정되고, 산업용화약의 경우에는 폭굉 진행구역 내에서 반응이 완료되는 것이 아니므로, 발파공 압력이 폭발속도만으로 결정되는 것은 아니다. 따라서 이 연구에서는 벌크형 혼합폭약의 비이상 폭발과 이로 인한 폭발속도를 발파공에서 직접 측정하고, 벌크형 혼합폭약의 밀도, 장약의 직경, 발파공 외벽 특성에 따른 구속정도 등이 폭발현상에 영향을 미치는 특성을 규명하였다. 또한, 현장시험 측정 자료를 바탕으로 벌크형 혼합폭약의 비이상 폭발정도를 평가하고, 폭발반응에 의해 발생하는 가스량과 폭발온도를 기반으로 반응에너지의 관점에서 벌크형 혼합폭약의 특성을 규명함으로써 발파공의 압력을 추정할 수 있는 식을 도출하였다. 그리고, 현재까지 측정보다는 계산에 의존하여 왔던 발파공 압력을 직접 측정하였고 벌크형 혼합폭약의 비이상 발파공 압력 추정식으로 대상암반의 특성과 현장조건에 대해 목적하는 발파효과를 얻을 수 있는 폭약의 조건, 혼합비와 밀도를 산정할 수 있는 기법을 도출하였다. 이를 바탕으로 시험발파를 시행하여 발파공의 폭발속도와 압력을 측정하고 추정식의 유효성을 검증하였다. 마지막으로, 본 연구에서 제시한 추정식과 기법은 다양한 국내외 작업환경 및 주변조건에 따른 평가를 통하여 데이터베이스화하고 보완․수정 할 수 있을 것으로 판단되며, 향후, 발파설계단계에서 암반조건, 작업조건에 따른 발파결과를 보다 정확하게 예측하고, 목적하는 발파효과를 향상시킬 수 있는 방안을 수립하는 자료로 활용될 수 있을 것으로 판단된다.목차 List of Figures ............................... ⅳ List of Tables ............................... ⅵ 1. 서론 ............................... 1 1.1 연구배경 ............................... 2 1.2 연구목적 ............................... 3 1.3 연구내용 및 기대성과 ............................... 4 2. 기존문헌 고찰 ............................... 6 2.1 충격파 ............................... 6 2.1.1 충격파 이론 ............................... 6 2.1.2 고체매질에서의 충격파 ............................... 7 2.1.3 충격파속도와 입자속도 ............................... 8 2.2 충격파의 전파 ............................... 9 2.2.1 매질내의 응력파 ............................... 10 2.2.2 응력파의 전파속도 ............................... 11 2.2.3 응력파의 반사 ............................... 11 2.3 폭굉이론 ............................... 12 2.3.1 폭굉파 ............................... 12 2.3.2 폭굉파의 에너지 ............................... 13 2.3.3 폭굉압 ............................... 14 2.3.4 C-J면 압력 ............................... 14 2.3.5 이상 폭굉 ............................... 16 2.3.6 산업용화약의 폭굉특성 ............................... 17 2.4 발파공 압력 ............................... 18 2.4.1 발파공 압력의 의미 ............................... 18 2.4.2 이상 폭발의 발파공 압력 ............................... 18 2.4.3 비이상 폭발의 발파공 압력 ............................... 20 2.4.4 폭굉이론에 기초한 동압력 ............................... 22 2.4.5 동적환경에서의 비이상 폭발 ............................... 23 2.5 폭발생성 가스 에너지 ............................... 26 2.5.1 폭발반응 상태방정식 ............................... 26 2.5.2 가스 에너지 ............................... 27 3. 혼합폭약의 폭발현상 ............................... 31 3.1 벌크형 혼합폭약 ............................... 31 3.1.1 임계밀도 ............................... 33 3.1.2 폭약의 예감화 ............................... 33 3.1.3 혼합비 ............................... 35 3.1.4 혼합폭약의 에너지 ............................... 36 3.2 혼합폭약의 폭발속도 ............................... 39 3.2.1 실험실 시험 ............................... 39 3.2.2 현장시험 ............................... 44 3.3 폭발속도 영향인자 ............................... 48 3.3.1 장약경 ............................... 49 3.3.2 장약경에 대한 결과분석 ............................... 51 3.3.3 장약밀도 ............................... 54 3.3.4 장약밀도에 대한 결과분석 ............................... 55 3.3.5 구속정도 ............................... 58 3.3.6 구속정도에 대한 결과분석 ............................... 59 3.4 폭발속도와 발파공 압력 ............................... 63 3.4.1 위고니오 상태방정식 ............................... 63 3.3.2 Win-CalQ 코드 ............................... 68 3.5 혼합폭약의 발파공 압력 ............................... 72 3.4.1 발파공 압력의 측정 ............................... 72 3.4.2 혼합폭약의 동압력 ............................... 76 3.4.3 혼합폭약의 정압력 ............................... 77 4. 비이상 폭발의 발파공 압력 추정식 ............................... 80 4.1 비이상 폭발 정도 ............................... 80 4.2 비이상 발파공 압력 추정식 ............................... 82 4.2.1 발파공 동압력 ............................... 82 4.2.2 발파공 정압력 ............................... 84 4.2.2 비이상 발파공 압력 추정식 ............................... 85 4.2.4 발파공 압력 비교 ............................... 85 4.3 비이상 발파공 압력 추정식의 검증과 활용 ............................... 88 4.3.1 암반강도 판정기준 ............................... 88 4.3.2 추정식의 검증과 피드백 절차 ............................... 89 4.3.3 추정식의 검증시험 ............................... 90 4.3.4 검증시험 평가 ...............................101 5. 결론 ...............................102 참고문헌 ...............................105 Appendix A. 용어정리 ...............................111 Appendix B. 현장시험 기록표 ...............................113 Abstract ...............................168Docto

안정성 동위원소 분석을 이용한 저어새(Platalea minor) 유조 먹이원의 시공간적 변화

학위논문 (석사)-- 서울대학교 대학원 : 산림과학부(산림환경학전공), 2014. 2. 이우신.우리나라는 국제적 멸종위기종인 저어새(Platalea minor)의 주요 번식지로서 종 보전에 있어 핵심적인 지역임에도 불구하고, 번식 현황과 번식기 채식 습성에 대한 구체적인 정보는 부족한 실정이다. 본 연구는 저어새의 번식현황, 육추기 먹이자원의 지역별 차이 및 연령별 변화 등을 파악하고, 이를 바탕으로 번식기 저어새 서식지의 보전을 위한 관리방안을 제시하고자 하였다. 2013년 3월에서 7월까지 구지도, 수하암, 남동유수지, 칠산도의 저어새 번식 둥지 수를 파악하고, 저어새 유조의 첫째날개깃 끝부분(10일령 성장부분)과 안쪽 깃판(22일령 성장부분)을 수집하였다. 수집한 깃털에서 탄소 안정성 동위원소(δ13C)와 질소 안정성 동위원소(δ15N)를 분석하여 유조 먹이원의 지역 및 연령별 차이를 비교하였으며, 탄소 안정성 동위원소의 이원 혼합 모형을 통해 지역 및 연령별 담수성 먹이 및 해수성 먹이 비율을 추정하였다. 본 연구 결과 구지도 149쌍, 수하암 40쌍, 남동유수지 93쌍, 칠산도 27쌍 등, 총 309쌍의 저어새의 번식이 확인되었다. δ13C분석 결과에 따르면 저어새는 네 번식지에서 공통적으로 담수습지와 갯벌을 모두 채식지로 이용했으며, 특히 담수습지로부터 가장 멀리 떨어져 있는 구지도에서 오히려 담수성 먹이원의 비율이 가장 높게 나타났다. 또한 22일령에 성장한 깃털의 δ13C가 10일령에 비해 높게 나타나, 성조가 유조가 성장함에 따라 담수성 먹이의 급이량을 점차 줄이는 것으로 나타났다. 이러한 먹이원 이용의 변화는 갓 부화한 저어새 유조의 염선 발달이 불완전할 것이라는 가설을 지지하는 것이며, 번식 초기에 담수습지에 대한 의존도가 높을 수 있다는 점을 의미한다. δ15N는 칠산도, 남동유수지 및 수하암, 구지도 순으로 높게 나타났으며, 이는 칠산도의 저어새가 가장 높은 영양단계 또는 질소오염원에 노출된 먹이원을 이용한 반면, 구지도에서는 낮은 영양단계의 먹이원을 이용했거나 질소 오염원의 영향이 적은 채식지를 이용하였기 때문으로 생각된다. 또한 구지도, 수하암, 남동유수지에서는 연령에 따른 변화가 나타나지 않았으나, 칠산도에서는 유조가 성장함에 따라 δ15N가 낮아지는 것으로 나타났다. 이는 육추 초기에 칠산도에서 일시적으로 비료와 같은 질소 오염원의 유입에 따른 결과일 수 있으며, 높은 영양단계에 해당하는 먹이원의 가용성이 일시적으로 증가하였을 가능성도 배제할 수 없다고 생각된다. 앞으로 시기 및 지역별 먹이원 δ15N의 차이에 영향을 미치는 요인을 파악하기 위해서는 먹이원의 안정성 동위원소 분석과 함께 질소 오염원에 대한 분석이 필요할 것으로 생각된다. 연구 대상지인 구지도, 수하암, 남동유수지, 칠산도는 전체 개체군의 22.67% 이상을 차지하여 저어새 보전에 있어 중요한 번식지임이 확인되었다. 따라서 신규 보호구역 설정, 이해당사자간의 협력, 오염 통제, 대중 인식 증진 프로그램 등을 통해 본 연구에서 확인된 주요 번식지에 대한 보호와 인접한 담수습지와 갯벌 채식지의 적응적 관리가 필요할 것으로 생각된다.Ⅰ. 서 론 1 Ⅱ. 연구사 6 1. 저어새에 대한 연구 6 2. 고염분의 먹이가 번식에 미치는 영향에 대한 연구 8 3. 안정성 동위원소 분석을 이용한 먹이원 연구 11 Ⅲ. 연구 대상 지역 13 1. 구지도 13 2. 수하암 14 3. 남동유수지 15 4. 칠산도 16 Ⅳ. 연구 방법 19 1. 번식 현황 19 2. 저어새 유조 깃털 채집 21 3. 깃털 전처리 및 안정성 동위원소 분석 24 4. 혼합 모형을 이용한 먹이원 비율 추정 26 5. 통계분석 28 Ⅴ. 연구 결과 29 1. 번식 현황 29 2. 안정성 동위원소 분석 31 2.1 탄소 안정성 동위원소의 시공간적 변화 33 2.2 혼합 모델을 이용한 먹이원 비율 추정 35 2.3 질소 안정성 동위원소의 시공간적 변화 37 Ⅵ. 고 찰 40 1. 번식 현황 40 2. 저어새 유조의 담수 및 해수성 먹이원 비율 42 3. 질소 안정성 동위원소비의 시공간적 변화 45 Ⅶ. 종합 고찰 및 관리 방안 48 Ⅷ. 결 론 53 Ⅸ. 인용문헌 55Maste

(The) Relations of Community Capacity and Health Behavior Among Seoul Residents

학위논문 (박사)-- 서울대학교 대학원 : 보건학과, 2011.2. 조병희.Docto

Grand Alliance VSB HDTV 수신기를 위한 동기 시스템 설계 및 for Grand Alliance VSB HDTV Receiver

Maste

과학 영재 수업에서 언어적 상호작용을 통하여 본 교사의 발문과 피드백 유형 분석

학위논문(석사) --서울대학교 대학원 :과학교육과(화학전공),2008. 8.Maste

Influences of social capital on collective health status from the perspective of community capacity : a case study of health effects of voluntary organizations networks in Dobong-gu community : using multilevel analysis based on a methodological individualism model

학위논문(석사)--서울대학교 보건대학원 :보건학과 보건정책관리학전공,2007.Maste

GARCH bootstrap method and its properties

Thesis (master`s)--서울대학교 대학원 :경제학부,2003.Maste