• 한국대기환경학회지
• /
• 제28권3호
• /
• pp.233-248
• /
• 2012

The national emission from energy sector accounted for 84.7% of all domestic emissions in 2007. Of the energy-use emissions, the emission from mobile source as one of key categories accounted for 19.4% and further the road transport emission occupied the most dominant portion in the category. The road transport emissions can be estimated on the basis of either the fuel consumed (Tier 1) or the distance travelled by the vehicle types and road types (higher Tiers). The latter approach must be suitable for simultaneously estimating $CO_2$, $CH_4$, and $N_2O$ emissions in local administrative districts. The objective of this study was to estimate 31 municipal GHG emissions from road transportation in Gyeonggi Province, Korea. In 2008, the municipalities were consisted of 2,014 towns expressed as Dong and Ri, the smallest administrative district unit. Since mobile sources are moving across other city and province borders, the emission estimated by fuel sold is in fact impossible to ensure consistency between neighbouring cities and provinces. On the other hand, the emission estimated by distance travelled is also impossible to acquire key activity data such as traffic volume, vehicle type and model, and road type in small towns. To solve the problem, we applied a hierarchical cluster analysis to separate town-by-town road patterns (clusters) based on a priori activity information including traffic volume, population, area, and branch road length obtained from small 151 towns. After identifying 10 road patterns, a rule building expert system was developed by visual basic application (VBA) to assort various unknown road patterns into one of 10 known patterns. The expert system was self-verified with original reference information and then objects in each homogeneous pattern were used to regress traffic volume based on the variables of population, area, and branch road length. The program was then applied to assign all the unknown towns into a known pattern and to automatically estimate traffic volumes by regression equations for each town. Further VKT (vehicle kilometer travelled) for each vehicle type in each town was calculated to be mapped by GIS (geological information system) and road transport emission on the corresponding road section was estimated by multiplying emission factors for each vehicle type. Finally all emissions from local branch roads in Gyeonggi Province could be estimated by summing up emissions from 1,902 towns where road information was registered. As a result of the study, the GHG average emission rate by the branch road transport was 6,101 kilotons of $CO_2$ equivalent per year (kt-$CO_2$ Eq/yr) and the total emissions from both main and branch roads was 24,152 kt-$CO_2$ Eq/yr in Gyeonggi Province. The ratio of branch roads emission to the total was 0.28 in 2008.

• Journal of Preventive Medicine and Public Health
• /
• 제31권3호
• /
• pp.564-578
• /
• 1998

본 연구는 우리나라 노인요양시설의 종사자 배치현황을 검토하고, 노인요양시설에 수용되어 있는 노인들의 건강상태와 돌봄 필요시간을 측정함으로써 건강상태와 돌봄필요도에 따라 필요한 돌봄인력의 추정에 기초자료로 활용될 수 있도록 하고자 하였다. 본 연구의 대상은 강원도 2개 요양시설과 서울 2개 요양시설에 수용되어 있는 187명의 노인에 대하여 직접면접조사를 실시하였고, 159 명에 대하여 면접이 이루어져 85%의 조사율을 보였다. 각 시설의 시설장 및 총무를 대상으로 시설인력현황을 조사표를 통하여 직접 작성하게 하였다. 강원도와 서울 일부 요양원에 입소하고 있는 노인들의 인구사회학적 특성에 대해서는 본 대학원 석사과정 학생 6명이 노인들을 대상으로 직접 면접조사를 실시하였고, 돌봄 소요시간의 경우, 입소노인의 돌봄을 담당하는 간호(조무)사 및 생활보조원을 대상으로 조사표를 통하여 직접 작성하게 하였다. 이상의 연구방법을 통하여 다음과 같은 결론을 얻었다. 1. 조사대상 4개 요양시설 모두 노인복지시설의 인력기준에 미달하고 있다. 2. 노인의 건강평가도구 중 ADL이 노인의 총 돌봄 필요시간과 가장 높은 상관성을 보여주고 있고, 다중 회귀분석의 결과, 돌봄 소요시간에 영향을 주는 유의한 변수는 단지 일상생활 수행능력 하나였고, 일상생활 수행능력 점수가 클수록, 즉 노인의 의존도가 클수록 돌봄 소요시간이 커지는 경향이 있었다. 3. 개별 입소노인의 ADL 총 점수를 네 그룹으로 나누었을 때, 그룹별 평균 돌봄 소요시간은 I군에서 IV군이 각각 하루에 15.3분, 21.1분, 36.7분, 88.8분의 돌봄서비스를 요구하고 있다. 4. 노인요양시설별로 ADL 그룹별 환자의 분포를 보면, 전체적으로 II군에 속하는 노인은 40.0%로 가장 많았으며, I군, III군, IV군의 순을 보여주고 있다. 시설별로 ADL 그룹별 환자의 수는 유의한 차이를 보여주지 못하고 있다. 5. 모든 요양시설의 간호사와 보조원 수는 입소노인의 건강상태를 반영하여 추정한 적정인력에 미치지 못하고 있다. 그러나, 법정기준에 의한 간호사와 보조원수는 노인의 돌봄필요도를 반영하여 추정한 수 보다 높게 나타나고 있어 본 연구의 조사대상 요양시설과 비슷한 건강분포를 보이는 노인요양시설의 경우에는 법정기준에 부합하다고 할 수 있다. 그리고, 노인복지사업지침에 근거한 예산지원기준을 적용시킬 때는 모든 요양시설에서 적정인력에 미치지 못하고 있음을 알 수 있다. 6. 무료노인요양시설의 전문인력을 확보하기 위해서는, 무엇보다도 사회복지 시설 종사자에 대한 처우 개선방향과 상응하는 방향으로, 급여수준을 상향조정하여야 할 것이다. 7. 노인들의 요양필요도를 실제적으로 파악함에 있어 ADL을 이용한 건강지표가 장기적으로는 노인요양 시설의 단순 입소판정기준으로 사용될 수 있을 것이다.

