DAFTAR PUSTAKA Albrecht, J.C. 1971. A Climate Model of Agriculture Productivity in Missouri River Basin. Publication in Climatology. Laboratory of Climatology, Vo1.24, N0.2. Amien, I., P. Rejekiningmm, A. Pramudia dan E. Susanti. 1996. Effects of Interannual Climate Variability and Climate Change on Rice Yield in Java, Indonesia. Water, Air and Soil Pollution, 92: 29 - 39 Arsyad, S., A. Priyanto dan L. I. Nasoetion. 1985. Pengembangan Daerah Aliran Sungai Lokakarya Pengembangan Program Studi Pengelolaan DAS. FakuItas P-ana IPB, Bogor.
Ayoade, J.O. 1983. Introduction to Climatology for The Tropics. John Wiley and sons. New York.
Bach, W. 1980. Impact of Increasing Atmospheric COz Concentration on Global Climate: Potential Consequeness and Corrective Measure. In: Environment International: Carbon Dioxide and Climate. J. of Science, Techwloay. --Health, Monitoring and Policy. (Ed). A.A. Moghissi. pp. 215 - 228.
. 1989. Growing Consensus and Challenge Regarding a Green House C b a t e . In Climate and Food Security. IRRI-AAAS. pp. 289 - 3 10. Bachelet, D., G.A. King daa J. Chaney. 1995. Climate Change Scenarios. In:: Modeling the Impact of Climate Chroge on Rice Production in Asia. (Ed.%). Matthews, R.B., M.J. Krope I).Backlet dan H.H. van Laat. IRRI.pp. 6782.
Badan Meteomlogi daa Geofisika. t994. Tipe Hujan di Indonesia (1%1 Depaaemenferhubungan.s-.
-
1990).
Bahi Rebabilitasi Lahan clan Koaservasi T d Wilayah IX. 1996. Pola Rehabiiitasi Lahan dan Konservasi Tanah DAS Walanae. BUKU I (Laporan Utama). Ujung PandangBey, A. 1981. Meso Climatonomy Modelling of Four Restoration Stages Foliowing Khkatau's 1883 Destruction (Pumice Ash Desert-Savana-Young -Forest-Rahhrest). PbD Thesis Dept. of Meteorology, Univ. of Wwonsin, Madison. Bey, A. Irsal Las dan Y. Koesmaryono. 1990. Laporan Akhir Pernodelan Evapoklimatonomi untuk Prediksi Lengas Tanah clan Evapotranspirasi Sub DAS Konto @AS Brantas) Jawa Timur Sebagai Studi Kasus. FMIPA IPB, Bogor.
Blrmey, H.F. 1955. C h t e as an Index of Irrigation Needs in Water. The Yearbook of Agriculture. The United States Depart. Of Agriculture, Oxford & IBH.
Boer, R 1 m . Analisis Resiko Kekeringan untuk-P dan Produksi Kcdelai di Flores, Nusa Tenggara T i . Lapom Riset Uoggulan Terpadu IV f 1996- 1998). Kantor Memteri Negara Riset dan Teknologi. Dewan Riset Nasional. Jaksrta 199sb. Keragaan dan Dinamika Sumberdaya IWim di Indonesia Serta Disajikan pada PeLatbn Analisis dan Pemantauan Faktor Iklim untuk Pertanian. Diseknggarakan oleh Direktorat Bina PerIindungan Tunaman, Departemen Pertanian, Cipayung 30 Skptember 1998.
Kaaannya dengan Pertanian.
,Kahpuddia, M.A. Ratag dan A. Bey. 1999. Impact of Doubling C 9 on Forest Productivity. Paper presented at International Symposium, Hotel Millenium, Jakarta. Boer, G.J., K. Arpe, M. Blacburn, M. m u & , W.L. Gates, T.L. Hart,H. LeTreut, E. Roeckner, D.A Sheiain, I. Simmonds, RN.B. Smith, T. Tokioka, R T . Wethaakl and D. Wi&amsom. 1991. An Int~comparisonof the Climates Simhkd by 14 Atmospheric General Chilation Models. WCRP Report No. 58, WMOfID-425, Word Meteor. Organ., Geneva, 156 p.
Chalnaborty, 3.dan M. U.1994. Moonson Climate and Its Change in a Doubled COz Atmosphere as Simulated by CSIRO9 Model. TAO, Voi. 5, No. 4. Taipei Trtiwaa, R0.C. pp. 5 15 - 536. Chow, V.T. (Ed.). 1964. Handbook of Hydrology. A Compendium of Water Resource T e c h ~ o l o ~ yMcGraw-Hill . Book Co. New York. Cooley, K. R 1990. E h t of Ca-induced Climate Change on Snowpack and Streamftow. HydroL Sci. J., 35: 51 1-534. de Rozari, M.Bi. , Koesoebiow, N. Sinukaban, D. Murdiyarso dan K. Matcarim. 1990. Assessment of Socio Economic Impacts of Climate Change in Indonesia in the Potential Socio-Economic Effect of Climate Cbaage in South- East Asia. United Nation Environment Program (IWEP).
Direktorat Jenderd Reboisasi dan Rehabilitasi Lahan. 1991. Rencana Pengelohan DAS Terpadu DAS Saddang. Buku I1 (Laporan Utama). Kerjsama Proyek Pengembangan Rencana Pengelolaan DAS Terpadu dengan Lembaga Sumberdaya Infbrmasi In* Pertanian Bogor. Dooge, J.C.I. 1960. The Routing of Ground Water Recharge Through Typical Elements of Linear Storage. International Association of Scientific Hydrology, No. 52. Doorenbos, J. clan W.O. Pruitt. 1977. Guideline for Predicting Crop Water Requirements. FAO. Rome.
Fds, S.B., J.T. KiehI, A.A. Lacis aad M.D. S c h d o N . 11991. Idhmi Cooling Rate Caldation in O p e m i d Generai C i i n Me Comparison with &nchmark Computation. J. Geophys. Res., 96: 9105-9120. Fradce, O.L. and N.E. McClymonds. 1972. Summary of The Hydrologic Situation On Long Island, N.Y., as guide to Water Management =es, In: Hydrology and Some E&t of Urbanization on Long Island, N.Y. Geological Survey Professional Paper 627 F. Washiagton D.C. Gill, N. 1979. Watershed Development with Special Reference to Soil and Water Co~lservation.FA0 Soil Bun. No. 44.
Godal, O., L. Sygna, J. S. Fuglestvedt rfan T. Berntsen. 2000. Estimate of Future Climate Based on SRES Etnission Scenarios. Working Paper 2000:4 (http://www. cicero.ui0.m). Gordon, G. 1980. System Simulation. Prentice-Hall of India Private Ltd., New Delhi.
Haan, C.T. 1977. Statistical Methods in Hydrology. The Iowa State Udv. Press. Ames, Iowa. 378p. Hack, J.J. 1992. Climate System Simulation: Basic Numerical and Computational Concepts. In: Climate System Modeling (Ed.) Kevin E. Trenbert. pp. 283-318. Js., H. 1985. Perencane~an Pengelolaan D a e d Aliran Sungai, M. Dalam A. Mangundikoro, ed. Prosiding LokaIcarya Pengelolaan DAS Terpadu di Yogyakarta, 3-5 Oktober 1985. Departemen Kehutananl Universitas GajahWh, Yogyzhrta.
Haenunan,
422-43 I.
Hampicke, U. 1980. Source a d Sinks of Carbon -oxide in Terrestrial Ecosystems In: Environmental International: Carbon Dioxide and Climate- J. Science. Tech-nology, Health, Monitoring and Policy. (Ed). A.A. Moghissi. pp:
-
301 316.
Hay, L.E., W. A. Battaglin, R S. Parker dan G. H. Leavesley. 1993. Modeling the E*ts of Climate Change on Water Resources m the Gunnison River Basin, C o b d o . In: Environmr:nta.l Modeling With GIs. (Eds) M. F. Goodchild, B. 0.Parks clan L.T. Steyaert. New York Word. Oxford University Press pp. 173 - 181. Hengeveld, H. and P. Kertland. 1995. An Assessment of New Developments Relevant to the Science of Climate Change. Climate Change Newsletter Vol. 7 (3). m.1-13. Hillel, D. 1977. Computer Simulation of Soil-Water Dynamic. A Compendium of Recent Work. IDRC, Ottawa. Houghton, J.T. 1991. ScientSc Assessment of Climate Change: Summary of the IPCC Working Group I Report Ln: Climate Change: Science, Impacts and Policy. Proceedings of the Second WorId Climate Conference, J. Jager and H.L. Ferguson Wds.), p. 23 - 45. Intergovenm~ntalPanel on Climate Change (IPCC). 1990. Climate Change: The LPCC Impacts Assessment. In: W.J. McG. Tegart, G.W. Sheldon dan D.C. GdEths (Eds.), A u s t d h Government Publishing senice, Canp. (4-1) - (442).
. 1992. Climate Change 1992: The Supplementary report to the IPCC Scietsc Assessment. Ia: J.T. Houghton, B.A. Calbder d m S.K. Varney (Eds.), Cambridge University Press, Cambridge, UK., 200pp. . 1995. Climate Change 1994: Radiative Forcing of Climate Change and An Evahration of the rPCC IS92 Emission Scenarios. In: J.T. Houghton, L. G. M e b Fiio, J. Bruce, Hoesung Lee, 3.A. Calhuder, E. Haites, N. Harris dan K. Maskell (Eds). Cambridge University Press, Cambridge, 339 pp. . 1996.
Climate Change 1995:
The Science of Climate Cbange. Irx Houghton, J.T., L.G. Meiro F i , B.A. Callemlar, N. Harris., A. Kattenburg dan K. Maskell (Eds.). Cambridge University Press. Cambridge, U.K, 572p.
.
1999. Special Report on
Emissions Scenarios (SRES). http://sres.ciesin.ord. I
d Las. 1992. Pewilayahan Komoditi Pertanian Berdasarkan Model Ildh Kabupaten Sikka dan Kabupaten Ende NTT. Disertasi Program Pascasarjana IPB, Bogor. Tidak Dipublikasikan
Je&rs, 3.N.R 1978. An Introduction to System Anatysis with Application. Edward Arnold PubL Ltd., London. Kaimuddin. 1999. PanodeJan Dampak Perubabn IMim dan Tataguraa Lahan DAS Mamasa. Laporan Akhir Penehan Dosen Mu& Lembaga Penelitian UNHAS Makasar. Tidak DipuMikasikan.
..
,H. Pawitan, A. Bey, S. Hardjoarnidjojo, R. Boer and M. A. Ratag, 1999 Compmtive Analysis of General C i o n Madels tbr Indonesia. Paper presemted at Intemationsl Symposhnn on Land-Use Chaage and Forest Management h r Climate Cbenge Mitigation, Bogor, 18* - 20* October 1999. Kachmar%h, Z. and J. Napiorkowski. 1996. Water Resoutces -on Strategy in an Uncertain Envkmmnt. I.:Adapting to Climate Change. Awesunen& and Issues. Joel B. Smith, Neeloo Bhatti, G ~ M & V. M a Ron Beniolat; Max Campoq Bubu Jaltow, Frank Rijsbemran, Mikhail I. Budyko aad RK.Dbton. Sp-, New York. pp. 21 1 224.
m).
-
Kondratyev, K.Ya 1969. Radiation in The atmasphere. Intmmational Geophysics Series Vd. 12. Acad. Press. New York. 912 p.
Laciq kk and J.E. Haaren. 1974. A P a r a m . o n fot the Absurption of Solar Radiation in the Earth's Atmosphere. 1. Atmas. Sci., 3 1: 1 18-133. Laimeheriwa, S. 1994. Analisis Neraca Air MeModel Evapoklimatonomi dm Topog-Yield di Sub DAS Citere-Jawa Barat. Tesis Program Pascasarjana IPB. Tidak Dipubljkasikan.
LA,M., P. H. Whetton, A. B. P i o c k dan B. Chakraborty. 1998. The Greenhouse Gas-Induced Climate Change over the Indian Subcontinent as Projected by General Circulation Model Experiments. TAO, Vol. 9, No. 4, Taipei, Taiwan, R.O.C. pp 673-690.
Lawson, T.L. dan J.W. Pendletoa 1989. Climate variability and sustainability of crop yields in the humid tropics. In: C h t e and Food Security. IRRIAAAS. pp. 57-68.
Lee, R. 1980. Forest Hydrology. Colombia Univ. Press. New York. 4 3 lp. Lernon, E., D. W. Stewart and R. W. Shawcroft. Codield. Science 174, pp. 371-378.
1971. The Sun Work in a
Lettau,
H. 1969. Evapotranspiration Climatowmy I: A New Approach to Numerical Prediction of Monthly Evapotranspiration, Runoff and Soil Moistme Storage. Mon. Wea. Rev., Vol. 97, No. 10 pp. 691-699.
. 1977. Climatonomical Modelling of Temperature Response to Dust Contamination of Antartic Snow Surface. Boundary-Layer Meteorology, 12. pp. 213-229.
. and M. W. Baradas. 1973. Evapotranspiration Climatonomy II: Refinement of Parameterization, ExcepMed by Application to The Mabacan River Watershed Mon. WeaRev., VoL 101, No. 8, pp. 636-649. ., K Lettau and D.C.B. Molioa 1979. Amazonia's Hydrologic Cycie and The Role of Atmospheric Recycling in Assessing Deforestation Effect. Mon W e a Rev., Vol. 107, No. 3, pp. 227-238. Lieth, H. and E. Box. 1972. Evapotranspiration and l ? ? y Productivity. Paper on Thomthwaite Memorial Vol. 2, J.R. Selected Topic in Climatology, Mather ed., pp.3746. Linsley, RK. 1951. The Hydrobgic Cycie and Its Relation to Meteorology River Forecasting. Compendium of Meteorology, pp. 1048-1054. Manetsch, T.J. snd G.L. Park. 1977. System Analysis and Sixnufation with Application to Economic and Social Systems, Part I. Department of Electrid Engbxrkg and system Sciences, Michigan State University, Michigan.
Manabe, S. and R.T. Wetherald. 1987. Large-Scale Changes of Soil Wetness Induced by an Increase in Atmospheric Carbon Dioxide. J. Atmos. Sci, 44: 121 1-1235. Mather, J.R. 1978. The Climate Water Budget m Environmental Analysis Lexington Books, Toronto. McCabe, G.J. dan Wobck, D.M. 1992. Sensitivity of irrigation demand in humid temperate region to hypotetical climate change. J. Hydrol., 136: 353 - 380. McGregor, J.L., H.B. Gordon, I.G. Watterson, M.R Dix and L.D. Rotstayn. 1993. The CSIRO 9-level Atmospheric General Circulation Model. CSIRO Div. Atmos. Res. Technicai Paper N0.2 6 pp. 1-89 McKenney, M.S. d m Rosenberg, N.J. 1993. Sensitivity of Some Potential Evapotranspiration Estimation Methods to Climate Change. Agric.For.Meteorol., 6481-110.
M o u , M.A., Hadjisawa, P.S., Kouvopoulous, Y.S. dan Ahteos, H. 1991. Regional Climate Change Impacts: 11. Impact on Water Management Works. RydroL Sci. J., 36: 259 270.
-
hiiize, J.H. and J.G. Cox. 1963. Essential of Simulation. PrenticehiU, Inc. Englewood CNew Jersey. 234 p. Molion, L.C.B. 1975. A Climatonomic Study of The Energy and Moisture Fluxes of Amamne Basin with Consedaration of Defbrestation Effects. Ph.D Thesis, Dept. Of Meteorology, Univ. of Wisconsin, Medison. 133 p. More, RJ. 1969. The Basin Hydrological Cycle. In: Introduction to Physical Hydrology. (Ed.) Richards, J. Chorley. Methuen & Co. Ltd., Bungai Sufnok. pp. 245-264. Murray, D. L. and A. M.H. Gorgens. 1981. Storm Runoff Analysis on Three Semi-Arid Catchment. Water of South Mca. Vol. 7. No. 4. Nasendi, B.D. 1979. Peranan Analisis Sistem dalam Penelitian Sistern Hidrologi D M . Kebutanan Indonesia, No. 8 Th.V. pp. 35 - 37. Nieuwok, S. 1977. Tropical Climatology. An Introduction to the Climates of the Low Latitudes. John Wiley & Sons, New York. Oguntoyinbo, J. 1970. Deflection Coefficient of Natural Vegetation Crop and Urban Surf;ace inNigeria. Quart. J. Roy. Met. Soc. %: 430 - 441.
Parry, M., M.BL Rozari, A. L. Chong dan S. Panick 1992. Socio-economic Impacts of climate change in South-East A s k UNEP-Earthwatch: NairobA, Kenya Pawitan, H., J.S. Baharsjah, R. Boer, Le I. Amein dan B.D. Dasanto. 1996. Keseirnbangan Air Hidrologi Whyah Indonesia Menkabupaten. Laporan Penelitian. Lembaga Penelitian Institut Pertanian Bogor. Pittock, AB. 1993. Climate d o development. In: Modelling Change in Environmental Systems. (Eds). A.J. Jakeman, M.B. Beck clan M.J. McAleer. pp. 481 503.
-
Rahayu, S.R. 1993. Pendugaan Evapotranspirasi d m Lengas Tanah Daerah Aliran Sungai @AS) C i l i m g Huh dengan Model Evapoklimatonorni. Laporan Masalah Khusus, Jurusan Geomet FhiZrPA IPB, Bogor. Tidak Dipublikasikan. Rao, AR., J.W. Dellew and P.B.A. Sarma. 1972. Conceptual Hydrologic Models fbr Urbankhg B a s h . J. Hydr. Div., A. Soc. Civ. Engrs. pp.1205-1220. Ratag, M.A., Asiati, S., Nuryanto, A, Siswanto, B and Suyadhi 1998. Greenhouse gases, climate change and El-Nino. In Economic Assessment of Greenhouse Gas Abatement Options. ALGAS-FINAL National Workshop. Robock, A. 1989. Global climate change. In: The Potential of Global Climate Ctaange on The United States, Smith J.B. clan Tirpak D. (Eds)., U.S. Environmental Protection Agency, Washinton, DC, pp. 9 28.
-
Rosenkg, N.J. 1989. Potensial effect on crop production of carbon dioxide enrichment m the atmosphere and greenhouse induced climate change. In: Climate and Food Security. IRRI-AAAS. pp. 359 374.
-
Rutter, A.J. 1981. Concluding Remarks. In: Plant and Their Atmospheric Environment. (Eds.) Groce J.0, E.D. Ford and P.G. Jarvis. The 21mSymp. of British Ecological Soc., Blackwell Sci. Public., Odord. Schwab, G.O., RK. Frevert, T.W. Edmister and K.K. Barnes. 1981. Soil and Water Conservation Engineering. Third Ed. John Wdey and Sons, New York. Schlesinger, M.E. 1986. Equiliium and Transient Climate Warming Induced by Increased Atmospheric C a . Clhnate Dynamics, 1: 35-5 1. Sch-opE
M.D. and S.B. Fek. 1985. Improvements to the Algoritm fbr Computing COz Tmmmkivities and Cooling Rates. J. Geophys. Res., 90: 10541-10550.
. 1991. The Simplified Exchange Method Revisited: An Accurate, Rapid Method for Computation of Inhued CooRate and Fluxes. J. Geophys. Res., 96: 9075-9096. Sellers, W. D. 1972. Physical Climatology. Univ. of Chicago Press, Cbicago. Seyban, E. 1977. Dasar-Dasar Hidrologi. Edisi II (Terjemahan: S. Sub~igyo1983). Gajah Mada University Press, Yogyakarta. 388 p. Sitorus, RP. 1985. E v b i Sui&erdaya Lahan. Penerst PT. Tarsito Bandung.185 p.
Skelly, W.C., A. Henderson-Sellers d m A.J. Piiman. 1993. Land suface data: Global Climate Modeling Requirements. In: Environmental Modeling with GIS. M.F. GoodchiId, B. 0. Parks dan L.T. Steyaert. (Eds). Oxford University Press. pp. 135 141.
-
Smith, J.B. 1989. Methodology. In: The Potential of Global Climate Change on The United States, Smith J.B. dan Tirpak D. (Eds)., U.S. Environmental Protection Agency, Washinton, DC,pp. 9 28.
-
Soeriamgara, I. 1978. Pengelolaan Sumberdaya Alam, Bagian 11. Sekolah Pascasarjaua InstiN Pertaaian Bogor, Bogor. Sosrodarsono, S. dan K. Takeda. Par& Jakarta.
1980. Hidrologi untuk Pengakan. PT Pradnya
Suhatsow, H. 1994. I k h Tropika Dalam: Klimatologi Dasar. Landasan PemaCiran Fisika Atmosfk dan Unsur-Unsur Ikhn. (Ed.) Handoko. Pustaka Jaya. XI: 175-1992. Supomtana, S. 1971. Climatowmy of Bangkok. Dept. of Meteorology, Univ. of Wisconsia, Medison 49 p. Suppiph, R 1994. The Asian Monsoons: Simulation Erom Four GCMs and Likely Changes under Eahsnced Gxedmwe C o d i t i o ~ In: ~ . Climate fmpact Methods fbr Asia aad t k Pacific. (Eds.) A. 3. Jakeman and k B. Pittack. AID-, C a t n h
Tadjang, H.L., Kaimuddin dan A. Yassi. 1990. IClimatologi Pertanian. Bahan Pengajaran. Jurusan Agronomi, FakuItas Pertanian dan Kehutanan UNI-IAS. Ujung Pandang. Tidak Dipublhsikan.
Thornthwaite, C.W. and J.R Mather. 1957. Instructions and Tables for Computing Potential Evapotmmpkation and Water Balance. PubL In Climatology, Drexel Inst. Of Tech. Lab. of Climatology, 10, 185-311 . United Nation Environment Progranrme. 1992. The Impact of Climate Ghnlnge. UNEPIGEMS Environment Library No. 10. US-CS. 1994. Guidance for Vulnerability and Adaptation Assessments Version 1.O. US Country Studies Management Team. Washington D.C. Wanggai, F. 1975. Kemungkhan Penggunaan Amilkis Sistem dalam Penelitian Hidrobgi DAS Citanduy. Tesis Sarjaaa Kehutanan IPB. Tidak Dipublikasikan.
Watterson, I.G. 1993. Global Climate Modelling. In: Modelling Change in Environmental System. (Eds.) A.J. Jakeman, M.B. Beck and M.J. McAieer. John Wiley & Sons Ltd. Whetton, P.H., P.J. Rayner, A.B. Pittock d m M.R. Haylock. 1994. An -at of Posst'Me Climate Change in the Austdan Region Based on an Intercomparison of Generai C i a t i o n Mod* Results. Amxican Meteomlogical Society. pp. 44 1 462.
-
Wiersum, K.F. 1979. Introduction to Principles o f Forest Hydrology and Erosion with Special Reference to Indonesia Lembaga Ekologi Universitas Padjadjaran Bandung. Wilson, C . k and Mitchell, J.F.B. 1987. Simulated Climate and CO2-Induced Climate Change Over Western Europe. Climate Change, 10: 1142.
Tabel L a m p i 1. Diskripsi u m u m dan run model GClMs yang digunakan dalam penelitianini.
1. R u n Model Goddard Institute for Space Science (GISS) Tiga model iklim yang d ik-
pada GISS yang telah dirurming dan
digunakan sebagai penilai dampak oleh U.S. Couutry Studies Program, temmsuk model trausien. Run model pertama kali selesai tahun 1982 dan mas& relatif kompetitif (Hansen et al., 1983) dalam (US, CSP, 1994). &US
Model ini mempuuyai
diurnal dengan resolusi horizontal yang rendah (7.83" x
sembilan level vertikal.
10.OO) dengan
Kelompok GISS menggunakan prosedur Q-fluks mtuk
memperbaiki simulasi terhadap laut dan eslsalju. Suhu air laut dan penutupau es dihitMg bembarkan pertukamn energi dengan a t m o s k dari jam ke jam, pengangkutan tertentu b a b q laut (Q fluks), d m kapasitas bahang pencampuran lapisan laut. Model ini juga menggunakan yang disebut "slab-ocean" yang sederhana,
tidak lebih dari 65 m kedalam untuk merun 1 x COZdan 2 x C01, dengan beberapa keragaman keddamcm pencampwan siklus musiman. Kepekaan model adalah 4.2- C (yaitu pemamsm suhu permukaan global 4.2- C
untuk run 2 x COz d
i
i
dengan I x C a ) . Model kedua adalah simuIasi transien selesai tahun 1985. Kasus transien
mirip dengan kasus 1 x COz
dan 2 x COz, tetapi mempunyai kedalaman
pencampurau yang lebih dalam. Hasil simulasi transien tidak ditujukm untuk
digunakan didalam Program Studi Wilayah.
Model ketiga adalah versi baru dari run 1 x C& dan 2 x C a . Model ini mempunyai resohrsi yang terbaik (3.9" x 5.0") dan Q fluks ''slab ocean". Pennulaan
run dengan panjang data selama 70 tahun dengan hanya sembilan level vertikal yang selesai kira-kira Pebruari 1994 (kontrol run clan peningkatan gas sebesar 1 %). Sesudah 70 tahun suhu meningkat sebesar 1.5OC. Luaran model mirip dengan model yang lain, tetapi model ini mmpunyai runoff sungai kedalrun laut. Sesudah 70
tahun penkgkatan muka laut 11 cm; kira-kira 7 c m -tan termal. Penggunaan model ocean lebih maju &
ini kruena ekspansi
i
beberapa run model
sebelumnya Waktu yang dibutuhkan untuk komputasi yaitu dua bulan pada IBM 590 Workstasiun dengan simulasi selama 50 tahun (sepertiga waktu tersebut dig
'
untuk model laut) Model GISS menggumkm skema radiasi "k-distribution" berdasarkan Lack
dan Hansen (1974). 2. Run Model Geophysical Fluid Dynamics Laboratory (GFDL) GFDL mempunyai empat model nm. Versi pertama, menggunakan model
"shb ocean" kedalaman 64 m dan mempunyai resolusi 4.44" x 7-50" (R-15, model s p e w dengan 15 gelombang) (Manabe dan Wetherald, 1987). Run selesai tahun 1985 dan tidak mempunyai aliran atau arus laut atau fluks bahang-palsu untuk
memperbaiki suhu air untuk ikiim saat hi. Model run yang Lain pada GFDL teLah selesai Pebruari 1988. Model ini kira-kira sama dengan sebelumuya, kecuali model ini memasuWcan prosedur Q fluks untuk menghasiUcan kembali suhu permukaan laut
dan luas laut es dari iklim saat ini. GFDL menyebut model ini model "Q-Flux". kepekaan kedua model ini ad-
4.0" C (perubahau global dari lXCO2 ke 2XC02).
Run ketiga dilaksanakan Bulan Mei 1989 dengan sebagian besar model yang sama dengan resolusi yang t e r W yaitu 2.22O x 3.75O (R-30,model spektral dengan
30 gelombang). Belum mempunyai siklus d i d Awal simulasi melalui hasil dari run R-15 dan proses simulasi 10 tahun lebih (untuk 1 X C a dan 2XC02). Dengan resolusi lebih baik, model presipitasi lebih realistik saat d saat
ini.
i
i dengan iklim
Kadar air tanah juga lebih reaIistik, kernungkkmnya karena presipitasi
lebih baik. Kepekaau run R-30 sampai 4.0" C.Run ini memberkin perubahan tidak
realisUc didalam presip-
dari satu grid ke grid bdcutnya dan disini akan
dhduskm oleh penggma.
GFDL sedang meqpahzhn uutuk merun dengan 14 level vertikal, resolusi
R-30 dan mequnyai Q-flux "slab ocean", seperti sebelumuya. Siklus diurnal tidak dirun dalam versi ini
Adanya perubdmn dalam ketinggian awan dan albedo es,
kepekaan suhu lebih kecil dari hasil run sebelumnya yakni sebesar 4.0° C. GFDL transien dengan run 100 tahun selesai tahun 1990. Model atmosfer rnempunyai resolusi latitude-longitude 4.4O x 7.5O (R-15).
Model hut sama
resolusinya dengau latitude tetapi lebih baik dengan longitude (3.75O) yang
membantu menyiapan simulasi yaug lebih baik, dengan Iapisan perbatas diilahan bumi b a g k barat untuk kondisi iklim saat ini (current cIimate) pads kedalam laut 5
km.
Dinamika laut yang digudam pada model ini untuk kondisi ikIim saat ini
(current climate) Iebih
baik dari pada model "slab ocean" yang digunakan model
paling akbir. Model atmosfer secara esensial sama dengan model GFDL yang mempunyai kepekaan 4.0° untuk 2XCOz. Run yang paling akhir dilakukan, yaitu
pada tahun 1989 yang digunakan didalam studi IPCC. Tahapau waktu (time step) yang d i -
dalam mode1 laut menyebabkan beberapa masalah kecil. Mode1 di
run kembab tahun 1990 dengan sedikit lebih pendek tahapan waktu yang dig& dalam model laut. Model GFDL menggunakan skerna Mmabe-Stone untuk parameterisasi radiasi gelombang panjang (lihat Fels et al., 1991) dan skema Lacis dan Hansen
(1974) untuk radiasi gelombang pendek. 3. Run Model Canadian Climate Centre Model (CCCM)
CCCM tdah dikembangkan sebagai model dan mempunyai run model untuk kondisi 1 X C a dan 2XCOz
(Boer et al., 1991).Model ini mempunyai resolusi
tertinggi (3.75" x 3.7S0) d
i dengan kebanyakan model simuIasi iklim
dengan 10 level vertjkal.
Modei ini mempunyai kedalam 50 m "'slab ocean" dan
menggunakan prosedur Q-Flux.
Model ini mempunyai skema tertentu untuk
kelembaban udara untuk pengembangan awan dan lebih rendah &'banding -tan
beberapa
&US
diurnaL Kepekaan iklim
modei GCMs b y a ; Model menunjukkan
suhu sebesar 3.5O C dan peningkatan presipitasi sebesar 3.8 % saat
komentmsi 2XCOz. Data dilepas untuk umum bulan Mei 1991. 4. Run Model U.K. Meteorological Office (UKMO)
UKMO mempunyai dua model run. Pertama, run lXC02 dan 2XCOz t e U selesai bulan Juni 1986 ( W i n dan Mitchell 1987). Model ini mempuuyai resolusi
5.0° x 7.5O, 1 1 level vertikal @ads permukaan sigma) dan siklus diurnal. mempuqai "slab ocean'' dengan kedalaman 50
Model
m, seperti model yang lain Prosedur
Q-Fluks dengan perbaikm simuhsi pada suhu air hut dan es laut saat ini. Model menggmakan "spin-up" kira-kira 20 tahun clan meaguiqulkan data 15 &hun. Run
model menunjukkau peningkatsn sebesar 5.2" C untuk 2XC02 dan merupalan model yang paling tinggi penkgkatarmya. Urmmmya model GCMs yang lain hanya 4" C. Sectma global model ini mempunyai kepekaan terthggi terhadap kenaikan suhu,
sebab model ini mempunyai lapisan lebih banyak. Secara mum, pemanasan pada
tingkat ini dipercaya cukup tinggi.
Run kedua model UKMO telah selesai bulan Nopember 1989 dengan resolusi 2.50° x 3.75O. Run model t kumulus.
e
d hdungan air pada awan dan penetrasi awan
Model ini lebih realistik menggambarkan pen&ambura.u awan sirrus.
Mode1 ini masih mempunyai siklus d i d dan Q-Flux "slab ocean". kepekaan untuk 2 X C a adalah rnemikhn suhu global sebesar 3.5O C clan me-t
untuk presipitasi
menjadi 9 %
Lampiran 2. Prosedur Menjalankan Model CSIR09-R21
Model CSIR09 dirun dengan menggunakan SGI Power Challenge fbr UNXX. Prosedur menjalankan sebagai berikut: 1) Buat direktori clan subdirektori untuk mengecek jalannya model clan menyimpan
data luaran hasil skenariosasi. 2) Pada subdjrektori buat Be-file data, file batch, file log, file restrat, file input, file history dan file luaranbasil skenariosasi 3) Copy semua file GCM XSIR09 pada direktori kerja dengan w a : Cp /Widbaddclimod/csiro/r2 l/sdkode/* . 4) Lakukan kompilasi file-f2e surnber (subrutin-subrutin), yang s e b e l m y a telah
. .
d m d d a n sebuah jjle yang bemengkompilasi file-file sumber. Untuk melakukan kompilasi file sumber GCM CSlRO9 pada Silicon Graphics dengan rnulti prosessor, yaitu: make-P SGI. File sumber berextension *.f dan hasil kompitasi berextension *.o dan luaran kornpilasi disebut model y a q akan kita guuakan untuk melakukan skenariosasi
-
5) Periksa , apakah fie-file data yang diperlukan sudah tersedia Terdapat sejumlah file data yang sehlu diperlukan untuk suatu resolusi tertentu yaitu:
data albedo (albnew;?l) r data ozon (03-data.new9) data karbondioksida (co2-data.new9) data tekanan perm(psrk21f.dat) data t o p o d (newtop21) r data suhu ~ u k a air a laut ~ Cclim2f sst) r data unG NSIB scheme sib*idat aika NSIB = True(T)) data untuk sea-ice scheme: ocuv*.dat (arus laut), bbm*.dat ( k -
.
lapisan
-1.
Data-data tersebut di atas merupakan data klimatobgis (rata-rata selama 30 tahun). File data-datatersebut dapat diperoleh dengan mencopy sebagai berikut: cp /iklim/banks/clhnod/csiro/r2l/src/dinp/*. 6) Buka file input yang berisi informasi a d mtuk menjalankan model dart luaran stat&& yang kita inginkan dari model tersebut.
7) Sebehun rnenjalankan model, ada beberapa bsl yang perlu diperhatikan pada file input sebagai berikut: Apakah model yang akan dijalankan hanya model atmosfer saja atau model hutan atau kopel keduanya.
-
3erapa lama model akan dijalankan, pilihan yang tersedia: harian (kurang dari satu Mi,bulanan (satu sampai 31 hari), tahunan (satu sampai 12 bulan) dan dapat dijalankan selama beberapa tahun. Periksa, apakah file data yang tertulis sudah sesuai dengan file data yang tersedia dahm resohi yang d i p f i (default resolusi: 64 X 56). File restart yang digumkm harus sesuai dengan walctu atau periode yang
Jika model dijabkan untuk pertama kalinya, maka harus rmmdapatkan data qflux terlebih dahulu dengan cara mengubah flag di peubh saveqhx menjadi T (true). Sesudah data q-flux diperoleh, flag flux untuk peubah saveqqflux dikemenjadi F (&Is) dan untuk peubah qflux diubah menjadi T. Berikan extension pada file restart clan luaran model (file statistik). 8)
Tidak semua luaran model (file statistik) dijalankan, mengingat terbatasnya kapasitas hardisk. Untuk memilih luran model yang diinpinkan. maka pilihlah luarrtn model pada list input &statvars dengan merubah F menjadi T.
9) Setelah** file* input diperiksa dan dilakukan perubahan-perubahan sesuai prosedur yang dmgmkaq maka selanjutnya dapat menjalankan model dengan berbagai cara: Climatolos~icalRun. Pada waktu menjalankan model perhma kaljnya, model perlu penyesuaian dari eikk keadaan awal yang berbentuk file restart. File ini digunakan untuk menharga-harga dari parameter-parameter yang digunakan dahn model dan dirmnakan untuk melaniutkan suatu perialanan model (numing). Dalam menjaiankan model pert* &ya, data suhu permukaau alaut (SST) y m g digunakan adalah data SST klimatologi selama satu tahun. data SST klimatologi ini merupakan data SST rata-rata selama 30 tahun. Model dipaksa menggunakan data satu tahun SST yang sama untuk betahun, sehingga model tersebut mencapai suatu keadaan iklim ratarata. Dengan demikian perjalanan model ini disebut c ~ t o l o g i c a run. l Faktor q-flux dihitung dari masukan fluks panas dan data SST klimatologi. Untuk menghasilkau data q f h q flag pada saveqflux diubah menjadi T dan flag pada qflux d i i menjadi F. Setiap satu bulan data qflux disimpan sebuah file dan setelah climatological run selesai data q£hm yang d i h s j k m dirata-ratakan dan disimpan dalam satu file qflux.dat. Control Run. Setelah model mencapai keadaan iklirn rata-rata, model siap digunakan untuk berbagai simulasi perubahan iklim global. Tetapi sehlumuya model perlu dijalankan selama beberapa tahun dengan memakai data SST yang dihitung oleh model itu sendiri. karena model d i i k a n mengatur j w y a sen&-& maka disebut control run. Dalam control run huuan yang d h d k m akan bervariasi dari tahun ketahun sampai model mencapai keadaan stabil. Keadaan stabil dahm control run tidak berarti bahwa hasil hwwmya mempunyai variasi yang sarna dari tahun ke tahuu, tetapi variasi tersebut tidak
-
begitu &sat &banding waMu sebelum mencapai keadaan stabil. Hasil keluaran dari control run ini diguna kansebagai War pembanding untuk berbagai simulasi perubahm iklim global yang dhkukan ssudahuya. Anomali Run. Jika model &%arkan berjalan dengan mengguaakan data SST hasil perhhmgan sendiri, model menyerupai variasi perubahau iklim yang berubah dari waktu ke waktu. Dengan mengubah suatu masukan, misaLnya kenaikan konsentrasi C& di atmosfer atau kenaikan radiasi sinar matahari maka model dapat digunaksn untuk simuIasi perubahan iklim global yang diselxbkan oleh efkk gas rumah kaca ataupun aktivitas mntahari. S i yang dilaksanakan pada penelitian ini adalah perubahan iklim glow dengan s k d o 1 X (2% clan 2 X C a . Untuk me- simperubaban iklim dapat dilakukan &ogan dua cara yaitu: perubahan transien dan perubahan Iangsung. Pada penelitian ini cara yang diIakukan adalah cara perubabn langsung yaitu data 1 X C@ yang sekarang (current climate) hgsung diubab rnenjadi 2 X C& (peubah iklim pada saat 2 X COz).
Lampiran 3. Prosedur down s c a l k g resohisi spasil data observasi dan modei GCMs menggunakan Winsurfer V.5.0 dengan metode Kriging
1) Buat data spasial Bujur dan Lmtang (posistif untuk lintang utara dan negatif untuk lintang selatan) dengan Microsoft Excel. Data pada lembat worksheet dari microsoft exceI terdiri dark Kolom A = 'Bujur", Kolom B = '%intang'' dan Kolom C clan seterusnya = "Nilai Z" (nilai Z adalah data peubah yang akan di down stalling, yaitu curah hujan, suhu, radiasi dan lain-lain. baik harian, bulanan mupun tahunan). 2) Data yang tetah & i t dalam worksheet di atas, disimpan (save as) dengan extension .csv agar dapat dihca okh wbu&w tersebut. 3) Setelah data shp, maka aplikasi wimsurface d i i A p h i whsmfizr tterdiri dari menu utama, yaitu: File Edit Viw Draw Arrange Grid Map Window Help. 4)
Buat lembar kerja bmu dengan mengKWc File pada menu utama, sehjutnya sub menu New. kemudian klik Plot lalu enter. seDerti t e m t oada Garnbar
Gambon 1. Menu utama whsdkr. 4) Untuk analisis komparatif antara data observasi dengan data model GCMs, lakukan seperti langkah No. 1 dm 2. Untuk memudahkan analhis komparatif, maka &wan grid (resolusi horizontal la. x ion) data observasi yang menyesuaikan dengan model GCMs. Sebagai contoh model CS-9 mempunyai resolusi horizontal (lat x bn) 3.18" x 5.62O, sehingga data observasi (mhhya
suhu) harus diinterpolasi dengan ukuran grid yang sama dengan model CSIR09
dengan metode 'Kriging". 5) Langkah selanjutnya adalah Klik menu Grid pada menu utama, selanjutnya Klik Data pada sub menu @hat Gambar I), maka akan ada perintab untuk membuka nama file data yang telah disinP,m dengan format .csv. Kemudian pilih file yang akan diolah, sebagai contoh nama file d:\tizaldy\gcm\suhuUob*.csv dan Klik OK. Selanjutnya akan muncul tampiIan seperti pada Garnbar 2 dan B e d:\rizaldy\gcm\suhuUob*.csv berubah secafa otomatis menjadi d:kkaMy\gcm\suhuUob*.grd.
G a d x u 2. Menu scattered data interpohtion.
Gambar 3. Menu kriging option untuk memntukan tipe variogram atau semivarioj3rarn.
6) Surfk akan meminta batas-batas minimum dan maksimum bujur lintang serta skala grid (spacing) sesuai model GCMs yang dipiW Sebelunmya harus diisi &ata Columns Z Cpilihan Jan sampai Des). Kemudian klik metode interpolasi yang akan digunakan pada Gridding Method. Surfer menyedkkan beberapa metode interpoksi, diantaranya metode "Krigiig". Catatan: metode Kriging hanya dig& untuk interpolasi data, sehingga batas minimum dan m a k s b bujur htang harus sesuai dengan batas minimum dan maksimum data observasi, namum spacinnya tetap sesuai dengan model GCMs yang dipilih mtuk dikompamMan
7 ) Pada Gambar 2 terdapat pilihaa Qption. Klik a t i o n dan muncul tampilan seperti gambar 3. Yang harus diisi pada Gambar 3 adalah pada Variogram Model. Pilihan yang terdapat pada Type Variogram Model adalah: Gaussian, Hole E f k t , Linier @epault), Quadratic, Ftational Quadratic clan Sperical. Type variognun dapat ditentukan melahsi analisis &gram atau sembmiogram (lihat macro sxohrbgram dengan program Midab). Selanjutnya klik OK d m tampitan kembali ke Gambar 2. IUik OK, maka interpohi data selesai dan data disimpan dalam file dengan extension *.GRD.
7) Agar data yang telah tersimpan dengan extension *.GRD dapat &ha oleh Microsof Excel atau spreadsheet himya, maka data tersebut hams dikonversi dengan mengKLik Grid pada menu utama (Gaml#u l), kemudian Klik convert sehingga m m c d tampilan Opes Grld, M u Klik OK, sehjutnya akan muncul tamp* Save Grid seperti tampilan pada Gaml#rr 4. Pada Gambasr 4 pilih Save File as lalu Kfik ASCII XYZ (*.DAT), kemudian isi Fie Eame sesuai dengan nama yang diinginkan kemudian Ktik OK, sehingga data telab tershpau
w,
Gambar.4. Menu Save Grid untuk menconversi dari data *.GRD ke *.DAT. dengan extension *.DAT dan dapat d i i dengan Microsof Excel. Selanjutnya data
observasi yang telah diinterpolasi dengan metode kriging dapat d (analisis komparatit) dengan model GCMs yang dipilih.
i
i
i
9) Untuk membuat peta Isohyet, Isotherm, Isobar dan lain-lain, maka diperlukan peta dasar, mkahya peta Indonesia atau peta Pulau J a n Peta dasar dapat dibuat dengan berbagai l~lacamsoftware a y a ArcMo atau AutoCAD. Jika peta dasar diiuai dengan AutoCAD dengan extension *.Dm, maka file tersebut dapat
diinport dengan mengKEk Map pada menu utama kemudian KLik Load Base Map pada sub menu, sehingga numcui t a m p h Import File seperti pada Garnbar 5. File-file yang dapat diiaport adalab dengan extension *.GSB (Golden SoBoundary), *.BNA (Atlas Boundry), *.LGO (USGS D i g i t a l Line Graph), * . D m (AutoCAD DXF Drawing), *.PLT (Golden SofwarePlot Call), *.BLN (Golden Safware Blanking), *.CLP (Wimdows Metafile (Picture)), *.TIF (Tugged Image (TIFF)) dan *.BMP (Wimdows Bitmap). Setelah memilih extension yang t e d dan mewlis nama file, Wc OK, maka file tersebut telab selesai diinfort dan dapat dishpan dengan p e r h h s Save As dengan extension *. SRF.
Gambglr 5. Menu untuk Inport file.
Gambar 6. Peta dasar Indonesia. 10)Untuk membuat peta Isohyet, pertama-tama buka file (Open File) peta dasar dengau mengKKk w n pada menu Eile, kemudian pilih nama file yang berextension *.SRF. Misalnya peta dasar Indonesia seperti terlibat pada Gambar 6. 11) Misalnya kiia akan membuat peta Isotherm model CSIR09 dengan skenario 1 X COz, dimana data file tersebut t e a dishpan (Save) dengan extension *.csv Kemudian lakukan langkah s e p r t i No. 5 sampai dengan No. 7. Batas minimum dan ' bujur lintmg pada Gambar 2, harus sarna dengan batas minimum dan maksimum peta dasat. Pada akhir tahap ini (No. 7),data sudah dalam format *.GRD dan siap diolah menjadi konturflsohyet, dengan cara memenu utama Map kemudian diKlil sub menu Contour dan rnuncul perhtab untuk memilih Be yang akan diolah menjadi konturAsohyet dalam format *.GRD seperti terlihst pada Gambar 7. Setelah dipilih tile rsng dikeheddci akan mmcul periutah-perintah tain seperti t e r W pada Gambar 8.
Gambar 7. Menu untuk membuka file yaug telah dishpan dengan extension * .GRD
12) Pa&
gambar 8. Box Contour Level meminta batas maksimum clan minimum
dari data nilai (suhu udara). Interval disini merugilkaa interval dari kontur. Setelah semuanya selesai Klik tomb01 OK clan a h terbentuk kontur/Isotherm, seperti terlihat pada Gambar 9.
13) Setelah peta dasar dan peta kontur/Isotherm berada dalam satu lembar kerja, kemudian keduanya &%lok dengau cara WE dua kali peta dasar. Setelah itu pilih gap pa& menu utama dan Klik Q d y Maps, tnen&asi&m peta konturflsotherm (sebagai wntoh peta suhu udara bulan Januari model CSIR09 I X C a ) , seperti terlihat pada G a d m 10.
sews
-
-
- - - -- -
-
- --
Gambar 8. Menu Contour Map mtuk maenkbn interval kontur.
Gambar 9. Peta dasar (Indonesia) dan peta konturflsothenn yang akan dioverlay.
Gambar 10. H a d overlay antpeta dasar (Indonesia) dengan peta konturAsotherm modeL CSIR09 1 X COZ.
Program malcro minitab untuk menetukan nilai variograrn ''w-mtb'
Note'program menentukart nilaivariogram' Note'letak buju pada kolom c1' Note'letak lintang pada kolom c2' Note'data pada kobm c3' set cloo %
ewl Copy ctOO kl erase cloo Let kl=kl-1 Let k2=1 Let W=1 Let k4=2 Let k5=kl
Exec'c:\data\jarak l .mtbIc1 End 'jarakl.mtb' exec 'a:\data\jataKZ.mtb'k5 let k3=k3+1 let k4=k3+1 let W=k5- 1 end
"jarak2.mtb1 let C 4 ( k z ) = s n r t ( ( ( c l W ) - c l ( k 4 ) ) * * 2 ) + < ( ~ 2 C k 3 * * #jarak
let c5(k2)=(~3@3)-~3(k4))**2 let k2=k2+1 let k4=k4+1 end
#nilai variograrn
Lampiran 4. Listing Program V i Basic Tahap K a l i i (I)
Dim I, K As Integer Dim np, PN, j, KAI As Double Dim P(13), G(13), N(13). a(13), q13), pq13), e1(13), nl(13) As Double Dim ve(13), vn(13), p1(13), mi(13), m(13), y(13), V(13) As Double Dim pb, gb, nb, veb, vnb, vb As Double Dim Nm-Bulan(l2) As String Public ErlainForm As 6mMiiiu
Private Sub CopyGraph-Click0 MSCkrt 1.Editcopy End Sub Private Sub F~rm~ActivateO gsTahap = 1 End Sub
Private Sub Form_Load() Tahap1.Caption = "Tahap 1 " + gsDatabase Datal .DatabaseName = gsDatabase Datol2.DatabaseName = gsDatabase Data3,DatabaseNatne = gsDatabase Data4.Datab~seNatne= gsDatabase End Sub
Private Sub TbmRu~~ClickO On Error GoTo ErrorHandle: Datal .Re* Data2.Refbh
Data2.Recordset.MoveFirst np = Data2.Recordset.FieIds
Go = Datal.Recordset.Fields(3).Value
NQ) = Datal .Recordset.Fields(4).Value Tabel Lampiran 3 (Lanjutan)
Data1.Recordset.MoveNext P(13) = P(13) + P o q 1 3 ) = G(13) + G(I) N(13) = N(13) + N o Next I pb=P(13)/ 12 gb = P(13) 1 12 nb = N(13) / 12 q13) = 0: pff13)=0 F o r I = 1 To I2 = 0.154 0.05 * (Popb) I pb m = ( l -a01* GO q13) = q13) + ql) p m = P o * ~ / l o o pW3) = pR13) + pfm Next I Z = KAI * (P(13) - N(13)) I pq13) el(13) = 0 ForI= 1 To 12 el(I)=Z*pfm el(13) = el(13) + el(I) Next I ForI= 1 To 12
-
-
vnO=O n l 0 =N o Next I veb = 0.5927 * (P(13) - N(13)) / For I = 1 To 12 ve(I) = veb * qi) 1(q13) I 12) Next I ForI= 1 To 12
vm=veO+vnO ~ 1 =0PQ) - nYI) - el(I)) yo = EM-1 * VO)
* 12)
V(l>
Next I mi(1) =j For K = 1 To 20 F o r I = 2 To 13 mi(I)=y(I-I)*--l)+(l-y(r-l))*pl(I-1) Next I Tabel Lampiran 3 (Lanjutan)
If Abs(mi(l3) - mi(1)) < 0.000001 Then Exit For mi(1) = mi(13)
Next K For1=2 To 13 mi(Q=y(I-l)*mi(I-1)+(1-y(I-l))*pl(I-1) Next I m(13) = 0 F o r I = I T o 12 m O = p 1 0 + tW+1) -miO)/Log(yO) m(13) = m(13) + m(I) Next I m(13) = m(13) / 2 nl(13) = 0 For I = 1 To 12 A(I) = np * PN) If d(I) < 0 Then nl(1) = 0 nl(13) = nl(13) + ni(I) vnO = 0 d(I))/ mu) Next I Dim Rve As Double Rve = 0 ForI=lTol2 Vrn=veO+-(I) Rve = Rve + veO Next I Rve = Rve / 12 nl(13) = d(13) / 12 vnb = (-1) + vn(2) + vn(3) + ~ ( 4 +) ~ ( 5+)~ ( 6+) ~ ( 7+)vn(8)
(Pa-
M10) +
+ -9) +
-11) + vn(12)) / 12 vb = (V(1)-I-V(~)+V(~)+V(~)+V(~)+V(~)+V(~)+V(~)+V( 10) +V(11)+V( 12))/12 W& MSChartI .cOhCount =2 .ROwLabelCount = 12 .RowCount = 12 . C 0 l =~ 1 F o r I = 1 To 12 .Row = I .Column = 1 .Data = =(I); .Column = 2 Tabel Lampiran 3 &anjutan)
.Data = ve(I) Next I End With Data1.Recordset,MoveFirst
.
Data3.Recordset.MoveFirst For I = 1 To 12 Datal .Recordset.Edit Data3.Recordset.Edit Datal .Recordset.Fiel&(S)-Value = vn(l)
Data3.Recordset.FieI&(2).Value = P(1) Data3.Recordset.Fields(3).Value = G O Data3.Recordset.Fields(4).Value = Data3.Recordset.Fields(S).Value = fo Data3.Recod~&.Fieads(6).Value = v e O Data3.Recordset.Fields(7).Value = m i 0 Data3.Recordset.Fields(S).Value = m(l) Data3.Recordset.Fields(9) .Value = n l o Data3.Recordset.Fields(lO).Value = vnQJ Data3.Recordset.Fiilds(l l).Value = V(1)
No
Datal .Recordset,Update Dat&.Recordset.Update Datal .Recordset.MoveNext Data3.Recordset.MoveNext Next I
Data4.Recordset.MoveFirst Data4.Recordset.Edit Data4.Recordset.Fields(3).Value = Rve Data4.Recordset.Update Datal X e M Data2.Rehsh Data3.Re&& Data4.Refresh Exit Sub ErrorHandle: -ox "Terjadi kesalahanperhhngan, Cek Data Anda", vbCritical+ vbOKOnly, "Kesalahan" End Sub Private Sub TbmTahap-Click0 gsTahap = 2 Unload Me Tabel Lampiran 3 (Lanjutan) Tahap2.Slaow
End Sub Private Sub TmbKeluarluarClickO Datal ,Recordset.Close Data2.Recordset.Cfose Data3.Recordset.CIose
Data4.Recordset.Close =
""
fMai~FonnnmuFileClose~Click End Sub
Lampiran 5 . Liitng Program Visual Basic Mmghhmg Luaran Model (Tahap Il).
Private Sub CopyGraphhClickO MSCbart 1.Editcopy End Sub Private Sub Form-IdialkeO gsTahap = 2 End Sub Private Sub Form-Load0 Tahap2.Caption = "Tahap 2 " + gsDatabase Datal .DatabseName= gsDatatmse
End Sub
-
Private Sub TbrnRuu-CIickO On Error GoTo ErrorHandle: D i m I, L As Integer Dim np, PN, ep,veb As Double DimP(13), G(13), Nx(13), a(13), 413), ve(12), m(13), V(13) As Double Dimy(12), P1(13), nri(l3), n$13), e1(13), e2(13), E(13), n1(13), nZ(13), N(13), Px(13), K(t3), n3(I3) As Double Data2.Refresh
Data2.Recordset.MoveFirst np = Data2.Recordset.Fields(O).VahK PN = Data2.RecorW.Fields(I).Value ep = Data2.Recordset.Fields(2)-Vdue veb = DataZ.Remrdset.Fields(3).Vahre Datal .Re*& Datal .Recordset.MoveFirst P(13) = 0: q 1 3 ) = 0: Nx(13) =O: vn(13)=0 For I = I To 12 P O = Datal .Recordset.Fields(2).Value G(I) = Datal .Recordset.Fields(3)-Value NxO = Datal .Recordset.Fields(4).Vdue vn(I) = Datal .Recordset.Fields(S).Value Datal .Recordset.MoveNext P(f 3) = P(13) + P O W13) = W13) + W) Nx(13) =Nx(13) + N x O Tabd lampiran 4 (Lanjutan)
~ ( 1 3=) -13) + vn(I) Next I pb = P(13) / 12: gb = G(13) / 12 f(13) = 0 F o r I = 1 To 12 If W~tbAlbedo= True Then a(I) = Albedo Else a(I) = 0.154 - 0.05 * pb) I pb End If qI) = (1 - a@)) * G(I): q13) = q13) + qI) Next I fb = f(13) / 12 F o r I = 1 To 12 ve@=veb* (qT)/fb) vm=veO+vnO Next I nl(13) = 0: el(13) = 0: Pl(13) = 0: V(13) = 0 F o r I = 1 To 12 n 1 0 = np * P(I)- PN) If nl(I) < 0 Then nl(I) = 0 nl(13) = nl(13) + n l o ) e l m = e p f P(I)*@)/ib el(13) = el(13) + e l 0 PxO = ( P O - n 1 0 -el@)) Px(13) = (P(13) nl(13) - el(13)) PIO=pxO/VO Pl(13) = Pl(13) + P 1 Q y(I)= EM-VO) V(13) = V(13) + V(I) Next I V(13) = V(13) / 12 mi(1) = 300 ForL= I To20 F o r I = 2 To 13 mirn=~-l)*mi~-~)+(l-~-I))*P~(I-l) Next I If Abs(mi(l3) -1)) < 0.00000 1 Then Exit For -1) = mi(13) Next L Tabd Lampiran 4 (Lanjutan)
(Po-
-
-
If L <= 20 Then F o r I = 2 T o 13
mi(I)=y(I- l ) * m i ( I - l ) + ( l - y ( r - l ) ) * P l ( I - 1 ) Next I m(13) = 0 For1 = 1 To 12 mO = PlQ) + (& +( 1).I &(I)) I Log(y(I)) m(13) = m(13) + mo Next I mb = m(13) / 12: K(13) = 0 F o r I = 1 To 13 K O = f i O - C V f l ) *mO) K(13) = K(13) + K(I) Next I e2(13) = 0: nZ(I3) = 0: E(13) = 0: N(13) =; 0: d(13) = 0 For1 = 1 To 12 e 2 0 = ve(I) * m(l): e2(13) = e2(13) + e 2 0 0 2 0 = v@J * d(13) = n2(13) + &(I) E(I)=elQ +e20: E(13)=E(13) + E m N O = n l o + 0 2 0 : N(13) = N(13) + N O = Nx(I) NO: n3(13) = d ( 1 3 ) + d o Next I EIse MsgBox "TidakKonvergen dalam 20 Items?', vbCriticaI + vbOKOnly, "Konvergensi" Exit Sub End If WithMSChartl .Columncount = 2 .RowLabelCount = 12 .Rowcount = I2 .Cohmm = 1 For I = 1 To 12 .Row = I .Column = 1 .Data = N x o .C o h = 2 .Data = N O Next I End with Data3.Recordset.MoveFirst Tabel Lampikin 4 (Lanjutan)
-
a:
do
-
Data3.Recordset.Fielas(3).Vdue
=G
O
Data3.Recordset.Fields(4).Value = Nx(I) Data3.Recordset.FieIds(5).Vdue = vn@) Data3.Recordset.FieIds(6).Value = rn(I) Data3.Records&.Fields(7).Value = JCO Data3.Recordset-Fklds(S).Vahre = E O Data3.Recordset.Fields(9).Value = N O Data3 .Recordset.FieIds(lO)-Value= PI ( I ) Data3.Re~ordset.~eIds(l l).Value = PxO Data3.Recordset.Fields(l2).Value = V O Data3.Recordset.Fields(l3).Value = yO Data3.Recordset.FieIds(14).Value = e l 0 Data3Xecordset.Fie~ds(I S).Value = e 2 0 Data3.Recordset.Fields(l6).Value = n 1 0 Data3.Recordset.Fielas(l7).V~= 1120 Data3.Recordset.Fields(l8)-Value = n 3 O Data3.Recordset.Update Data3.Recordset.MoveNext Next I Datal .Refi-esh;Data2.Refiresh Data3 .Refkesh
Exit Sub ErrorHandle: Msgl3ox "TerjadikesalatLan perbitungan, Cek Data Anda", vbCritical+ vbOKOdy, "Kesalahan" End Sub Private Sub TbmTahapPClickO
Unload M e Tabapl.Show
E d Sub Private Sub TmbKeluarhLarClickO Datal .RecordsetClose Data2.Recordset.Close Data3 .Recordset.CIose gsDatabase = "" Unload M e End Sub
Tabel Lampiran 6 . Nilai r clan R M S I E hujan (mdhai) dan suhu ("C) bulaoan model GCMs dengan data BMG
Tabel Lampiran 7+ Nilai r d m RMSIE hujan (mm/har) dan suhu ("C)bulanan model GCMs dengan data NCAR
Tabel Lampiran 8. Skroring Nilai r dan RMSIE Hujan dan Suhu model GCMs dengao Data BMG
Tabel Lampiran 9. Shoring Nilai r dan RMSIE Hujan dan Suhu model GCW dengan Data NCAR
Tabel Lampiran 10. Rerata Bulanan Komponen Iklim DAS Jeneberang, Walanae dan Saddang --
Tabel Lampiran 1 1. Penggunaan Iahan clan arahan penggunaan lahan DAS Walanae Hulu
5.6~wasauBudidaya TsnamanTahuaan 52367 TanamanSemusim 78537 273501 2846000 1 0 0 . 0 ~ ~ ,Jumlah Sumber:BRLKT Whyah M (1996) *) data telah diaaalisis ulang
pemukimsn
159091
19.1 28.7 100.0
KNA (-)
P m M P E' E" E N' N Nm NdI'
SdisibsktOzvrrSlren3(%) DASWdanaeHulu DASSaddang ' ~ ~ 1 ~GmlL 0 9 CSIR09 0 0 0 6 5 6 47 6 24 -7 1 5 -16
-9
0
'
0 0 24 0 10 -35 1
-7 -8
11 10 10 42 11 24 -58 10 -10
-6
GFDL 0 4 4 4 23 4
12 -32 4 -5 -6
Bnjur Gambar Lamph.an 1. Perbedaanagihen dan jeluk GIltah hujmtahunau Indonesia periode tahun 1931 -1%0 demgsn perbedaanjumlab stasiun hujan amhis: (A) 3540 stasitm hujm dm (B) 214 stasiuu hjau
Gambar
2. fklisih curah hqiau talnman antara 2 X C 9 dsn Ixca model GcMs ( C c W csIR09,GFDL, GISS dan UKMO).
Gambar L a m p i i 3. Selisih suhu rata-rata tahunan antara 2XC02 dengan 1XCOz model GCMs ( CCCh4, CSIRO9, GFDL, GISS clan UKMO).
