주요 검역해충에 속하는 오리엔탈과실파리와 오이과실파리 두 종의 국내에서의 정착평가를 위한 기초자료를 작성하기 위하여 기존에 발표된 온도반응 데이터를 수집하여 개체군 모의에 필요한 단위모형의 매개변수를 구하였다. 이 값들로 두 해충의 개체군모형을 구성하고 제주도내 과거 기후관측 데이터를 이용하여 이들 해충이 각각 유입되었을 때의 연간 발생세대수와 그 양상을 모의하였다. 오리엔탈과실파리와 오이과실파리는 각각 연 3세대, 2세대가 발생하는 것으로 나타났으나 오리엔탈과실파리의 경우 번데기 상태로 월동하는 것을 고려해 볼 때 마지막 세대 성충의 실제 발생여부는 번데기의 온도반응 특성을 통해 검토할 필요가 있다.

본 연구의 목적은 시스템 식별기술과 칼만필터링을 기반으로 계측된 가속도응답을 이용하여 구조물에 작용하는 모달하중을 역으로 추정하는 기법을 제시하는데 있다. 하중추정 과정에서 가속도응답은 모드응답으로 분해되며 각 모드응답에 대한 모달하중을 역으로 추정한다. 역해석과정에서 발생하는 ill-posed문제를 피하기 위하여 가상진동기에 의해 구조물의 고유진동수와 감쇠비를 보다 정밀하게 구하기 위한 시스템 식별법을 사용하였다. 본 연구에서 제안된 기법의 적용성을 위하여 40층의 철골구조물에 태풍이 작용 했을 때 계측된 가속도를 사용하였다. 적용결과에 의하면 본 연구에서 제안된 하중추정법에 의하여 모달하중 추정이 가능함을 알 수 있었다.

곤충 발육기간분포모형은 어느 특정발육단계가 완료되는 시기의 변이를 나타내는 것으로 시간에 따라 발육이 완료되는 개체들의 비율로 나타낸다. 개체군 모형 작성시 발육단계전이모형의 기본 구성요소로 사용되며, 발육기간 누적빈도분포 자료를 Weibull 함수에 적용시켜 사용하고 있다. 곤충 발육기간분포모형을 추정하려면 발육실험의 원자료가 필요하기 때문에 원자료가 없는 경우 대체할 수 있는 방법의 개발이 필요하다. 따라서 본 연구에서는 곤충의 발육기간 빈도분포가 정규분포를 따른다는 가정 하에 해당 곤충의 발육기간 평균 및 표준편차를 이용하여 가상의 발육기간 빈도분포를 조성하고 매개변수를 추정하는 방법을 제시하였다. 정규분포에 준한 발육기간 누적빈도분포로 Weibull 함수의 모형 매개변수를 추정하고, 이것을 실제 발육기간 빈도분포를 이용하여 추정한 모형 매개변수와 비교하였다. 총 8종 (18개 발육단계) 곤충을 대상으로 분석한 결과 곤충 발육기간 빈도분포가 정규분포를 따르지 않음에 따라 원자료 값과 편차가 발생하였다. 일부 매개변수를 보정함으로써 원자료 값에 가깝도록 조절 할 수 있었으며, 이 방법을 통해 어떤 곤충의 발육기간 빈도분포자료 없이도 발육기간 평균과 표준편차만을 이용하여 발육기간분포모형을 추정할 수 있을 것으로 기대되었다.