• 한국토양비료학회지
• /
• 제37권2호
• /
• pp.74-82
• /
• 2004

토양의 질산태 질소 함량이 46에서 $344mg\;kg^{-1}$으로 분포되는 시설재배지 12개 토양으로 시비수준을 무비구와 표준시비구로 구분하여 방울토마토를 공시 작물로 한 생산력 검정시험을 포트재배로 수행하였다. 무비구의 토마토 건물중은 최소 $28.9g\;plant^{-1}$에서 최대 $112.5g\;plant^{-1}$으로 토양 비옥도 차이에 의해 다양한 생산능력을 보였으며 무비구 토마토의 양분 흡수량도 동일한 경향을 나타냈다. 토양의 질산태 질소 함량은 무비구의 토마토 건물중 및 질소 흡수량과 고도로 유의한 정의 상관 ($r=0.83^{**}$ 및 $r=0.78^{**}$)을 보였고, 표준시비구와 무비구의 건물중과 질소흡수량 차이로 평가한 비료효과 및 질소시비효율과는 부의 상관을 보였다. 토양의 질산태 질소 수준에 따른 무비구 토마토의 건물중. 질소흡수량, 비료효과. 질소시비효율과의 상호관계 분석으로부터 질소 무시비 재배를 위한 토양의 질산태 질소 상한기준은 $280mg\;kg^{-1}$, 질소 표준시비 재배를 위한 토양의 질산태 질소 하한기준은 $50mg\;kg^{-1}$ 이하로 추정되었다. 이러한 질산태 질소 한계기준은 시설재배 토마토의 질소시비 추천식에서 전기전도도의 한계기준을 질산태질소로 환산한 기준과 거의 동일하였다. 따라서 상기의 질산태 질소 상한기준과 하한기준으로부터 토양의 질산태 질소 함량이 $50-280mg\;kg^{-1}$ 사이인 경우 질소시비량은 다음과 같은 직선모델식으로 추천하였다. Y (질소 표준시비량에 대한 시비비율, %) = -0.4348 X 토양 중 $NO_3-N$($mg\;kg^{-1}$) + 121.74

• 수산해양기술연구
• /
• 제27권1호
• /
• pp.21-30
• /
• 1991

