본 연구는 우리나라 곰솔 현실림의 생장특성을 반영한 지위지수곡선을 개발하고자 최근 산림자원조사를 통해 획득한 493 plots의 곰솔 임분 자료를 수집하여 진행되었다. 수고생장모델 개발을 위해 생장추정에 널리 쓰이고 있는 Chapman-Richards, Schumacher, Gompertz 및 Weibull 대수차분방정식을 적용한 결과, 적합도가 가장 높은 Weibull식을 수고 생장모델에 대한 최적합식으로 선정하였으며, 유도곡선법을 이용하여 지위지수 곡선을 개발하였다. 개발된 해안⋅내륙 각각의 곰솔 지위지수 분류곡선을 이용하여 임령에 따른 독립된 두 집단의 우세목 수고생장을 비교 분석한 결과, 임령 50년 이후에서 유의적 차이가 있는 것으로 나타났다(p<0.05). 따라서 본 연구는 변화된 곰솔 현실림의 생장특성을 반영한 지위지수를 통해 지역 간 생장특성 차이를 구명함으로써 지역별 맞춤형 임분관리계획 수립 시 활용가치가 더욱 기대되는 바이다.

본 연구는 주요 난대 상록활엽수종인 붉가시나무, 구실잣밤나무, 동백나무, 황칠나무 등 4종 임분과 곰솔림을 대상으 로 2019년 5월~2020년 1월까지 각 산림의 임분 특성 및 NVOCs 발산특성을 알아보고자 수행되었다. 분석 결과, 각 조사구의 계절별 기온 및 습도는 여름철에 고온다습하고 겨울철 저온건조한 일반적 우리나라의 기후적 특성을 보였고 기압은 대체적으로 겨울이 높았으며 봄, 여름보다 가을과 겨울의 기압이 높았다. 조사구 5개소의 총 NVOCs는 전체적으로 여름철에 가장 높았으며 봄, 겨울은 황칠나무림, 가을은 구실잣밤나무림, 여름은 붉가시나무림의 NVOCs량 이 상대적으로 높은 특성을 보였다. 수종별로 살펴보면, 침엽수종인 곰솔림보다 난대 상록활엽수종인 붉가시나무, 구실잣밤나무, 동백나무, 황칠나무의 NVOCs 발산량이 낮지 않거나 오히려 높은 것을 확인하였다. NVOCs 발산에 기온(r=0.590, P=0.000), 습도(r=0.655, P=0.000)는 양(+)의 상관관계, 기압(r=-0.384, P=0.000)과 풍속(r=-0.263, P=0.018)은 낮은 음(-)의 영향이 있었다. 미기상인자 중 습도(β=0.507, P=0.000)가 NVOC 발산량에 가장 크게 영향을 미치는 것으로 나타났으며 기온, 기압, 풍속 순으로 영향을 미치는 것으로 분석되었다.

공주지역에 조림된 리기다소나무 군락에서 국내 대표 수종인 소나무와 방풍림으로 주로 조림되는 곰솔 낙엽의 분해율 및 분해과정에 따른 영양염류의 함량 변화를 파악하였다. 분해 60개월 경과 후 소나무 낙엽과 곰솔 낙엽의 잔존율은 각각 42.12±5.30과 24.79±1.98%로 소나무와 곰솔의 낙엽 분해율은 곰솔 낙엽의 분해가 소나무 낙엽의 분해에 비해 빠르게 일어났다. 60개월 경과 후 소나무 낙엽과 곰솔 낙엽의 분해상수 (k)는 각각 3.02과 3.59로 곰솔 낙엽의 분해상수가 다소 높게 나타났다. 소나무 낙엽의 분해과정에 따른 C/N, C/P 비율은 초기에 각각 14.4, 144.1 이었으나 60개월 경과 후에는 각각 2.26와 40.1로 점차 감소하였으며, 곰솔 낙엽의 경우 초기 C/N, C/P 비율은 각각 14.4와 111.3로 나타났고, 60개월 경과 후에는 각각 3.06와 45.8로 나타났다. 낙엽의 초기 N, P, K, Ca, Mg 함량은 소나무 낙엽에서 각각 3.07, 0.31, 1.51, 16.56, 2.03 mg g-1, 곰솔 낙엽에서 각각 3.02, 0.39, 0.99, 19.55, 1.48 mg g-1로 소나무 낙엽과 곰솔 낙엽의 질소와 인의 함량은 유사하였다. 60 개월 경과 후 N, P, K, Ca, Mg의 잔존율은 소나무 낙엽에서 각각 231.08, 130.13, 35.68, 48.58, 36.03%이었고, 곰솔 낙엽에서 각각 143.91, 74.02, 28.59, 45.08, 44.99%로 나타났다.