건전성 예측은 구조물의 고장이 발생될 때까지 남은 시간인 잔존유효수명을 예측하는 것으로, 이는 안전 및 정비 계획과 직접적으로 연관되기 때문에 매우 중요하다. 건전성 예측방법에는 물리모델 기반방법, 데이터 기반방법과 두 방법의 장점 을 통합하는 방법이 있으며, 본 연구에서는 잔존수명 예측의 정확도가 모델변수 추정과 직접적으로 관련되는 물리모델 기 반 건전성 예측에 초점을 맞춘다. 물리모델기반 건전성 예측에서는 모델변수 추정을 통해 시스템 상태의 장기 예측이 가능 하지만, 대부분의 실제 구조물들의 상태모델은 여러 개의 모델변수를 포함함은 물론이고, 그 변수들이 서로 상관되어 있기 때문에 모델변수를 추정하는 일은 간단한 문제가 아니다. 본 연구에서는 물리모델 기반 건전성 예측을 위한 세 가지 변수 추정방법들의 차이를 논한다. 이 세 가지 방법들은 파티클 필터, 전반적인 베이지안 접근법, 그리고 순차적인 베이지안 접 근법으로 모두 베이지안 추론이라는 하나의 이론적 바탕에 기반하지만, 샘플링 방법이나 갱신 절차 등에서 차이가 있다. 균열성장을 표현하는 Paris 모델의 변수 추정을 통해 세 가지 방법의 차이점이 논해지고, 건전성 예측 메트릭을 이용하여 정량적 차이를 표현한다. 파티클 필터방법이 건전성 예측 메트릭 측면에서 가장 높은 성능을 나타내었지만, 전반적인 베이 지안 방법은 파티클 필터방법과 근소한 차이를 보이면서도 데이터가 집단으로 존재할 때에는 가장 효율적인 방법으로 나 타났다.

 ,  , Temperature-related parameters of Panonychus citri (McGregor) (Acarina: Tetranychidae) development were estimated and a stage-structured matrix model was developed. The lower threshold temperatures were estimated as 8.4℃ for eggs, 9.9℃ for larvae, 9.2℃ for protonymphs, and 10.9℃ for deutonymphs. Thermal constants were 113.6, 29.1, 29.8, and 33.4 degree days for eggs, larvae, protonymphs, and deutonymphs, respectively. Non-linear development models were established for each stage of P. citri. In addition, temperature-dependent total fecundity, age-specific oviposition rate, and age-specific survival rate models were developed for the construction of an oviposition model. P. citri age was categorized into five stages to construct a matrix model: eggs, larvae, protonymphs, deutonymphs and adults. For the elements in the projection matrix, transition probabilities from an age class to the next age class or the probabilities of remaining in an age class were obtained from development rate function of each stage (age classes). Also, the fecundity coefficients of adult population were expressed as the products of adult longevity completion rate (liiongevity) by temperature-dependent total fecundity. To evaluate the predictability of the matrix model, model outputs were compared with actual field data in a cool early season and hot mid to late season in 2004. The model outputs closely matched the actual field patterns within 30 d after the model was run in both the early and mid to late seasons. Therefore, the developed matrix model can be used to estimate the population density of P. citri for a period of 30 d in citrus orchards.

Poisson model and Gamma-Poisson model are popularly used to analyze statistical behavior from defective data. The methods are based on binary criteria, that is, good or failure. However, manufacturing industries prefer polytomous criteria for classifying

본 연구는 아스팔트 콘크리트 포장의 소성변형 추정을 보다 쉽게 하기 위하여 새로운 시험기법인 김테스트의 변형강도 및 시험 장비를 개발함에 있어 적정규격을 선정하기 위한 것이나. 김테스트에서 하중봉의 직경(D) 및 하중봉 하단의 원형처리 반경(r)이 시험결과에 미치는 영향을 파악하기 위하여 다양한 혼합물에 대하여 D와 r을 변화시켜가며 실험을 수행하였다. r에 따른 변형강도 및 변형하중과의 상관관계 분석결과 r=0.5와 1.0cm로 원형 처리한 측정값에서 소성변형과 높은 상관성을 보여주었다. 공시체의 직경(S)은 중요 변인이 아니었으며, 하중봉은 직경 4cm에 반드시 하단을 원형처리를 해야하고 이때의 절삭 반경 r은 1.0cm가 가장 좋은 것으로 나타났다. 통계프로그램 SAS의 STEPWISE 를 이용하여 골재별로 변형강도로부터 소성변형 깊이 및 동적 안정도를 추정하기 위해 모델을 개발하였으며 $R^2$은 0.95이상이 얻어졌다. 향후 보다 많은 실험을 통해 이 시험법의 표준화 연구가 지속적으로 이루어져야 할 것이며, 새로운 배합설계 방법을 개발한다면 본 실험법의 적용을 검토해 볼 필요가 있을 것이다.