어군탐지기를 이용하여 어군량을 추정할 때 필요한 어체의 평균반사강도를 정량적으로 계측하기 위한 디지털측정시스템을 제작하고, 이의 성능을 ogive와 장타원체 표적을 대상으로 검정하는 한편, 실제 8개 어종에 대한 체장별, 체중별, 주파수별 반사강도를 측정하여 반사강도의 평균화에 대하여 검토한 결과를 요약하면 다음과 같다. 1. 본 연구에서 제작한 어체의 반사강도 디지털측정시스템은 ogive와 장타원체 표적을 대상으로 이론 및 실험적인 산란특성을 비교, 검정한 결과, 양자가 서로 잘 일치하여 정도가 매우 높다고 판단되었다. 2. 8개 어종에 대한 어체의 평균반사강도는 말쥐치, 병어, 갈전갱이, 전갱이, 민어, 강달이의 경우, 평균체장이 각각 16.1cm, 28.8cm, 17.9cm, 17.6cm, 32.0cm, 15.5cm이고, 평균체중이 각각 62.0g, 475.2g, 110.0g, 58.8g, 319.5g, 39.0g였을 때 25kHz에서 각각 -40.8dB, -39.9dB, -41dB, -36.9dB, -35.7dB, -41.1dB였고, 100kHz에서는 각각 -43.0dB, -45.8dB, -45.9dB, -41.9dB, -39.1dB, -45.2dB였다. 또, 50kHz에 대한 말쥐치, 민어, 노래미, 넙치의 평균반사강도는 각각 -44.5dB, -40.6dB, -42.0dB, -42.2dB이었다. 3. 어종별 어체의 체장(L, cm)과 평균반사강도(TS, dB)와의 관계를 TS=A+20Log(N)의 회귀직선식에 근사시켰을 때, 어종별에 대한 회귀직선식의 계수 A 값은, 말쥐치, 병어, 갈전갱이, 전갱이, 민어, 강달이의 경우, 25kHz에서 각각 -65.1dB, -69.4dB, -67.5dB, -62.7dB, 65.7dB, -65.5dB였고, 100kHz에서는 각각 -67.2dB, -75.9dB, -71.0dB, -67.9dB, -69.5dB, -69.0dB였다. 또, 50kHz에 대한 말쥐치, 민어, 노래미, 넙치에 대한 A 값은 각각 -69.1dB, -70.1dB, -65.5dB, -68.6dB이었다. 4. 8개 어종에 대하여 어체의 중량 1kg당에 대한 평균반사강도는 말쥐치, 병어, 갈전갱이, 전갱이, 민어, 강달이의 경우, 25kHz에서 각각 -32.1dB, -37.7dB, -35.9dB, -29.3dB, -32.2dB, -32.0dB였고, 100kHz에서는 각각 -34.3dB, -44.2dB, -39.4dB, -34.5dB, -35.9dB, -35.6dB였다. 또, 50kHz에서 말쥐치, 민어, 노래미, 넙치에 대한 그 값은 각각 -37.2dB, -37.3dB, -31.9dB, -36.0dB였다. 한편, 어체의 평균반사강도의 추정정도는 자원량추정치의 신뢰도를 좌우하는 직접적인 요인이 되므로, 앞으로 이 문제에 대해서는 보다 체계적인 연구가 필요하다고 판단된다.

• 대한원격탐사학회지
• /
• 제30권2호
• /
• pp.275-292
• /
• 2014

이 연구에서는 Geostationary Ocean Color Imager(GOCI) 센서에 적용할 수 있는 고유의 Tasseled Cap Transformation(TCT) 계수를 제시하고 있다. TCT는 다중밴드 센서 자료로부터 지표의 특성을 분석하는 전통적인 영상변환 방법 중 하나로 새로운 다중밴드 광학센서가 관측을 시작하는 경우 센서의 특성 차이로 인하여 각각의 육상관측 위성센서에 적합한 TCT 계수들이 장기 분석을 통하여 수립되어야 한다. GOCI 센서는 해양관측이 주 목적으로 개발되었으나 영상의 상당 부분은 육지를 관측하고 있으며 밴드 구성은 육지관측에도 일반적으로 이용되는 Visible-Near InfraRed(VNIR) 영역의 정보를 포함하고 있다. 또한 GOCI 센서의 높은 시간 해상도는 지표의 일별 변화의 관측에도 유용하게 사용될 수 있다. 이러한 장점을 이용하여 GOCI 센서에 대한 고유한 TCT가 제공된다면 GOCI 센서의 관측범위 내에서 준 실시간으로 지표변화에 대한 분석과 해석이 가능할 것이다. TCT는 일반적으로 "Brightness", "Greenness", "Wetness"의 세 가지 정보를 포함하지만, ShortWave InfraRed(SWIR) 파장대역이 없는 GOCI 센서의 경우에는 "Wetness"의 정보를 얻을 수 없다. GOCI 센서의 높은 시간 해상도의 활용을 극대화하기 위해서는 "Wetness"의 정보가 제공되어야 한다. "Wetness"의 정보를 얻기 위해 GOCI 주성분 분석(Principal Component Analysis: PCA) 공간을 MODIS TCT 공간에 선형 회귀하는 방법이 사용되었다. 이 연구에서 산출된 GOCI TCT 계수는 정지궤도의 특성에 의해 관측 시간대별로 다른 변환계수를 가질 수 있다. 이 차이를 알아보기 위하여 GOCI TCT 자료와 MODIS TCT 자료 사이의 상관관계가 비교되었다. 그 결과, "Brightness"와 "Greenness"는 4시 자료, "Wetness"는 2시 자료의 변환계수가 선택되었다. 최종적으로 산출된 변환계수의 적절성을 평가하기 위하여 GOCI TCT 자료는 MODIS TCT 영상 및 여러 육상 파라미터들과 비교되었다. GOCI TCT 영상은 MODIS TCT 영상보다 지표 피복의 분류가 더 세밀하게 표현되었으며, GOCI TCT 공간의 지표 피복 분포도 유의미한 결과를 보여줬다. 또한 GOCI TCT의 "Brightness", "Greenness", "Wetness" 자료는 Albedo($R^2$ = 0.75), Normalized Difference Vegetation Index(NDVI) ($R^2$ = 0.97), Normalized Difference Moisture Index(NDMI) ($R^2$ = 0.77)와 각각 비교적 높은 상관관계가 나타났다. 이러한 결과들은 적절한 TCT 계수의 산출이 이루어졌다는 것을 보여준다.