한반도해송숲의 지형환경, 기상조건, 토양 이화학적 성질과 곰솔의 생육상태를 조사해, 이곳 곰솔의 수세약화 원인을 찾아내고자 했다. 곰솔의 수세가 불량하거나 고사한 지역(수세불량지)과 대조구로서 곰솔의 생육상황이 양호한 지역(수세양호지)으로 나눠 조사했다. 그 결과, 수세불량지의 입목밀도는 수세양호지보다 유의적으로 높았지만, 수관폭은 반대로 수세양호지에서 더 넓었다. 고사목 개체수와 입목밀도는 서로 부(-)의 상관관계를 보였고, 또 입목밀도와 흉고직 경·수고·수관고(crown height)·수관폭 사이에서는 부(-)의 상관관계를 보였다. 즉, 입목밀도가 높으면 개체의 수관이 겹쳐 하부 가지가 죽어 수관고와 수관폭이 작아졌다. 이에 엽량 감소에 따른 광합성량이 줄어들어 수고와 직경의 생장량이 적어졌다고 본다. 이로 인해 수세가 약해지면서 해안사구라서 수분·양분의 부족과 강한 조풍 등의 환경압이 가중되어 해안가 근처의 곰솔이 대규모로 죽은 것으로 판단된다. 한반도해송숲의 고사목 발생과 수세약화 원인은 복합적이지만, 곰솔 생장에 따른 적정한 입목밀도 관리 부족이라는 내부요인을 제1 요인으로 꼽을 수 있다. 또 곰솔 생육에 불리한 토양조건과 조풍, 강풍 등의 환경압을 제2 요인이라고 볼 수 있다.

본 연구는 부산광역시 해안림 곰솔군락의 식생구조 및 생태적 특성을 파악하여 해안림의 식생구조를 파악하고 향후 부산광역시 해안림 관리방안을 위한 기초자료를 구축하고자 수행하였다. 부산광역시 해안림 지역에 100m² 크기의 조사구 97개를 설치하여 조사·분석하였다. TWINSPAN과 DAC기법을 사용하여 군락을 분리한 결과, 곰솔-졸참나무 군락, 곰솔-사스레피나 무(1)군락, 곰솔-사스레피나무(2)군락, 곰솔-갈참나무군락, 곰솔-동백나무(1)군락, 곰솔-동백나무(2)군락, 곰솔-사스레피나무- 동백나무군락으로 분리되었다. 조사결과 부산광역시 해안림 곰솔군락의 교목층에서는 곰솔이 주요 우점종으로 나타났으며, 하층식생에서 사스레피나무와 동백나무가 우점하였다. 곰솔-졸참나무군락에서 졸참나무가 일부 교목층과 아교목층에서 곰솔과 경쟁관계를 이루고 있으며 변화를 살피기 위해서는 지속적인 모니터링이 필요할 것이다. 해안림의 환경 특성상 곰솔을 비롯한 내염성이 강한 수종들이 군락을 이루고 있기 때문에 곰솔-졸참나무 군락을 제외한 군락들은 급격한 환경변화가 발생하지 않으면 현 상태의 군락 구조를 당분간 유지할 것으로 예상된다.

