The ecosystem provides a diverse array of environmental conditions for organisms, and only those that are capable of successfully adapting to these conditions within their habitats can endure, thrive, and proliferate. Further, the environmental conditions within these habitats can significantly affect the bioavailability of chemicals that are introduced therein, thus resulting in varied adverse impacts on the organisms. The present study aims to evaluate the sensitivity of Yuukianura szeptyckii - a species adapted to riparian - to heavy metals following ISO guideline 11276, with the objective of assessing its potential as an indicator species for ecotoxicological evaluations in riparian habitats. The findings revealed that cadmium and copper both had significant toxic effects depending on their concentrations. For cadmium, the LC50 was 280 mg kg-1, EC50 was 66 mg kg-1, and NOEC and LOEC were 25 and 50 mg kg-1, respectively. For copper, the LC50 was 911 mg kg-1, EC50 was 151 mg kg-1, and LOEC was 50 mg kg-1. Comparative analysis with previous results for the international standard species Folsomia candida and the domestic standard species Allonychiurus kimi indicated that Y. szeptyckii exhibited even greater sensitivity to toxicity values. The adverse effects on survival and reproduction were closely associated with the influx concentration of heavy metals in their bodies. Altogether, the results suggest that Y. szeptyckii is a sensitive species for ecotoxicological assessments in riparian habitats, thus making it suitable as an indicator species, particularly in riparian ecosystems that are characterized by relatively high humidity conditions.

This study conducted an investigation into the effects of fruit type and cultivation practices (organic and conventional) on soil characteristics and soil arthropod communities within apple, blueberry, grape, peach, and pear orchards. The collection of soil arthropod communities was achieved through the utilization of pitfall traps, with concurrent measurements taken for soil moisture content, electrical conductivity, and temperature. The findings of this study unveiled substantial impacts attributed to fruit type and cultivation practices on soil characteristics. Specifically, within organic apple orchards, discernibly higher levels of soil moisture content, electrical conductivity, and temperature were observed when compared to their conventional counterparts. The investigation into soil arthropod communities yielded a total of 1,527 individuals, classified in to five phyla and 15 orders. The range of abundance, species richness, and diversity indices varied across conventional and organic orchards. Cultivation practices were found not to exert a significant influence on soil arthropod community characteristics. However, Non-metric Multidimensional Scaling (NMDS) analysis indicated a significant differentiation in soil arthropod community structure based on cultivation practices. This study underscores the importance of considering vegetation structure and environmental characteristics that may influence soil arthropod communities comprehensively when assessing the impact of cultivation practices on soil arthropods. Furthermore, it emphasizes the need to account for both the characteristics and structure of soil arthropod communities in understanding the implications of cultivation practices on these organisms.

The key to invasive pest management lies in preemptive action. However, most current research using species distribution models is conducted after an invasion has occurred. This study modeled the potential distribution of the globally notorious sweet potato pest, the sweet potato weevil (Cylas formicarius), that has not yet invaded Korea using MaxEnt. Using global occurrence data, bioclimatic variables, and topsoil characteristics, MaxEnt showed high explanatory power as both the training and test areas under the curve exceeded 0.9. Among the environmental variables used in this study, minimum temperature in the coldest month (BIO06), precipitation in the driest month (BIO14), mean diurnal range (BIO02), and bulk density (BDOD) were identified as key variables. The predicted global distribution showed high values in most countries where the species is currently present, with a significant potential invasion risk in most South American countries where C. formicarius is not yet present. In Korea, Jeju Island and the southwestern coasts of Jeollanam-do showed very high probabilities. The impact of climate change under shared socioeconomic pathway (SSP) scenarios indicated an expansion along coasts as climate change progresses. By applying the 10th percentile minimum training presence rule, the potential area of occurrence was estimated at 1,439 km2 under current climate conditions and could expand up to 9,485 km2 under the SSP585 scenario. However, the model predicted that an inland invasion would not be serious. The results of this study suggest a need to focus on the risk of invasion in islands and coastal areas.