• 생태와환경
• /
• 제35권1호통권97호
• /
• pp.28-35
• /
• 2002

이상으로 농업기반공사에서 운영중인 193개소의 농업용 저수지에 있어서 저수지의 물리적 인자가 저수지의 수질에 미치는 영향에 대하여 고찰해 본 결과, 우리나라의 농업용 저수지는 외국의 자연호수에 비하여 유역면적이 상대적으로 크기 때문에 많은 영양염류의 유입가능성이 큰 것으로 나타났으며, 유역면적에 대한 호소 표면적 비율이 비슷한 외국호소와 비교하여 볼 때 체류기간이 작은 것으로 나타났다. 유역면적에 대한 저수지 면적 비율(DA/SA)과 수리학적 체류기간(td)과의 회귀분석 결과 관계식은 $t_d\;=\;42.21(TA/ST)^{-1.0}\;(R^2\;=\;0.89)$로써 역비례의 관계를 나타내었다. 유효저수량을 $t-d$로 나누어 일유출량을 구한 결과 관계식은 $Q_d\;TA/42.21$로써 일 평균 유입량은 총 면적에 선형적인 비례관계를 나타내었으며, Modified Tank 모형에 의한 예측치와 매우 유사한 값을 나타내었다. 연간 표면적당 인 부하율은$4\;gT-Pm^{-2}yr^{-1}$ 이하가 전체 농업용 저수지의 52%를 차지하여 가장 일반적인 연간 표면적당 인 부하율이었으나, DA/SA의 비율이 100 이상의 경우에는 $4\;gT-Pm^{-2}yr^{-1}$이하가 차지하는 비율이 0%, $10\;gT-Pm^{-2}yr^{-1}$ 이상이 차지하는 비율이 100%로서, DA/SA의 비율이 커질수록 연간표면적당 인 부하율은 일반적으로 증가하는 것으로 나타났다. 수심은 Chl-a의 농도에 영향을 주는 물리적인자인데 일반적으로 수심이 깊을수록 Chl-a의 농도는　감소하는 경향을 나타내었다. 이러한 경향은 Chl-a의　농도가 $10\;{\mu}g/l$이하와 $40\;{\mu}g/l$이상에서 뚜렷하게 구분이　되었으며, 일반적으로 우리 나라의 농업용 저수지 계획시 수심을 깊게 할수록 저수지 부영양화 관리에는 유리할 것으로 나타났다. 호소내에서의 총인 농도를 종속변수로 하고 배출부하량과 유효저수량,　저수지 표면적,　유역면적에 대한 저수지 표면적을 독립변수로 하여 회귀분석한 결과 관계식은 $C\;=\;6.03{\times}W^{0.68}{\times}V^{-0.53}{\times}SA^{-0.21}{\times}(DA/SA)^{-0.50}$이며 $R^2$가 61%의 설명력을 나타내어 농업용 저수지의 인 농도 경향을 예측하는데 사용 가능할　것으로 생각된다. 이상의 결과로서 우리 나라 농업용 저수지의 수질예측에 있어는 복잡한 유출모형이나 유역모형, 호소수질모형 등을 사용하기보다는 본 연구에서와　같이 기존자료를 분석하여 도출한 경험적 관계식을 사용하면 비교적 적은 노력으로도 상당히 신뢰성있게 예측하여 수질관리에 도움이 될 수 있을 것으로 판단된다.