본 연구는 새만금 간척지에서 방풍시설의 설치에 따른 수림대의 초기 생존율을 분석하기 위하여 생존 분석 기법 중의 하나인 생명표법을 적용하였다. 생명표법은 의학, 사회과학 등 많은 분야에서 활용하고 있는 방법 이며, 본 연구에서는 방풍책의 설치에 따른 초기 수림대 수목의 생존 기간별 생존율, 중위 시간 등을 생명표 를 통해 알아보고, 생존함수 그래프와 위험함수 그래프 그리고 수림대별 비교를 통해 생존율의 차이를 분석 하고자 하였다. 주요 결과로는 곰솔 수림대의 중위 시간은 100% P-tb와 25% P-tb에서는 48.00개월 이었 고, 50% P-tb는 20.85개월로 짧았다. 25% P-tb에서 누적 생존율이 가장 높았으며, 위험률은 50% P-tb가 가장 높았다. 졸참나무 수림대의 경우 모든 수림대에서 중위 시간은 48.00개월 이었으며, 50% P-tb에서 누 적 생존율이 가장 높았고, 위험률은 100% P-tb에서 가장 높았다. 생존율 비교 결과, 곰솔 수림대는 전체적 으로 방풍책의 설치가 생존율에 영향을 미쳤으며, 특히 50% P-tb와 25% P-tb에서 유의미한 차이가 나타났 는데, 25% P-tb가 생존율과 위험률 측면에서 우수한 결과를 나타내었다. 반면에 졸참나무 수림대의 경우 전 체 수림대의 생존율이나 수림대별 생존율은 방풍책의 설치에 영향을 받지 않는 것으로 분석되었다.

곰솔의 줄기밀도, 현존량 확장계수, 상대생장식 및 임분 현존량을 파악하기 위하여 전라남도 여수지역에 위치하고 있는 15년생, 29년생, 45년생 등 3개 곰솔 천연임분을 대상으로 임분별 5주씩 총 15주의 표본목을 선정 벌목하여 조사하였다. 줄기밀도는 0.440-0.457g/cm³이었으며, 임분간 유의적인 차이는 없었다. 지상부와 뿌리를 포함한 임목 전체의 현존량 확장계수는 모두 임령이 증가할수록 감소하였으며, 임분 1은 지상부 현존량 확장계수에서 임분 2, 3과 유의적인 차이가 있었고, 임목 전체 현존량 확장계수에서 임분 3과 유의적인 차이가 있었다. 이것은 임령이 증가할수록 줄기의 건중량 구성비가 높아지기 때문이며, 임분 1의 경우 유령림의 생장특성을 보이기 때문으로 판단되었다. 흉고직경(D) 또는 흉고직경과 수고(D2H)를 독립변수로 하고 부위별 건중량(Wt)을 종속변수로 하는 2개 유형의 상대생장식(Wt=aDb, Wt=a(D2H)b)을 유도하고 적합도를 검정하였다. 임분 1, 2, 3의 지상부 현존량은 각각 50.72t/ha, 89.92t/ha, 194.07t/ha이었으며, 뿌리를 포함한 임목 전체 현존량은 각각 61.62t/ha, 113.12t/ha, 248.36t/ha이었다.

기후변화로 인한 곰솔림과 소나무림의 세력은 약해질 것으로 예상되고 있으나, 곰솔과 소나무에 대한 식재요구도는 여전히 높은 실정이다. 이에 본 연구는 한반도 난온대림지역에 속하는 부산일대, 남동부지역에 주요 조경수목인 소나무 림과 곰솔림을 대상으로 자연림의 생육밀도를 파악한 후 난온대림에 적합한 생태적 식재밀도를 제시하고자 하였다. 연구는 한반도 남동부에 해당하는 부산광역시를 중심으로 소나무림 22개소, 곰솔림 60개소를 선정하여 100m² 방형구 를 설치하여 조사구 내 출현하는 교목, 아교목층 수종의 흉고직경 및 개체수를 조사한 후, 회귀분석을 통해 평균흉고직 경급별 적정 개체수 산정을 위한 회귀식을 도출하였다. 회귀분석 결과, 흉고직경의 생장에 따른 교목층 생육개체수 예측모형의 결정계수는 곰솔림 0.700, 소나무림 0.533으로 단일 요인에 의한 설명력으로는 매우 높은 것으로 분석되었 다. 예측 회귀모형은 곰솔림의 경우 Ｙ= 31.176ｅ -0.055χ , 소나무림의 경우 Ｙ= 38.351ｅ -0.059χ (X=흉고직경, Y=개체수 /100m²)였다. 남동부지역의 곰솔림과 소나무림의 생육밀도는 중부지방보다 다소 높았다.