본 연구는 아스팔트 콘크리트 포장의 소성변형 추정을 보다 쉽게 하기 위하여 새로운 시험기법인 김테스트의 변형강도 및 시험 장비를 개발함에 있어 적정규격을 선정하기 위한 것이나. 김테스트에서 하중봉의 직경(D) 및 하중봉 하단의 원형처리 반경(r)이 시험결과에 미치는 영향을 파악하기 위하여 다양한 혼합물에 대하여 D와 r을 변화시켜가며 실험을 수행하였다. r에 따른 변형강도 및 변형하중과의 상관관계 분석결과 r=0.5와 1.0cm로 원형 처리한 측정값에서 소성변형과 높은 상관성을 보여주었다. 공시체의 직경(S)은 중요 변인이 아니었으며, 하중봉은 직경 4cm에 반드시 하단을 원형처리를 해야하고 이때의 절삭 반경 r은 1.0cm가 가장 좋은 것으로 나타났다. 통계프로그램 SAS의 STEPWISE 를 이용하여 골재별로 변형강도로부터 소성변형 깊이 및 동적 안정도를 추정하기 위해 모델을 개발하였으며 R2은 0.95이상이 얻어졌다. 향후 보다 많은 실험을 통해 이 시험법의 표준화 연구가 지속적으로 이루어져야 할 것이며, 새로운 배합설계 방법을 개발한다면 본 실험법의 적용을 검토해 볼 필요가 있을 것이다.

This study investigates the performance of four Bayesian methods, Random Walk Metropolis (RWM), Hit-And-Run Metropolis (HARM), Adaptive Mixture Metropolis (AMM), and Population Monte Carlo (PMC), for estimating the parameters and uncertainties of probability rainfall distribution, and the results are compared with those of conventional parameter estimation methods; namely, the Method Of Moment (MOM), Maximum Likelihood Method (MLM), and Probability Weighted Method (PWM). As a result, Bayesian methods yield similar or slightly better results in parameter estimations compared with conventional methods. In particular, PMC can reduce parameter uncertainty greatly compared with RWM, HARM, and AMM methods although the Bayesian methods produce similar results in parameter estimations. Overall, the Bayesian methods produce better accuracy for scale parameters compared with the conventional methods and this characteristic improves the accuracy of probability rainfall. Therefore, Bayesian methods can be effective tools for estimating the parameters and uncertainties of probability rainfall distribution in hydrological practices, flood risk assessment, and decision-making support.

수문모형의 매개변수 추정에 필요한 유량 관측 자료의 수집은 시·공간적으로 제한이 있어 우리나라도 아직 상당수의 미계측유역이 존재하며, 이를 보완하고자 주변 유역의 정보를 활용하는 지역화 방법들이 연구되어 왔다. 그러나 지역적 특성이나 기후 조건에 따라 지역화 방법의 결과가 상이하여 어느 지역에 어떠한 지역화 방법이 가장 우수하다고 판단하기 어렵다. 본 연구에서는 보편적으로 사용되는 지역화 방법인 지역회귀모형의 설명 변수에 공간근접모형으로 추정한 수문모형의 매개변수를 추가하여 회귀모형의 적합성을 향상시켰으며, 이를 하이브리드 지역화모형이라 정의하고 기존 방법들과 비교하였다. 계측유역으로는 관측 자료가 충분한 남한의 37개 유역을 선정하였고, 수문모형은 개념적 수문모형인 GR4J를 사용하였으며, 계측유역에 대한 수문모형의 매개변수 산정은 Shuffled complex evolution 알고리즘을 사용하였다. 유역 특성변수들 간 다중공선성을 고려하기 위해 Variation inflation factor를 사용하였고, Stepwise regression을 통해 회귀모형의 최적 설명변수를 선택하였다. 통계 값을 통해 모형의 적합성을 비교한 결과, 하이브리드 지역화모형에서 가장 작은 RMSE 값을 나타내었으며, 유역별 모의 값의 변동성이 줄어들어 결과의 불확실성 또한 낮아짐을 확인할 수 있었다. 따라서 하이브리드 모형이 미계측유역의 유출량 산정을 위한 하나의 대안이 될 수 있음을 확인하였다.