• 한국사회복지학
• /
• 제57권3호
• /
• pp.223-249
• /
• 2005

1990년대 후반 경제위기를 경험하며, 우리나라에서 실업의 문제에 대한 사회적 관심과 더불어 실업자의 재취업 가능성에 대한 연구가 지속되고 있다. 재취업과 관련된 기존의 연구들은 인적자본이론을 토대로 연구를 진행해 왔는데, 본 연구는 기존의 연구와는 다르게 Sen에 의해 제시된 인간능력관점을 기초로 하면서 이러한 요인의 실업자 재취업에의 영향력을 인적자본론 요인과 비교론적 차원에서 분석을 시도해 보았다. 아울러, 본 연구는 인간능력관점에서 논하는 사회적 능력 요인 중 핵심적 요인으로 간주되는 성별 차이라는 요인에 주목하면서 성별에 따라 인적자본요인과 인간능력요인이 재취업에 어떠한 영향을 주는가를 비교 검토하였다. 본 연구는 한국노동패널 5차년도 자료 중 연구목적에 적합한 1,871사례를 대상으로 선정하여 재취업의 위험도를 콕스회귀분석을 통해 분석하였다. 분석 결과, 학력, 기술수준 적합도, 근속기간 등의 인적자본이론 요인들은 기대와는 다르게 재취업의 가능성에 부적인 영향을 갖는 것으로 나타난 반면, 기초교육능력, 의료보험가입과 사회보험수혜 및 수도권거주여부 등의 인간능력요인들은 기대한 바에 일치하며 재취업에 정적인 영향을 주는 것으로 나타났다. 성별에 따른 분석결과를 보면, 성별에 따라 각 요인들이 상이하게 재취업에 영향을 주는데 그러한 경향성이 인간능력요인에서 더욱 부각되는 것으로 나타났다. 본 연구의 결과는 재취업에 있어서 인적자본관점에서 제시되는 요인들에 대한 고려도 중요하지만, 동시에 인간능력관점에서 강조하는 사회 제도적 요인이 역시 매우 중요한 역할을 할 수 있음을 시사해 준다. 그리고 성별에 따른 비교연구 결과는 재취업 현상에 대해서 인간능력관점과 더불어 성분절적 노동시장구조에 대한 관점과 함께 이에 대해서 다양한 접근이 필요함을 시사해 준다. 본 연구는 기존 연구에서 실증적으로는 주로 국가단위에서 분석되었던 인간능력관점을 개인단위차원에서 분석하면서 재취업에의 영향력을 다각적인 시각에서 검증하여 인간능력개발과 관련된 사회복지적 개입의 의미를 제고시킨다는 점에서 연구의 의의를 갖는다.

• 한국농림기상학회지
• /
• 제8권3호
• /
• pp.190-198
• /
• 2006

본 연구는 벼의 질소영양 상태가 광 이용효율에 미치는 영향을 검토하여 벼 생육모델을 구축하기 위한 기초자료를 얻기 위해 수행되었다. 1998년부터 2000년까지 3개년에 걸쳐서 일반계 및 통일계 품종을 공시하여 다양한 질소시비 조건에서 건물중, 군락의 흡광량 및 질소영양 상태를 조사하였고, 이들 자료를 이용하여 군락의 광 이용효율 및 광 이용효율과 질소영양 상태와의 관련성을 검토하였다. 벼 군락의 흡광계수(K)는 최고분얼기까지 벼의 생장과 함께 증가하다가, 최고분얼기에서 출수기까지는 0.4정도로 일정하게 유지되고, 출수기 이후부터 다시 급격하게 증가하였다. 출수전 벼 군락의 누적 PAR interception과 건물중과는 직선 회귀관계가 성립하였고, 이 직선 회귀계수를 이용하여 군락의 출수전 평균 광 이용효율(RUE, g/MJ of PAR)을 추정하였는데, 평균 RUE는 질소시비량이 많을수록 높아지는 경향이었다. 벼에서 출수전의 RUE는 질소영양 상태를 나타내는 질소영양지수(NNI, nitrogen nutrition index) 및 비엽 질소농도(SLN, specific leaf nitrogen concentration; $g/m^2$ leaf area)가 증대됨에 따라 최대 RUE에 점근하는 다음과 같은 지수함수로 잘 표현할 수 있었다. $$RUE=3.13\{1-exp(-4.33NNNI+1.26)\}$$ $$RUE=3.17\{1-exp(-1.33SLN+0.04)\}$$ 위의 식은 질소영양 조건에 따른 RUE의 변이를 각각 80%와 75% 정도를 설명할 수 있다. 여기에서 구해진 RUE와 NNI 및 SLN 관계는 출수전 전 기간의 평균적인 관계에 근거한 것으로서, 벼 생육모델의 건물중 추정의 구성모델(component model)로서 활용하는 데는 한계가 있으며, 건물중 추정 구성모델로 이용하기 위해서는 보다 짧은 기간을 대상으로 이와 같은 분석이 이루어져야 할 것으로 판단된다.

