The physical-based and lumped-parameter hydrologic groundwater flow model for predicting the rainfall-triggered rise of groundwater levels in hillside slopes is developed in this paper to assess the risk of landslides. The developed model consists of a vertical infiltration model for unsaturated zone linked to a linear storage reservoir model(LSRM) for saturated zone. The groundwater flow model has uncertain constants like soil depttL slope angle, saturated permeability, and potential evapotranspiration and four free model parameters like a, b, c, and K. The free model parameters could be estimated from known input-output records. The BARD algorithm is uses as the parameter estimation technique which is based on a linearization of the proposed model by Gauss -Newton method and Taylor series expansion. The application to examine the capacity of prediction shows that the developed model has a potential of use in forecast systems of predicting landslides and that the optimal estimate of potential 'a' in infiltration model is the most important in the global optimum analysis because small variation of it results in the large change of the objective function, the sum of squares of deviations of the observed and computed groundwater levels. 본 논문에서는 가파른 산사면에서 산사태의 발생을 예측하기 위한 수문학적 인 지하수 흐름 모델을 개발하였다. 이 모델은 물리적인 개념에 기본하였으며, Lumped-parameter를 이용하였다. 개발된 지하수 흐름 모델은 두 모델을 조합하여 구성되어 있으며, 비포화대 흐름을 위해서는 수정된 abcd 모델을, 포화대 흐름에 대해서는 시간 지체 효과를 고려할 수 있는 선형 저수지 모델을 이용하였다. 지하수 흐름 모델은 토층의 두께, 산사면의 경사각, 포화투수계수, 잠재 증발산 량과 같은 불확실한 상수들과 a, b, c, 그리고 K와 같은 자유모델변수들을 가진다. 자유모델변수들은 유입-유출 자료들로부터 평가할 수 있으며, 이를 위해서 본 논문에서는 Gauss-Newton 방법을 이용한 Bard 알고리즘을 사용하였다. 서울 구로구 시흥동 산사태 발생 지역의 산사면에 대하여 개발된 모델을 적용하여 예제 해석을 수행함으로써, 지하수 흐름 모델이 산사태 발생 예측을 위하여 이용할 수 있음을 입증하였다. 또한, 매개변수분석 연구를 통하여, 변수 a값은 작은 변화에 대하여 목적함수값에 큰 변화를 일으키므로 a의 값에 대한 최적값을 구하는 것이 가장 중요한 요소라는 결론을 얻었다.