본 연구는 제주도 저지대에 분포하고 있는 곰솔군락 하층 식생의 구성종을 파악하여, 한국의 난온대지역인 제주도에서 상록침엽수림이 천이가 진행되면 원래의 잠재식생인 상록활엽수림으로 바뀌는지를 알아보고자하였다. 식생조사 는 곰솔우점군락 65지점의 조사결과를 이용하였으며, 하층 식생(T2, S, H)을 구성하는 목본종의 계층별 상대피도(RFi, %), 상대빈도(RCi, %) 그리고 중요치(importance value)를 활용하였다. 상대피도는 아교목층이 까마귀쪽나무, 곰솔, 후박나무, 사스레피나무, 참식나무 순으로, 관목층은 까마귀쪽나무, 사스레피나무, 후박나무, 예덕나무, 참식나무 순으로, 초본 층은 까마귀쪽나무, 참식나무, 예덕나무, 우묵사스레피나 무, 후박나무 순으로 높았다. 상대빈도는 아교목층이 후박나무, 사스레피나무, 참식나무, 예덕나무, 까마귀쪽나무 순으로, 관목층은 까마귀쪽나 무, 예덕나무, 참식나무, 후박나무, 팽나무 순으로, 초본층은 까마귀쪽나무, 참식나무, 예덕나무, 우묵사스레피, 후박나무 순으로 높았다. 중요치는 아교목층이 후박나무, 사스레피나무, 곰솔, 까 마귀쪽나무, 참식나무 순으로, 관목층은 까마뀌쪽나무, 사스레피나무, 후박나무, 예덕나무, 천선과나무 순으로, 초본 층은 까마귀쪽나무, 후박나무, 참식나무, 사스레피나무, 예덕나무 순으로 높았다. TWINSPAN에 의한 군락구별은 Level 1에서 사스레피 나무(+), 까마귀쪽나무·우묵사스레피(-)로, Level 2은 까마 귀쪽나무(+)와 천선과나무·예덕나무·돈나무(-) 및 종가시시 나무(+)로, Level 3은 우묵사스레피(+), 구실잣밤나무(+), 예덕나무(+)와 팽나무(-), 후박나무(-)의 출현에 의해 총 8개 군집으로 구별되었으며, CCA 분석을 통해 해발고도, 경사, 해안으로부터의 거리 그리고 교목층의 피도가 곰솔군락의 하층식생 구성종에 영향을 미치는 것으로 생각된다. 향후, 곰솔 군락은 앞으로 난온대 상록활엽수인 후박나 무, 까마귀쪽나무, 참식나무, 사스레피나무 등으로의 천이 가 예측된다.

재선충속(Bursaphelenchus) 선충의 식물 기주체내 증식 능력 및 병원성을 비교하기 위하여 선충 3종 (이하 Bursaphelenchus xylophilus: Bx, B. mucronatus: Bm, B. thailandae: Bt)과 해송(Pinus thunbergii), 잣나무(Pinus koraiensis) 2종을 실험에 사용하였다. 5가지 선충 처리 조합 (①Bx+Bt, ②Bx+Bm, ③Bx+Bt+Bm, ④Bt+Bm, ⑤Bm)으로 한 본당 10,000마리씩 접종 한 후 7주 동안 외부 병징 및 선충 밀도를 조사하였다. 그 결과, Bx가 포함된 대부분 처리구에서는 5주 차에 전체 고사가 진행 되었고, Bx가 포함되지 않은 처리구는 병징이 나타나지 않았다. 예외적으로, Bx+Bm, Bx+Bt+Bm 해송 처리구에서는 7주 후에도 전체고사가 나타나지 않았다. 이 경우 선충 밀도 조사 결과, Bm이 약 42%를 차지하였으며, 이는 소나무재선충병의 병징 발현을 약화를 유도하는 역할을 하는 것으로 추정된다. 모든 처리구에서 Bx의 밀도가 가장 높았으며, Bt는 전혀 검출되지 않았다.

In Mt.Baekun in Yeongjong island, there is Pinus thunbergii community, so there are many insect species. We caughtinsects to investigate entomofauna in dead pine tree. The tree is about 6 months later after dead, so fungis grow there.And it contains much moisture and it is located about 100 altitude. Caught insects were treated by ethanol, and larvaeswere made into immersion specimen, and the others were made into exsiccata. Also, we divided the tree in 5 parts inheight, and 3 parts in radius. The most frequently appeared species were cerambycidaes. Many species appeared suchas Peseudo japonicum, insularis Motschulsky. Moreover, we confirmed many insects inhabit there, such as Pentatomidaes,Elateridaes. And insect distribution changes by height and radius. This shows even just one dead tree forms a small ecosystem,and suggests necessity of environmental preservation.