본 연구에서는 미계측 유역에서 오염부하량 모의를 위해 LOADEST (LOAD Estimator) 기반 회귀모형의 최적 매개변수를 추정하고, 다중회귀분석 기법을 이용하여 유역특성에 따른 회귀 모의 모형의 매개변수 추정 방법을 개발하였으며, 개발된 모형의 적용성을 평가하였다. 오염부하량 모의 모형으로, T-N (Total-Nitrogen)은 LOADEST의 5번 회귀모형을, T-P (Total-Phosphorous)는 3번 회귀모형을 선택하였다. 모의결과, T-N, T-P 모두 선택된 회귀모형이 실측치를 잘 반영하였으나, 두 물질 모두 오염부하량이 과소 모의되어 실측치와 편의가 발생하는 것으로 나타나, 분위사상법을 이용하여 모의치의 편의보정을 실시하였다. 보정결과, 모형의 정확도는 크게 변하지 않았으나, 오염부하량이 과소 모의 되는 경향이 감소하는 것으로 나타났다. 다중회귀분석을 이용하여 회귀모형 매개변수와 유역특성간의 회귀식을 개발하였으며, 개발된 식을 평가한 결과, 실측치를 잘 반영하여 모의할 수 있는 것으로 나타났으며, 기존 매개변수에 의한 모의치와 유사한 모의능력을 갖는 것으로 나타났다. 본 연구에서 개발된 매개변수 추정방법은 실측자료가 확보되지 않은 소유역에 대한 오염부하량 모의와 정책결정을 위한 스크린 모델로서 활용할 수 있을 것으로 사료된다.

As hydrogeologic parameters such as hydraulic conductivity and specific yield are estimated by aquifer test , these are dependent on specific points at which field test was conducted. To overcome these site-specific limitations, a method of estimating transmissivity of aquifer using distribution features for parameters in Water table fluctuation model is newly suggested. Distribution features in reaction factor, specific yield and transmissivity having the function of pore space in aquifer are used to derive empirical equation for estimating transmissivity. From the result for applying the equation for 10 groundwater stations in Northeast Jeju Island, this equation is available for estimating transmissivity compared to the value estimated by existing equations. The estimated transmissivity ranged from 14.2 to 3,716.9 m2/day, and its average was 821.8 m2/day.

