석회암 공동이 발달한 도심지역에서 지하수위 하강에 수반되어 발생한 함몰형 지반침하의 원인 규명 및 공동의 분포 특성 파악을 위하여 시추공을 이용한 전기비저항 토모그래피탐사를 실시하였다. 이때 지하수 수리지질 특성을 파악하기 위하여 시추코아의 비저항 측정, 지하수위 측정 및 수리전도도 해석을 병행하였다. 연구 지역에서의 완만한 지하수위 분포 특성과 0.8-9.3×10-4 cm/s 범위의 수리전도도 분포로 부터 연구지역의 수리지질 특성은 불균질성이 크지 않은 것으로 나타났다. 시추코아를 이용한 전기비저항 측정 결과 연구지역의 석회암은 파쇄가 많은 경우, 변질이 심한 경우 및 신선한 경우로 나눠지며, 전기비저항은 각각 103-161, 218-277 및 597-662 ohm-m의 범위로 나타났다. 시추결과 점토로 충전된 석회암 공동 지점은 토모그래피 탐사자료의 역산 결과 50 ohm-m 이하의 낮은 비저항으로 나타났으며, 각 시추공 간 비저항 영상 단면으로부터 연구지역 전체적으로 지표 하부 심도 약 10-20 m 구간까지 파쇄대 또는 석회암 공동 구간이 분포하는 것으로 나타났다. 또한 석회암 공동의 직경은 약 4-6 m 규모로, 대부분 점토질로 충전된 것으로 판단된다.

본 연구에서는 죽산천 유역에 대해 지표수와 지하수의 통합거동을 유역 스케일로 장기간 모의할 수 있고 양수정의 공간분포 및 저수지 운영을 복합적으로 고려할 수 있는 유역단위 통합수문해석모형 SWAT-MODFLOW를 적용하여 지하수 양수 및 농업용 저수지가 하천유량에 미치는 영향을 평가하였다. 지하수 양수 및 농업용 저수지 고려 유무에 따라 4가지 시나리오((1) 현 지하수 이용량 고려, 저수지 고려, (2) 현 지하수 이용량 고려, 저수지 미고려, (3) 지하수 이용량 미고려, 저수지 고려), (4) 지하수 이용량 미고려, 저수지 미고려(자연상태)) 를 구성하고 각 시나리오별로 하천유량의 변화를 모의한 결과, 죽산천 유역 출구부를 기준으로 지하수 양수로 인한 영향이 농업용 저수지에 의한 영향 보다 상대적으로 크게 발생하였으며, 갈수량은 용수이용이 없는 자연 상태(natural flow)에 비해서 지하수 양수만의 영향으로 약 17%, 저수지만의 영향으로 약 11%, 양수 및 저수지의 복합영향으로 인해 약 23% 만큼 감소하는 것으로 분석되었다. 농업용 저수지로부터 떨어진 거리에 따른 하천유량의 변화를 평가하기 위해서 죽산천 본류를 따라 비갈수량을 산정한 결과, 저수지 관개용수의 회귀수 영향으로 하류로 내려갈수록 비갈수량이 증가하였으며, 죽산천 상류에 위치한 덕산 저수지로부터 약 6.5km까지는 저수지의 영향이 상당히 미치는 것으로 분석되었다.

본 연구는 SWAT모형을 이용하여 기상, 지하수 취수량, 토지이용 변화에 대한 건기 총유출량의 민감도를 제시하였으며 더 나아가 건기의 총 유출량을 추정하기 위해 건기 총강우량, 전 우기 총강우량, 평균 일 최대온도, 일평균 태양복사량과 같은 기상 변수들과 지하수 취수량 및 도시면적 비율을 이용하여 회귀식을 도출하였다. 도출된 식을 이용하여 기후변화에 대한 건기 총유출량 변화를 살펴보기 위하여 온도와 강우량의 변화에 대한 건기 총유출량의 변화율을 제시하였