본 연구는 해안림의 식생특성을 분석함으로써 다층구조의 혼효림으로 발전할 수 있는 가능성을 파악하고, 관리방안을 제시할 목적으로 수행되었다. 본 연구에서는 서해안에 위치한 기지포, 춘장대, 장항 및 구시포 등 4개 지역의 해안림을 대상으로 식생구조와 식물상을 조사 분석하여 다음과 같은 결과를 얻었다. 4개의 조사지역 교목층에서 곰솔이 우점하고 있으며 조사지역에 따라 일부 교목층과 아교목층에서 소나무, 아까시나무, 졸참나무, 산벚나무, 상수리나무 및 밤나무가 우점하고 있었다. 4개 지역 해안림의 곰솔군락은 전반적으로 어린 개체와 큰 개체의 밀도보다 중간 개체의 밀도가 높아 곰솔군락으로 계속 유지될 것으로 판단된다. 구시포 해안림의 경우는 지형적인 입지여건에 따라 장기적으로 참나무류의 군락으로 천이가 진행될 가능성이 있다. 조사대상지에서 총 205종의 식물 종이 조사되어 다양성이 높은 것으로 나타났다. 다층구조의 해안림으로 유도하기 위해서는 곰솔림의 간벌, 통행로 설치, 이용구역 제한 등의 관리를 통해 다양한 식생이 침입하여 생육할 수 있는 환경을 조성할 필요가 있다.

Nondestructive evaluation (NDE) technique method using a resonance frequency by free vibration mode was carried out for damaged and sound black pine woods( Pinus thunbergii ) by Matstucoccus thunbergianae in order to utilize damaged trees. There were a close relationship of dynamic modulus of elasticity and static bending modulus of elasticity to MOR. It was found that there were a high correlation at 1% level between dynamic modulus of elasticity and MOR, and static modulus of elasticity and MOR in all conditions. However, the result indicated that correlation coefficient is higher in dynamic modulus of elasticity to MOR than that in static modulus of elasticity to MOR. Therefore, the dynamic modulus of elasticity using resonance frequency by free vibration mode is more useful as a nondestructive evaluation method for predicting the MOR of sound and damaged woods by Matstucoccus thunbergianae.

본 연구의 목적은 우리나라 남부지역 해송에 적합한 지위지수를 개발하는데 있다. 남부지역에 분포하는 해송림을 대상으로 우세목과 준우세목의 수고와 임령을 측정하였다. 5개의 생장모형을 비교하고 수집된 자료에서 40%는 모형의 개발, 나머지는 모형의 검정에 이용하였다. 지위지수식의 개발을 위한 생장모형의 적합도는 적합도 지수, 편차, 평균오차제곱근을 이용하였다. 남부 지방 해송의 지위지수 추정을 위하여 5종의 생장모형을 비교한 결과, Gompertz 식에서 높은 적합도를 보였으며 적합도 지수, 편차, 평균오차제곱근은 각각 0.618, -0.281, 2.266으로 나타났다. Gompertz 모형을 이용하여 지우지수를 유도하고 기준임령은 해송 벌기령을 고려하여 50년으로 정하였다. 남부지방 해송의 지위지수(기준임령 50년)는 16에서 24의 범위를 나타내었다. 본 연구는 남부지역 해송의 경영에 기초 정보를 제공할 것으로 판단된다.

우리나라 남부권역에 분포하는 곰솔 임분의 식생, 생장, 형질을 권역 및 건전도에 따라 임분 특성을 분석하여 시업적 기초 자료를 제공하고자 본 연구를 수행하였다. 아교목층 이하에서는 교란임지가 건전임지 보다 많은 종이 출현 하였으며, 종 다양도 지수는 교목층을 제외한 모든 층위에서 건전임지(0.097 ∼1.520)와 교란임지(1.507∼1.669) 사이에는 유의적인 차이를 보였다. 교란임지와 지역별 유사도분석에서, 경상권과 전라권이 층위별로 0.280∼0.607, 제주권역은 0.189∼0.318로 지역 간 차이를 나타내었다. 연년생장량은 정상임분이(0.39∼0.48 ㎝/year) 교란임분(0.25∼0.37 ㎝/year)보다 높게 나타났으며, 생장율도 건전임지(7.74∼8.51%)와 교란임지(3.79∼5.52%)간에는 유의적 차이를 보였다. 임목 형질은 경상권과 전라권에서 B급 이하, 제주권은 B급 이상, 교란지에서는 C급의 비율이 상대적으로 높게 나타는 것으로 조사되었다. 따라서 곰솔 임분은 현재의 식생, 생장, 형질적 특성을 고려하여 차별화된 맞춤형 임분 관리가 필요할 것으로 판단되었다.