본 연구에서는 확률 대응법의 적용성을 평가하고, 이를 적용하여 실제 호우사상에 적합한 레이더 강우 추정관계식을 결정하였다. 먼저, 확률 대응법의 과정을 새로 정리하여 zT나 rT 등의 임계값을 결정하였다. 둘 째, 확률 대응법에 적용할 자료의 수가 매개변수 결정에 미치는 영향에 대하여 모의 자료를 이용하여 민감도 분석을 수행하였다. 셋 째, 확률을 나타낼 반사도와 강우강도 히스토그램의 구간 간격이 매개변수 결정에 미치는 영향에 대해 민감도 분석을 수행하였다. 그 결과, 확률대응법으로 매개변수를 결정할 경우 확률대응법에 2시간 정도 누적된 반사도와 강우 쌍을 적용하는 것이 적합함을 확인하였다. 아울러 확률대응에는 반사도와 강우강도의 누가확률 차이보다 1차 모멘트 차이를 이용하는 것이 양호한 결과를 얻을 수 있고, 30~100% 구간을 대응하는 것이 적합하다는 것을 확인하였다. 본 연구에서는 이러한 결과를 바탕으로 비슬산 레이더로 관측한 실제 호우사상에 대해 확률대응법으로 매개변수를 결정하였다. 또한, 매개변수를 이용하여 레이더 강우를 추정하였고, 이를 지상강우와 비교하였다. 전체적으로 레이더 강우가 지상강우보다 크게 산정되었지만, 레이더 강우와 지상 강우를 쉽게 대응하기 어려운 기존의 연구 결과를 고려하여 본 연구에서 평가한 확률대응법이 비교적 호우사상에 적합한 레이더 강우를 산정할 수 있다고 판단하였다.

본 연구에서는 극한강우를 고려한 산사태 및 토석류 재해 실시간 예측 및 대응을 위하여 필요한 원천기술들 중, 토석류 전이규준 개발 및 이를 통한 발생부 토사량 추정 기법, 그리고 전이 이후 토석류 확산의 정량적 예측 분석에 결정적으로 중요한 인자들인 연행 증가율 및 유변학적 인자들에 대한 입력값을 결정하는 기법에 대하여 연구하였다. 또한 본 원천기술들을 실시간으로 연계 및 적용 가능하도록, GIS(Geographic Information System) 및 Web 기반의 산사태 및 토석류 관련 D/B 시스템을 경기도 용인지역에 시범 구축하였다. 이 시스템은 클라우드 환경에서 산사태 및 토석류 유발인자 관련 지리 공간적 컨텐츠를 저장하고 게시하여 Web 상에서 D/B를 공동으로 활용하고 분석할 수 있다. 토석류 전이규준에서 발생부 토사량 추정으로 이어지는 연구는 두 가지 접근법을 통해 이루어졌다. 이는 지형학적 인자들을 바탕으로 통계학적인 규준값을 설정하는 접근법, 그리고 국내의 전형적인 사면파괴 형태인 강우침투로 인한 침윤전선(wetting depth)의 하강과 같은 실제 메커니즘을 고려하는 물리적 기반의 접근법 등이다. 한 편, 국내 과거 토석류 발생 지역 데이터들을 바탕으로 다중 회귀분석을 실시하여 DAN3D 모델을 통한 토석류 확산 예측 분석에서의 연행 증가율에 대한 추정식을 개발하였다. 본 연구결과들은 궁극적으로 국내의 실시간 산사태 재해 예·경보 및 종합적 리스크 평가 시스템 구축을 위한 핵심 구성 원천기술들로 활용될 수 있을 것으로 기대된다.

본 연구에서는 범용 매개변수 최적화 모형인 PEST를 이용하여 분포형 수문모형인 GRM(grid based rainfall-runoff model) 모형의매개변수및불확실성범위를추정하였다. 특히, 레이더강우및지상관측강우를각각적용하여, 입력자료 차이가매개변수추정에미치는영향을분석하였다. 자동보정모형은GUI (graphic user interface)에대한접근없이모형구동이가능하도록개선된GRM-MP (multiple projects) 버전과병렬PEST버전을결합하여매개변수추정에소요되는시간을단축시켰다. 이를낙동강수계금호강유역과감천유역에대해적용하여, 초기포화도, 지표면조도계수및토양투수계수의보정계수에 대해 매개변수 최적화 및 불확실성 추정을 수행하였다. 강우자료 분석 결과, 레이더와 지상 강우의 유역평균누적시계열은비슷하거나지상강우가조금큰경향을보였으나, 공간분포에있어서는지상강우에비해레이더강우에서큰 변동성이 확인되었다. 보정된 수문모의 결과는 레이더 강우 적용 시, 지상 강우에 비해 비슷하거나 더 나은 정확도를보였다. 추정된매개변수는레이더강우적용시, 토양투수계수의보정계수가일관되게1보다작은경향을보였으며, 이는강우강도가 강한 격자가 상당수 존재하기 때문으로 판단되었다. 초기 포화도 및 지표면 조도계수의 보정계수는 레이더및 지상 강우에서 일정한 경향성을 보이지 않았다. 본 연구의 대상 유역 및 호우사상에 대한 PEST의 최적화 모의 결과,동일유역및호우사상에대해서도강우추정방법에따라서로다른최적매개변수값을갖는것을알수있었으며, 이는향후 레이더 강우 자료의 수문 모의 활용 시 유의해야할 점으로 판단된다.