• 한국농림기상학회지
• /
• 제11권2호
• /
• pp.72-78
• /
• 2009

관측밀도가 동일한 조건에서 단위격자점의 크기를 줄일 경우 PRISM 방식에 의해 추정된 강수량 분포 가 단위격자점의 크기를 줄이기 전에 비해 개선되는지 확인하기 위해 PRISM 코드를 수정하여 $270m{\times}270m$ 격자점 단위로 구동할 수 있도록 하였다. 남한 전역의 지형자료를 270m DEM으로부터 준비하고 432개 기상청 자동기상관측소의 2007년 월별 적산강수량 자료를 입력자료로 하여 각 격자점의 PRISM 회귀식을 도출하였다. 회귀모형과 DEM 고도에 의해 각 격자점의 월별 적산강수량을 추정한 다음, 추정된 강수량분포도로부터 한국수자원공사 우량관측소 166개소에 해당하는 격자점의 자료를 추출하여 해당관측소의 실측값과 비교하였다. 동일한 강수자료를 이용하되 이번에는 5km 격자점의 PRISM 회귀모형을 유도하여 강수량 분포도를 작성하고 166개 지점 추정강수량을 추출하여 실측자료와의 차이를 RMSE로 표현하였다. 5km 대신 270m 분해능의 DEM을 사용할 경우 월 강수량이 100mm 이상인 경우 평균 10%의 오차 감소효과가 확인되었다.

• 한국농림기상학회지
• /
• 제2권2호
• /
• pp.47-61
• /
• 2000

북한 시ㆍ군별 벼 생육모의결과를 토대로 벼 재배 적합성 여부를 판정하였다. 생육모의에 필요한 시ㆍ군 별 일 기상자료는 지형기후학적 공간내삽기법을 근거로 한 3단계 과정을 통해 생산하였다. 우선 기온의 경우 51개 남북한 표준관측소의 14년간(1981~1994) 월평균값을 관측지점 위도, 해발고도, 해안거리, 경사도, 개방도 등 지리지형변수에 회귀시켜 얻은 통계모형(RMSE=0.4~1.6$^{\circ}C$)을 북한전역에 적용시켜 1 km$\times$1 km수평 격자점 단위로 월별 평균값을 추정하였다. 강수량의 경우 상대적으로 자료가 풍부한 남한의 지형-강수 관계를 도출하여 이를 북한지방에 적용한 윤 (2000)의 방법에 의해 월별 강수량 분포도를 작성하였다. 일사량의 경우 남한 19개 관측소의 14년간(1984~1997) 월 평균 수평면 전천일사량 관측값의 추정식([일사량, MJ m$^{-2}$ day$^{-1}$)=0.344+0.4756[대기외 일사량]+0.0299[남쪽 개방도]-1.307[운량]-0.01[상대습도], 결정계수 0.92, RMS error 0.95)에 의해 북한 지방 27개 지점의 일사량 자료를 복원하였다. 이를 거리역산가중법으로 공간내삽하여 북한전역의 월별 일사량 분포도를 작성하였다. 두 번째 단계에서는 얻어진 1 km$\times$1km 격자점 기후값을 183개 북한 시ㆍ군별로 공간평균값을 취했다. 마지막으로 시ㆍ군 단위 월별 기후값을 이용하여 통계적인 방법 (Pickering et al., 1994)에 의해 30년간의 일별 기상자료를 생성하였다. 북한의 대표적인 벼 품종 생육조사자료를 토대로 CERES-rice 모형의 유전적 모수를 조정하고, 준비된 기상자료를 입력시켜 183개 시ㆍ군별 벼의 생육을 30년치씩 모의하였다. 생육모의결과 중 성숙기와 수량 관련 특성을 점수화 하여 각 시ㆍ군의 벼 재배용 농업기후학적 잠재력을 정량적으로 표현하였다.