Many changes in the scale and structure of the Korean rice cropping system have been made over the past few decades. Still, insufficient research has been conducted on the sustainability of this system. This study analyzed changes in the Korean rice cropping system’s sustainability from a system ecology perspective using an emergy approach. For this purpose, an emergy table was created for the Korean rice cropping system in 2011, 2016, and 202, and an emergy-based indicator analysis was performed. The emergy analysis showed that the total emergy input to the rice cropping system decreased from 10,744E+18 sej year-1 to 8,342E+18 sej year-1 due to decreases in paddy field areas from 2011 to 2021, and the proportion of renewable resources decreased by 1.4%. The emergy input per area (ha) was found to have decreased from 13.13E+15 sej ha-1 year-1 in 2011 to 11.89E+15 sej ha-1 year-1 in 2021, and the leading cause was a decrease in nitrogen fertilizer usage and working hours. The amount of emergy used to grow 1 g of rice stayed the same between 2016 and 2021 (specific emergy: 13.3E+09 sej g-1), but the sustainability of the rice cropping system (emergy sustainability index, ESI) continued to decrease (2011: 0.107, 2016: 0.088, and 2021: 0.086). This study provides quantitative information on the emergy input structure and characteristics of Korean rice cropping systems. The results of this study can be used as a valuable reference in establishing measures to improve the ecological sustainability of the Korean rice cropping system.

The process of biological invasion is led by the dynamics of a population as a demographic and evolutionary unit. Spatial structure can affect the population dynamics, and it is worth being considered in research on biological invasion which is always accompanied by dispersal. Metapopulation theory is a representative approach to spatially structured populations, which is chiefly applied in the field of ecology and evolutionary biology despite the controversy about its definition. In this study, metapopulation was considered as a spatially structured population that includes at least one subpopulation with significant extinction probability. The early phase of the invasion is suitable to be analyzed in aspects of the metapopulation concept because the introduced population usually has a high extinction probability, and their ecological·genetic traits determining the invasiveness can be affected by the metapopulation structure. Although it is important in the explanation of the prediction of the invasion probability, the metapopulation concept is rarely used in ecological research about biological invasion in Korea. It is expected that applying the metapopulation theory can supply a more detailed investigation of the invasion process at the population level, which is relatively inadequate in Korea. In this study, a framework dividing the invasive metapopulation into long- and middle-distance scales by the relative distance of movement to the natural dispersal range of species is proposed to easily analyze the effect of a metapopulation in real cases. Increased understanding of the mechanisms underlying invasions and improved prediction of future invasion risk are expected with the metapopulation concept and this framework.

공간 샘플링은 공간모델링 연구에 활용되어 샘플링 비용을 줄이면서 모델링의 효율성을 높이는 역할을 한다. 농업분야에서는 기후변화 영향을 예측하고 평가하기 위한 고해상도 공간자료 기반 모델링에 대한 연구 수요가 빠르게 증가하고 있으며, 이에 따라 공간 샘플링의 필요성과 중요성이 증가하고 있다. 본 연구는 국내 농지 공간샘플링 연구를 통해 농업분야 기후변화연구의 공간자료 활용의 효율성을 제고하고자 하였다. 본 연구는 층화랜덤샘플링 을 기반으로 하였으며, 1 km 해상도의 농지 공간격자자료 모집단 (11,386개 격자)에 대해서 RCP 시나리오별 (RCP 4.5/8.5) 연대별 (2030/2050/2080년대) 공간샘플링을 설 계하였다. 국내 농지는 기상 및 토양 특성에 따라 계층화 되었으며, 샘플링 효율 극대화를 위해 최적 층화 및 샘플 배정 최적화를 수행하였다. 최적화는 작물수량, 온실가스 배출량, 해충 분포 확률을 포함하는 16개 목표 변수에 대해 주어진 정밀도 제한 내에서 샘플 수를 최소화하는 방향으로 진행되었다. 샘플링의 정밀도와 정확도 평가는 각각 변동계수 (CV)와 상대적 편향을 기반으로 하였다. 국내 농지 공간격자 모집단 계층화 및 샘플 배정 및 샘플 수 최적화 결과, 전체 농지는 5~21개 계층, 46~69개 샘플 수 수준에서 최적화되었다. 본 연구결과물들은 국내 농업시스템 대표 공간격자로써 널리 활용될 수 있을 것으로 기대된다. 또한, 기후변화 영향예측 공간모델링 연구들에 활용되어 샘플링 비용 및 계산 시간을 줄이면서도 모델의 효율성을 높이는 데에 기여할 수 있다.