본 연구에서는 Clark 모형의 시간-면적곡선의 구성 방법과 적용성을 검토하고 모멘트 원리에 의한 도달시간, 저류상수를 합리적으로 산정하기 위한 방법론을 고찰해 보았다. 격자 기반으로 폭 함수를 구성하고 운동과정을 순수 이류현상으로 가정하여 시간-면적곡선으로 사용하였다. 또한 도달시간과 저류상수는 모멘트 법의 원리에 따라 Clark 모형 구조에 적용하여 해석적으로 산정할 수 있는 방법을 제시하였다. 적용성 검토를 위해 (1) HEC-1에서 기본적으로 제공하는 좌우 대칭형상인 무차원 시간-면적곡선을 적용하고 매개변수 산정은 관측유출수문곡선과 계산된 유출수문곡선의 오차를 최소화하는 HEC-1의 최적화 기법 사용, (2) HEC-1에 폭 함수 기반의 시간-면적곡선을 적용하고 매개변수 산정은 HEC-1의 최적화 기법 사용, (3) 폭 함수 기반의 시간-면적곡선을 이용하여 모멘트 원리에 따라 매개변수를 직접 산정하는 방법을 적용하였다. 방법별로 산정된 Clark 모형의 매개변수들을 HEC-1을 이용하여 직접유출량을 산정하고 관측 직접유출량과 비교하여 얻은 결과는 다음과 같다. (1) 정량적으로 비교하기 위해 산정한 첨두유량과 첨두발생 시간의 상대오차 및 효율계수 E(Efficiency Coefficient)를 비교한 결과, 시간-면적곡선을 폭 함수로 대체하여 HEC-1으로부터 추정된 매개변수가 관측값을 잘 반영하였다. (2) Clark 모형의 올바른 적용을 위해서는 HEC-1에서 기본적으로 제공하는 좌우 대칭형상인 무차원 시간-면적곡선보다는 적용 대상유역의 배수구조가 적절하게 반영된 시간-면적곡선의 사용이 합리적일 것으로 판단된다. (3) 본 연구 방법은 첨두유량과 첨두시간의 상대오차 범위와 재현정도를 나타내는 효율계수를 비교하여 볼 때 대체로 양호하게 모의되었고, 대상유역별 유량측정성과인 하천평균유속과 비교했을 때 본 연구 방법이 다소 실제 유속에 접근하고 있음을 확인하였다. (4) 본 연구에서 모멘트 원리를 기반으로 제안한 매개변수 추정을 위한 방법은 유역의 이류현상과 저류현상을 정량적으로 계량할 수 있는 효율적인 관계식으로 사용할 수 있음을 확인하였다. (5) 본 방법에 의해 계산된 수문곡선이 대부분 관측수문곡선의 우측으로 왜곡되고 첨두유량은 과소평가 되는 것을 보이고 있다. 이것은 평균과 분산만을 고려하여 유역을 하나의 평균이송속도로 모의한 본 연구의 한계점으로 판단된다. 만약 모멘트의 왜곡도를 고려하고 유역을 지표면과 하천으로 나누어 평균이송속도를 모의한다면 물리적인 특성을 충분히 반영하여 매개변수를 추정 할 수 있을 것으로 판단된다.