본 연구의 목적은 6개의 수간곡선을 비교하여 우리나라 남부지역 해송에 적합한 수간곡선식을 개발 하는데 있다. 59본의 입목을 벌채하여 1년 단위 수간석해를 통해 자료의 전체 수고, 단면고별 수피 외직경, 단면고별 수피 내 직경, 임령을 측정하였다. 수간곡선식의 적합도는 적합도 지수(fit index, FI), 편차(bias), 평균오차제곱근(root mean square error, RMSE)을 이용하였다. 남부 지방 해송의 수간곡선 추정을 위하여 6종의 수간곡선식을 비교한 결과, Kozak(1988) 식에서 높은 적합도를 보였으며 FI, Bias, RMSE는 각각 0.9802, 0.0284, 1.2165로 나타났다. 본 연구의 주요한 결과로 추정된 수간곡선식은 남부지방의 주요 수종인 해송의 생장 정보 제공이 가능한 것으로 판단된다.

본 연구는 관광객의 증가 및 개발로 인해 식생의 파편화가 진행되고 있는 해안 곰솔림을 대상으로 효율 적인 보전과 관리를 위한 기초자료를 제공하기 위하여 곰솔 우점림에 조사구 70개(20×20m)를 설치하여 수행하였다. 경상북도 동부 해안가 곰솔림의 군집분석 결과 곰솔림은 3개의 군집(곰솔-소나무군집, 곰솔-아까시나무군집, 곰솔-소나무-팽나무군집)으로 구분되었으며, 다수응답순열절차 검정 결과, 군집은 유의성 있게 나타났다. 지표종 분석 결과, 유의한 종은 소나무, 작살나무, 노간주나무, 산철쭉 등 19종이 분석되었다. 곰솔림의 개체군 구조는 12∼26 ㎝의 직경급이 우점하는 것으로 나타났으며 고봉형 양상(hump-shape pattern)으로 나타났다. 종다양도(H')는 목본층에서 1.033±0.234∼1.629±0.226, 초본층에서 2.448±0.232∼2.545±0.318 사이로 나타났다.

본 연구는 대기오염으로 인한 SO2와 NO2 그리고 산성비 등의 건성 및 습성강하물의 영향을 받고 있는 공단지역 주변 산림과 대조지역의 해송림에서 해송림의 쇠퇴도 사이에 상관분석과 단계적 회귀분석을 통하여 해송 잎의 쇠퇴와 엽록소 함량에 가장 영향을 미치는 인자를 구명하였다. 수관통과우의 pH는 도시공단지역이 봄철 4.65, 가을철 4.72이었고, 대조지역이 봄철 5.32과 가을철 5.34였다. 음이온 성분 중에서 NO3-의 농도는 도시공단지역이 52.13 ㎎/ℓ, 대조지역은 37.85 ㎎/ℓ이었다. SO42- 함량은 도시공단지역이 57.89 ㎎/ℓ, 대조지역이 36.21 ㎎/ℓ이었다. Ca2+의 농도는 도시공단지역이 27.27 ㎎/ℓ, 대조지역이 9.48 ㎎/ℓ이었다. 해송 잎의 엽록소 a의 함량은 도시공단지역이 0.2378이었고, 대조지역은 0.4378이었다. 해송의 쇠퇴도는 도시공단지역 2.97, 대조지역 1.20이었다. 해송 임분의 쇠퇴도와 강우의 이온성분, SO2, NO2 농도와 상관분석 결과에서 강우산도 상호간에는 부의 상관(r=-0.8672)이었고, 대기 중 SO2 농도는 정의 상관(r=0.8924)이었으며, NO2 농도도 정의 상관(r=0.8428)이었다. 이러한 산성강하물 상호간의 상관관계는 pH(r=-0.8672), NO3-(r=0.6996), SO42-(r=0.8497), SO2 (r=0.8924), NO2(r=0.8428)는 1% 수준에서 상관이 있었다.