생물다양성협약 (1992년)과 나고야의정서 (2010년)의 체결 이후 우리나라는 한반도의 생물다양성 보전과 생물자원 확보를 위한 자생생물의 조사·발굴 연구에 박차를 가하였다. 이러한 계기로 2007년에 설립된 환경부 소속 국립 생물자원관의 주도로 “한반도 자생생물 조사·발굴 연구사업”이 진행되었다. 본 사업은 2006년 이후 현재까지 15년 동안 5단계 (1단계 2006~2008년, 2단계 2009~2011년, 3단계 2012~2014년, 4단계 2015~2017년, 5단계 2018~2020 년)로 나누어 진행되었다. 연구의 결과, 본 사업의 이전에 29,916종 (2006년)이던 한반도 자생생물이 본 사업의 각 단계가 마무리되는 시점에서 누계로 집계하여 볼 때, 1단계 33,253종 (2008년), 2단계 38,011종 (2011년), 3단계 42,756종 (2014년), 4단계 49,027종 (2017년), 그리고 5단계 54,428종 (2020년)으로 급속히 증가하여 본 사업 기간 동안 한반도 자생생물 기록종이 약 1.8배 증가하였다. 이 통계자료는 이 기간 동안 연평균 2,320종의 한반도 미기록종이 새로이 기록된 것을 보여준다. 또한 전체 발굴종 중에서 총 5,242 종의 신종을 기록하는 학술적 큰 성과를 거두었다. 분류군 별로는 총 연구 기간 동안 곤충 4,440종 (신종 988종 포함), 무척추동물 (곤충 제외) 4,333종 (신종 1,492종 포함), 척추 동물 (어류) 98종 (신종 9종 포함), 식물 (관속식물과 선태식물) 309종 (관속식물 176종, 선태식물 133종, 신종 39종 포함), 조류 (algae) 1,916종 (신종 178종 포함), 균류와 지의류 1,716종 (신종 309종 포함), 그리고 원핵생물 4,812종 (신종 2,226종 포함)이 한반도에서 새로이 기록되었다. 생물표본은 각 단계별로 집계하여 볼 때 1단계 247,226점 (2008년), 2단계 207,827점 (2011년), 3단계 287,133점 (2014년), 4단계 244,920점 (2017년), 그리고 5단계 144,333점 (2020년)이 수집되어 연평균 75,429점, 총 1,131,439점의 생물표본이 채집되었다. 그중에서 곤충 281,054점, 곤충 이 외의 무척추동물 194,667점, 척추동물 (어류) 40,100점, 식물 378,251점, 조류 (algae) 140,490점, 균류 61,695점, 그리고 원핵생물 35,182점이 채집되었다. 본 사업에 참여한 각 단계별 연구원/보조연구원 (주로 대학원생)의 수는 1단계 597/268명, 2단계 522/191명, 3단계 939/292명, 4단계 575/852명, 그리고 5단계 601/1,097 명으로 전체년도의 참여연구자는 연평균 395명, 총 연인원 약 5,000명이 참여하여 전국의 거의 모든 분류학자와 분류학 전공의 대학원생이 참여하였다. 본 사업 기간 동안 전문 학술지 논문 3,488편 (국내학술지 논문 2,320편, SCI급 국제 학술지 논문 1,168편 포함)이 출판되었다. 본 사업 기간 중 자생생물 조사·발굴 사업 및 생물표본 확보 사업에 투입된 예산은 총 833억원 (연평균 55억원)이다. 본 사업은 국가 주도의 대형 연구 프로젝트로서 전국의 거의 모든 분류학자가 참여하고 대규모 예산이 투입되어 단기간에 이루어 낸 한국식 압축성장의 한 성공 사례로 볼 수 있다. 본 사업의 종발굴 성과는 최근의 생물분류 체 계로 분류되어 국가생물종목록으로 만들어졌으며, 전문가와 학생 및 일반 시민에게 제공되고 있다 (https://species. nibr.go.kr/index.do). 본 사업에서 파생된 기재문, DNA 염기 서열, 서식처, 분포, 생태, 이미지, 멀티미디어 등 각 종의 정보는 디지털화되어 생물의 계통, 진화 연구 등 학문적 발전에 기여하였고, 기후변화에 따른 지표종의 변화 같은 생물 분포 모니터링 사업과 바이오산업의 생물소재를 탐색하는 기반이 되었다. 본 사업을 통하여 젊은 분류인력 (주로 대학 원생)의 양성을 지원할 수 있었던 것은 본 사업이 가져온 가장 의미 있는 성과라고 할 수 있다. 과거 15년간 숨 가쁘게 달려온 본 사업은 아직 진행 중이다. 그동안 발굴된 종 들에 대한 이명 (synonym)과 오동정 등을 바로잡아 학문적인 완성도를 높이고, 한반도에 존재하리라 예상되는 약 10만 종의 자생생물 중에서 남겨진 5만 종에 대한 조사·발굴 연구가 지속되어야 한다.

본 연구의 목적은 질산, 황산, 암모니아수, 과산화수소에 대한 생태독성평가를 통해서 사고대비물질들에 대한 기초 독성 데이터베이스를 구축하여, 향후 화학사고 발생시 환경 피해에 관한 의사결정에 과학적 근거를 제공하는 데 있다. 이를 위해 본 연구에서는 사고대비물질 중 토양의 물리·화학적 성질을 변화시킬 수 있는 질산, 황산, 암모니아수, 과산화수소를 대상으로 국내 토착 절지동물인 김어리 톡토기 (Paronychiurus kimi)를 이용한 생태독성평가를 수행하였다. 7일간의 급성독성평가와 28일간의 만성독성평가를 수행하였으며, 시험물질 농도에 따른 토양의 pH 변화를 관찰하였다. 토양의 pH는 질산, 황산, 과산화수소, 암모니아수의 농도가 10,000 mg kg-1 soil dry wt.일 때, 각각 2.86, 2.72, 7.18, 9.69이었다. 질산, 황산, 과산화수소, 암모니아수에 대한 만성독성평가 결과, LC50 값은 각각 2,703, 5,414, 3,158, 859 mg kg-1 soil dry wt.이었으며, P. kimi의 산란 수에 대한 EC50 값은 각각 587, 2,148, 1,300, 216 mg kg-1 soil dry wt.이었다. 비록 본 연구에서는 사고대비물질들의 유입에 따른 토양 pH의 변화만이 조사되었지만, 본 연구의 결과는 P. kimi가 사고대비물질에 의해 변화된 토양의 pH뿐만 아니라 사고대비물질의 유입에 의해 감소된 유기물 함량과 생성된 반응 산물에 의해서도 사망률과 산란수에 영향을 받을 수 있음을 의미한다. 대부분의 사고대비물질들이 토양의 특성을 변화시킬 수 있다는 점을 감안할 때, 토양의 특성 변화와 이에 따른 생물 영향을 고려한 화학사고 후 평가 및 복원 방법이 필요하다.

It is important to consider interaction between species in understanding structure and function of the biological community. Current ecological issues such as climate change and habitat loss emphasize the significance of the concept of species interaction in that varying species’ interaction across environmental gradients may lead to altered ecological function and services. However, most community studies have focused on species diversity through analysis of quantitative indices based on species composition and abundance data without considering species interactions in the community. ‘Ecological network analysis’ based on network theory enables exploration of structural and functional properties of ecosystems composed of various species and their interactions. In this paper, network analysis of Cheonggye stream as a case study was presented to promote uses of network analysis on ecological studies in Korea. Cheonggye stream has a simple biological structure with link density of 1.48, connectance 0.07, generality 4.43, and vulnerability 1.94. The ecological network analysis can be used to provide ecological interpretations of domestic long-term monitoring data and can contribute to conserving and managing species diversity in ecosystems.

Acute and chronic toxicities of methyl ethyl ketone and methanol were investigated on Paronychiurus kimi (Collembola), for evaluating the potential effects of accidental exposures of these chemical substances on the terrestrial environments. This study was undertaken to establish a toxicity database for these chemical substances, which was required for the preparation of the response compensation and liability act for agricultural production and environmental damage. The 7-d acute toxicity and 28-d chronic toxicity were conducted using the OECD artificial soil spiked with varying, serially diluted concentrations of methyl ethyl ketone and methanol. Mortality was recorded after 7-d and 28-d of exposures, and the number of juveniles were determined after 28-d of exposure in the chronic toxicity test. In both assessments, methanol was more toxic than methyl ethyl ketone in terms of mortality (LC50) and reproduction (EC50). The 7-d LC50 of methyl ethyl ketone and methanol were 762 and 2378 mg kg-1 soil dry wt., respectively, and the 28-d LC50s were 6063 and 1857 mg kg-1 soil dry wt., respectively. The 28-d EC50 of methyl ethyl ketone and methanol were 265 and 602 mg kg-1 soil dry wt., respectively. Comparison of results obtained in this study with literature data revealed that P. kimi was more sensitive to methanol than other soil invertebrates. However, given the high volatility of the chemicals tested in this study, further studies are necessary to improve the current test guideline, or to develop new test guidelines for an accurate assessment of chemicals that require toxicity databases for chemical accidents.

Greenhouse gas emission from agricultural land is recognized as an important factor influencing climatic change. In this study, the national CO2 emission was estimated for paddy soils, using soil GHG emission model (DNDC) with 1 km2 scale. To evaluate the applicability of the model in Korea, verification was carried out based on field measurement data using a closed chamber. The total national CO2 emission in 2015 was estimated at 5,314 kt CO2-eq, with the emission per unit area ranging from 2.2~10.0 t CO2-eq ha-1. Geographically, the emission of Jeju province was particularly high, and the emission from the southern region was generally high. The result of the model verification analysis with the field data collected in this study (n=16) indicates that the relation between the field measurement and the model prediction was statistically similar (RMSE=22.2, ME=0.28, and r2=0.53). More field measurements under various climate conditions, and subsequent model verification with extended data sets, are further required.

Much of the data used in the analysis of environmental ecological data is being obtained over time. If the number of time points is small, the data will not be given enough information, so repeated measurements or multiple survey points data should be used to perform a comprehensive analysis. The method used for that case is longitudinal data analysis or mixed model analysis. However, if the amount of information is sufficient due to the large number of time points, repetitive data are not needed and these data are analyzed using time series analysis technique. In particular, with a large number of data points in the current situation, when we want to predict how each variable affects each other, or what trends will be expected in the future, we should analyze the data using time series analysis techniques. In this study, we introduce univariate time series analysis, intervention time series model, transfer function model, and multivariate time series model and review research papers studied in Korea. We also introduce an error correction model, which can be used to analyze environmental ecological data.

We developed a new ecotoxicological test method using native test species, eggs of Ephemera orientalis, and five kinds of industrial wastewater were tested to validate the test method. The water samples were collected in Jun 2006 from the following industries: pesticide, metal plating, PCB, leather1, and leather2. Wastewater and effluent were diluted by distilled water, respectively, to prepare various concentrations, 100, 50, 25, 12.5, 6.3, 3.1, and 0%. For the egg bioassay, 20 freshly laid eggs (<24 h old) were exposed to test solutions in a Petri dish (52×12 mm) at 20℃ with photoperiod of 16 h light and 8 h dark for 14 days. The median egg hatching concentrations (EHC50) were estimated using Probit analysis. All EHC50s of wastewater were less than 3.1%, which meant very high ecotoxicity except for the wastewater of PCB industry having 6.1% of EHC50. Among the effluents, the least toxic effluent was from pesticide industry having 58% of EHC50, while the effluent of leather2 was the most toxic having 7.3% of EHC50.

기후변화로 인한 온도 증가는 농업생태계에서 작물에 대한 해충-천적 시스템의 상호작용을 변화시켜 해충의 발생과 천적의 생물 방제 효과에 영향을 미친다. 고추의 대표적인 해충 복숭아혹진딧물과 목화진딧물, 천적인 칠성무당벌레를 대상으로 각각 기존의 Rosenzweig-Macathur predator-prey모델을 이용하여 244일 동안 평균 온도 상승에 대한 개체군의 밀도 변동을 모의하였다. 기후변화의 온도 영향을 알아보기 위해 모델을 구성하는 각각의 생물 파라미터에 대해 온도 의존 함수가 추가되었고, 수정된 모델의 결과를 바탕으로 dynamic index를 이용하여 상호작용 강도를 산출하였다. 시뮬레이션 결과, 두 시스템 모두 평균온도가 증가함에 따라(+1°C, +3°C) 해충과 천적의 밀도 변동 주기와 해충의 최대 발생 밀도는 감소하고, dynamic index는 증가하였다. 또한, 온도가 5°C 증가할 경우, 목화진딧물-칠성무당벌레 시스템에서는 강한 상호작용 에 의한 포식자 및 피식자의 과도한 밀도 감소의 영향으로 dynamic index가 급격히 상승하는 구간을 보여주었다. 기후변화에 의해 평균 온도가 증가함에 따라 해충의 발생 빈도 증가 및 일부 발생 밀도의 증가가 나타날 수 있지만, 종간 상호작용 증가에 의해 그 영향이 상쇄됨에 따라 생물적 방제 효과는 증대될 것으로 예측된다.

기후변화에 따른 해충개체군 증감모형은 해충방제를 위 한 초발생예찰과 연속적 해충변동 양상의 파악에 매우 중 요하다. 이러한 예측은 농약사용의 효율성을 높이고, 환경 에 적은 영향을 줄 수 있으므로, 현대 해충방제전략의 화두 로 볼 수 있다. 본 연구는 온도변화에 따른 해충의 농약효과 에 따른 사충률의 변화를 개체군 모형과 결합시켜 모의했다. 감수성 점박이응애를 강낭콩을 기주로 20, 25, 30, 35℃에서 Acrinathrin-Spiromesifen 혼합제와 Azocyclotin 유기주석계 농약에 노출시켰다. 생물검정 결과 점박이응애의 사충률은 온도와 농약의 종류에 따라 유의한 차이가 발생했다. 점박이 응애의 개체군 밀도변동 모의는 DYMEX를 이용했으며, 모 의결과 농약의 종류별로 기후변화에 따른 초기방제 시기와 방제횟수에 차이가 나타날 것으로 예측됐다. 본 연구결과는 미래의 기후변화에 대응한 해충방제 전략과 농약 선발에 있 어 중요한 시사점을 제공할 것으로 사료된다.

기후변화로 인한 환경 요인 변화는 생물 간 상호작용에 영향을 주지만 기후변화 관련 곤충 영향 연구 대부분은 이를 배제한 채 비생물 환경인자와의 관계 분석에 초점을 맞추고 있다. 개체군 동태 또는 종 분포 모델을 중심으로 한 단일 종에 대한 연구가 주를 이루며, 영양단계 간 생물요인의 영향을 고려한 생태계시스템 측면의 연구는 부족하다. 농생태계의 경우 지역적 규모로 기후요인 뿐만아니라 생물요인의 영향을 특히 많이 받는다. 따라서, 기후변화의 영향요인으로 종간 최적 생존 온도 범위 및 온도 반응 차이, 생물 파라미터의 온도에 대한 비선형성으로 인한 생물 시스템 간 연결성 변화 및 이에 따른 멸종이 고려되어야 한다. 이를 위해 작물, 해충, 천적상을 포함하는 tritrophic 수준에서 시스템적 접근을 통해 기후변화에 따른 개체군 및 상호작용 dynamic을 구현할 필요가 있다. 생물요인을 포함하는 모델화 작업에 대한 연구 방향으로는 종의 공간적 분포를 고려하기 위해 작물재배적지 모형을 기반으로 재배적합도에 따른 작물의 미래 분포가능 지역을 결정하며, 작물에 대한 해충과 천적상의 공간적 synchrony를 전제로 농생태계 종 구성 시나리오를 작성한다. 이에 tritrophic dynamic system 모형을 적용하여 기후변화에 따른 생물 간 상호작용 강도 변화에 따른 종 구성 변화 및 시스템 변동을 예측하고, 농생태계 작물별 key trophic chain 및 해충압 변화를 시뮬레이션하여 기후변화의 통합 영향을 분석한다.

오늘날 통계적 사고와 통계적 기법의 활용은 학문분야에 있어 보편화되고 있는 추세이다. 또한, 국내 학술지가 국제 학술지로 격상되고 있는 현재, 논문에서의 통 계 기법 사용의 적정성을 점검해 볼 필요가 있을 것이다. 본 논문에서는 응용곤충 학 관련 학술논문집간 통계적 기법의 사용빈도를 비교해봄으로써 국내논문의 통 계 기법 사용의 현재 위치에 대해 고찰하고, 국내논문이 나아갈 방향에 대해 모색 해보았다. 본 연구에는 통계 빈도에 대한 조사를 비교하는 데에는 2005년에서 2009년까지 5년간 국내·외 학술논문에 기재된 논문 대상으로 Fisher's exact test를 사용하였다. 우선 국내와 일본 학술논문집간 비교 결과, 2007년까지는 통계적 처 리기법의 사용빈도에서 유의한 차이를 보였으나 그 이후에는 유의한 차이가 나타 나지 않았다. 국내와 네덜란드의 학술논문집간 비교에서는 2009년을 제외하고는 통계적 처리기법의 사용빈도에서 모두 차이가 있는 것으로 나왔다. 참고적으로 일 본과 네덜란드의 비교에서는 통계적 처리기법의 사용빈도에서 차이가 없었다. 이 러한 결과를 토대로, 국내 학술논문집에서의 통계적 기법의 사용빈도가 점차 국외 학술논문집들과 비슷해지는 과정에 있다는 것을 알 수 있었다. 따라서 논문에서의 통계적 기법을 사용하는 것에 더욱 신중하고 노력을 기울일 필요가 있겠다.

우리나라의 도시하천에 흔히 서식하는 조각깔따구(Glyptotendipes tokunagai Sasa)를 수서생태계의 새로운 수질지표종 또는 독성실험종으로 개발하기 위해 실온에서 계대사육하였다. 사육을 위해 우화한 성충을 교미･산란실(￠120×H210 mm)로 옮겨주어 교미와 산란을 유도하였으며, 산란된 알집은 항온기(온도 25±0.5℃; 광주기 16L:8D)로 옮겨 유충기간 동안 사육하였다. 우화할 시기가 되면 항온기에서 우화실(W550×L420×H415 mm)에 있는 수조(W300×L220×H220 mm)로 옮겨주어 우화할 수 있는 공간을 넓혀 주었다. 유충의 먹이는 Tetramin® (TetraWerke, Melle, Germany)을 0.2 mm이하로 갈아서 제공하였고, 배양액으로는 하루이상 폭기시킨 1차 증류수를 사용하였다. 배양용기의 바닥에는 멸균된 가는 모래(<0.5 ㎜)를 넣어주었다. 우화한 성충과 알집의 수는 24시간마다 기록하였다. 또한 각 세대별 알집의 수, 알집당 알의 수 및 부화율을 측정하였으며, 우화된 성충은 두폭, 가슴폭, 날개폭, 날개길이 및 체장을 측정하였다. 조각깔따구의 사육결과 2007년 4월부터 2008년 9월까지 11세대 이상의 계대사육이 진행되고 있으며, 이번 발표에서는 2008년 7월까지 진행된 6세대까지의 사육결과를 수집･측정하여 제시하였다. 우화한 성충의 수는 제 3세대에서 성충이 1,634개체로 가장 많이 우화하였고, 이 세대의 성충이 산란한 제 4세대의 알집이 486개로 가장 많았지만, 다음 세대부터 우화한 성충의 개체수와 알집 수는 점점 감소하였다. 우화한 성충의 성비(암컷 개체수/전체 개체수)가 0.24~0.46 범위를 나타냄으로써 암컷보다 수컷의 우화한 개체수가 모든 세대에서 더 많았다. 각 세대에 따른 알은 제 1세대에서 제 6세대로 진행될수록 첫 번째의 알집과 두 번째의 알집 모두 알집당 알의 수가 점점 감소하는 경향을 보였다. 알의 부화율은 세대에 따른 특성과 연관성을 보이지 않았으며, 첫 번째의 알집은 전 세대에서 80% 이상의 높은 부화율을 보였다. 성충의 크기 측정값은 세대가 진행됨에 따라 전반적으로 감소되는 경향을 보였다. 전체적으로 보아 조각깔따구의 사육 조건은 공간의 크기에 민감하며, 11세대가 지나는 시점에서도 비교적 번식력이 양호한 개체군이 유지되고 있다. 본 연구의 결과는 우리나라에서 깔따구류의 야생종을 순화하여 실내 계대사육을 시도한 최장기간(11세대)의 기록이며, 지금까지의 연구결과로 보아 조각깔따구는 수질평가를 위한 독성실험 등 다양한 환경생물학적 연구에 이용될 수 있을 것으로 사료된다.