저류함수법의 최적 매개변수를 추정하기 위한 연구는 오랜 동안 여러 가지 방법으로 수행되어왔다. 그러나 여전히 최적 매개변수를 결정하는 것은 시간이 오래 걸리는 일이며, 유역의 물리적인 특성과 상관없는 매개변수가 결과로 제시되는 경우가 잦다는 인식이 팽배하다. 본 연구에서는 저류함수모형의 연속방정식과 저류함수식을 충실히 분석하고 민감도 분석을 수행하였다. 그 결과, 많은 수의 국지해 중에서 유일해를 결정하는 방법을 제안할 수 있었다. 또한 유역의 직접유출 시작 시간을 고려하여 저류함수법의 지체시간을 결정할 수 있다는 것을 보였으며, 매개변수의 민감도 분석 결과, 모형의 지체시간을 결정하는 것이 매우 중요하다는 것을 알 수 있었다. 지체시간을 결정한 후에는 유일한 해를 비교적 쉽게 찾을 수 있었다. 그러므로 제안된 방법은 기존의 최적화 방법과 같이 시간이 오래 걸리지 않으며, 강우사상별로 비교적 정확한 매개변수를 산정할 수 있다는 장점이 있다. 제안된 방법을 이용하여 기존의 저류함수법의 매개변수를 추정하기 위한 다양한 방법 중 상수고정법을 수정하였으며, 그 결과 실무에서 업무효율을 높일 수 있을 것으로 기대된다. 또한 제안된 방법은 기존의 유출수문곡선의 계산오차에만 의지하여 매개변수를 최적화하는 방법과는 다르게 유역의 특성을 고려할 수 있다는 점에서 그 의미가 있다.

본 연구는 우리나라 농업을 사례로 1993~2010년 기간 동안 전국 56개 시·군별 기후, 지리 및 토양, 사회·경제적 변수들에 대한 패널자료를 구축하고 리카디언 방법을 사용하여 기후변수가 논과 밭의 경작유형별 농지가치에 미친 경제적 효과를 추정하였다. 추정결과에 따르면, 기온변수가 논과 밭의 토지가치에 부정적인 경제적 영향을 미치지만 일조량은 토지가치에 긍정적인 영향을 미치는 것으로 나타났다. 이는 지구 온난화에서 강수량변수보다는 기온변수가 더욱더 중요하다는 기존의 연구결과를 재확인시켜 주고 있다. 다른 한편, 토지가치의 사회·경제적 결정요인들은 논과 밭의 경작유형에 따라 상이한 것으로 나타났다.

본 연구에서는 연최대치 독립 호우사상의 결정에 사용되는 Freund 이변량 지수분포의 매개변수 추정과정을 구체적으로 검토하였다. 먼저, 모멘트법을 이용하는 경우를 구체적으로 검토하고, 그 결과를 최우도법을 적용한 결과와 비교하였다. 두 방법을 1961～2010년 서울지점의 시강우 자료에 적용하여 연최대치 독립 호우사상을 선정하고, 그 결과를 비교·검토하였다. 이러한 과정을 통해 얻은 결과는 다음과 같다. 첫째, 매개변수 추정방법으로 모멘트법을 적용하는 경우에는 두 변량의 평균과 분산뿐만 아니라 상관계수도 고려해 주어야 하는 것으로 나타났다. 둘째, 최우도법은 두 변량의 평균에 대한 재현성이 우수하고, 모멘트법은 분산의 경년변동을 잘 나타내는 것으로 나타났다. 셋째, 모멘트법과 최우도법을 통해 선정한 연최대치 독립 호우사상들은 대체로 유사한 것으로 나타났다. 다르게 선정된 호우사상은 최우도법의 경우에는 총 강우량이 큰 것, 모멘트법의 경우에는 강우강도가 큰 것으로 나타났다.