Microsoft Word 6735 doc Bé n«ng nghiÖp vµ ph¸t triÓn n«ng th«n tr−êng ®¹i häc thñy lîi B¸o c¸o tæng kÕt ®Ò tµi hîp t¸c theo nghÞ ®Þnh th− víi italia hîp t¸c Nghiªn cøu x©y dùng c«ng nghÖ dù b¸o lò tru[.]
Bộ nông nghiệp phát triển nông thôn trờng đại học thủy lợi Báo cáo tổng kết đề tài hợp tác theo nghị định th với italia hợp tác Nghiên cứu xây dựng công nghệ dự báo lũ trung hạn kết nối với công nghệ điều hành hệ thống công trình phòng chống lũ cho đồng sông hồng - sông thái bình Chủ nhiệm đề tài: pgs, ts vũ minh cát 6735 24/2/2008 hà nội - 2007 B KHOA HỌC VÀ CÔNG NGHỆ ********** BỘ NÔNG NGHIỆP VÀ PHÁT TRIỂN NÔNG THÔN ********** TRƯỜNG ĐẠI HỌC THUỶ LỢI ============================================================= BÁO CÁO KẾT QUẢ ĐỀ TÀI NGHIÊN CỨU KHOA HỌC CÔNG NGHỆ CẤP NHÀ NƯỚC Hợp tác nghiên cứu xây dựng công nghệ dự báo lũ trung hạn kết nối với cơng nghệ điều hành hệ thống cơng trình phịng chống lũ cho đồng sông Hồng – sông Thái Bình BÁO CÁO TỔNG KẾT Cơ quan quản lý: Bộ Khoa học Công Nghệ Bộ nông nghiệp phát triển nơng thơn Cơ quan chủ trì: Trường Đại học Thuỷ lợi Chủ nhiệm đề tài: PGS.TS Vũ Minh Cát Hà nội tháng 11 - 2007 BỘ KHOA HỌC VÀ CÔNG NGHỆ ********** BỘ NÔNG NGHIỆP VÀ PHÁT TRIỂN NÔNG THÔN ********** TRƯỜNG ĐẠI HỌC THUỶ LỢI ============================================================ BÁO CÁO KẾT QUẢ ĐỀ TÀI NGHIÊN CỨU KHOA HỌC CÔNG NGHỆ CẤP NHÀ NƯỚC Hợp tác nghiên cứu xây dựng công nghệ dự báo lũ trung hạn kết nối với công nghệ điều hành hệ thống cơng trình phịng chống lũ cho đồng sơng Hồng – sơng Thái Bình GIAI ĐOẠN 2006 - 2007 BÁO CÁO TỔNG KẾT Cơ quan quản lý: Bộ nông nghiệp phát triển nông thôn Bộ Khoa học Cơng Nghệ Cơ quan chủ trì: Trường Đại học Thuỷ lợi Chủ nhiệm đề tài: PGS.TS Vũ Minh Cát Những người thực hiện: PGS.TS Vũ Minh Cát, ĐHTL PGS.TS Lê Văn Nghinh, ĐHTL TS Nguyễn Lan Châu, TTDBKTTVTW Ths Đỗ Lệ Thủy, TTDBKTTVTW Ths Hoàng Thanh Tùng Ths Nguyễn Hồng Sơn Ths Phạm Xn Hồ NCS Ngơ Lê An Và cộng tác viên khác Hà Nội Tháng 11 - 2007 MỤC LỤC CHƯƠNG I: MỞ ĐẦU .1 1.1 Mở đầu 1.2 Mục tiêu nội dung nghiên cứu 1.3 Phạm vi đề tài .3 1.4 Nội dung phương pháp nghiên cứu 1.4.1 Dự báo mưa trung hạn 1.4.2 Dự báo lũ trung hạn 1.4.3 Tích hợp dự báo mưa lũ trung hạn điều hành hệ thống hồ chứa 1.5 Tổng quan tình hình nghiên cứu thuộc lĩnh vực đề tài 1.5.1 Tình hình nghiên cứu nước 1.5.2 Tình hình nghiên cứu ngồi nước 11 1.6 Đóng góp bên đề tài hợp tác nghiên cứu 13 1.6.1 Phía đối tác Italy 14 1.6.2 Phía đối tác Việt Nam 15 CHƯƠNG II: TỔNG QUAN VỀ LƯU VỰC SƠNG HỒNG –THÁI BÌNH THU THẬP, PHÂN TÍCH VÀ QUẢN LÝ CÁC TÀI LIỆU CƠ BẢN 19 2.1 Tổng quan lưu vực sông Hồng – Thái Bình 19 2.1.1 Tổng quan chung 19 2.1.2 Đặc điểm địa hình 20 2.1.3 Đặc điểm địa chất, thổ nhưỡng thực vật .21 2.1.4 Đặc điểm khí hậu .22 2.1.5 Đặc điểm hệ thống sơng ngịi 25 2.1.6 Đặc điểm dân sinh kinh tế 28 2.2 Thu thập phân tích tài liệu 28 2.2.1 Số liệu khí tượng 28 2.2.1.1 Mạng lưới trạm đo khí tượng 28 2.2.1.2 Đánh giá chung tài liệu khí tượng 29 ii 2.2.1.3 Một số đặc điểm chung khí tượng hệ thống 30 2.2.2 Số liệu thủy văn 31 2.2.2.1 Mạng lưới trạm đo đạc thủy văn .31 2.2.2.2 Đánh giá chung tài liệu khí tượng thuỷ văn .31 2.2.2.3 Đặc điểm thủy văn lưu vực sơng Hồng – Thái Bình 32 2.2.3 Số liệu thủy văn hồ chứa 35 2.2.4 Tài liệu địa hình hệ thống sơng Hồng -Sơng Thái Bình 36 2.2.4.1 Tài liệu địa hình hồ chứa 36 2.2.4.2 Bản đồ số DEM thuộc lưu vực sông Hồng 36 2.2.4.3 Tài liệu địa hình hệ thống sông Hồng .36 2.2.4.4 Tài liệu địa hình khu phân chậm lũ 36 2.2.5 Tài liệu thảm phủ thực vật đồ đất lưu vực sông Hồng .37 2.3 Xây dựng sở liệu 37 2.3.1 Cơ sở số liệu đồ 38 2.3.2 Cơ sở số liệu KTTV 40 2.3.2.1 Giới thiệu phần mềm WRDB 41 2.3.2.2 Xây dựng sở liệu KTTV sử dụng phần mềm WRDB 42 2.3.2.3 Khai thác sở liệu KTTV xây dựng 58 CHƯƠNG III: NGHIÊN CỨU XÂY DỰNG CÔNG NGHỆ DỰ BÁO MƯA TRUNG HẠN 65 3.1 Dự báo mưa Việt Nam 65 3.1.1 Đặc điểm hình thành mưa tổ hợp kiểu hình thời tiết gây mưa lớn lưu vực sông Hồng 65 3.1.1.1 Các hình điển hình 65 3.1.1.2 Đặc điểm hình thời tiết gây mưa lũ lớn lưu vực sông Hồng66 3.1.2 Đánh giá trạng phương pháp dự báo mưa Việt nam 69 3.2 Thử nghiệm mơ hình BOLAM Italia 70 3.2.1 Giới thiệu mơ hình BOLAM 70 iii 3.2.1.1 Động lực học phương pháp số 71 3.2.1.2 Vật lý sơ đồ tham số hóa vật lý 73 3.2.1.3 Điều kiện ban đầu điều kiện biên .74 3.2.1.4 Tiền xử lý hậu xử lý mơ hình BoLAM 76 3.2.2 Nghiên cứu ứng dụng mơ hình BoLAM cho Việt Nam 76 3.2.2.1 Thay đổi cấu trúc thư mục mơ hình BoLAM 79 3.2.3 Dự báo thử nghiệm mơ hình BOLAM mùa lũ 2007 83 3.3 Một số kết qủa so sánh mơ hình BoLAM mơ hình HRM 95 3.4 Kết luận kiến nghị 102 CHƯƠNG IV: NGHIÊN CỨU XÂY DỰNG CÔNG NGHỆ DỰ BÁO LŨ 103 4.1 Các phương pháp dự báo thủy văn thường dùng Việt Nam 103 a Các phương pháp thường dùng dự báo thủy văn 103 Phương pháp xu 103 2.Phương pháp lưu lượng mực nước tương ứng 104 Phương pháp lượng trữ 104 Nhóm phương pháp phân tích thống kê 105 Phương pháp mạng Nơ ron thần kinh – Mạng trí tuệ nhân tạo (Artificial Neural Network) 105 Phương pháp sử dụng mơ hình thủy văn thủy lực .106 Phương pháp kỹ thuật số .106 b.Tổng kết phương pháp thường dùng dự báo ngắn hạn, trung hạn dài hạn Việt Nam 107 Các phương pháp dự báo ngắn hạn 108 Các phương pháp dự báo trung hạn .108 Các phương pháp dự báo dài hạn 109 4.2 Phân tích đặc điểm mưa lũ hệ thống sông Hồng 109 4.2.1 Đặc điểm mưa gây lũ lưu vực hệ thống sông Hồng .109 4.2.2 Vai trò trữ lượng nước lưu vực hình thành lũ 113 iv 4.2.3 Tổ hợp lũ sông Đà, Lô, Thao 116 Đặc điểm mưa gây lũ lưu vực sông Thao: .117 Dịng chảy lũ sơng Lơ .121 4.3 Xây dựng công nghệ dự báo lũ ngắn hạn hệ thống sông Hồng 125 4.3.1 Phương pháp hồi quy đa biến 128 a Cơ sở phương pháp Hồi quy đa biến 128 b Các giả thiết phương pháp quy trình thực phương pháp 128 c Xây dựng ma trận tương quan .129 d Phương trình dự báo phép thử 130 a Cấu trúc mạng ANN 138 b Q trình qt xi 139 c Hàm kích hoạt 139 d Chuẩn hoá (Normalization) 140 e Lựa chọn sơ đồ mạng nơ ron thần kinh ban đầu 141 f Phương pháp quét ngược (Back Propagation Method) 141 g Một số yếu tố ảnh hưởng đến trình học theo phương pháp lan truyền ngược sai số .143 h Vấn đề học chưa đủ học thuộc mạng 144 i Một số giải pháp cho vấn đề học mạng: 144 j Giới thiệu phần mềm neuroSolutions .145 4.3.3 Xây dựng phương án dự báo lũ ngắn hạn cho vị trí khác hệ thống sông Hồng 147 4.3.3.2 Kết dự báo lũ cho sông Thao: 150 4.3.3.3 Kết dự báo lũ cho sông Lô - Gâm: 160 4.3.3.2 Kết dự báo lũ cho sông Đà: 166 4.4 Xây dựng công nghệ dự báo lũ trung hạn hệ thống sông Hồng 176 4.4.1 Ứng dụng mơ hình đối tác đề nghị DIMOSHONG 176 4.4.1.1 Giới thiệu mơ hình DIMOSHONG 176 v 4.4.1.2 Xây dựng mơ hình DIMOSHONG cho lưu vực hệ thống sông Hồng 178 4.4.1.3 Kết mô mô hình DIMOSHONG cho trận lũ lớn cho hệ thống sông Hồng 186 4.4.2 Ứng dụng mơ hình HEC-HMS .191 4.4.2.1 Giới thiệu mơ hình HEC-HMS 191 4.4.2.2 Xây dựng mơ hình HEC-HMS cho lưu vực hệ thống sông Hồng.192 b Kết nối liệu liệu mưa vào dự báo dòng chảy 193 c Tóm tắt số kết tính tốn dự báo dịng chảy từ mưa .196 4.5 Xây dựng công nghệ dự báo lũ cho hệ thống sông Hồng 197 4.5.1 Lập trình kết nối xây dựng cơng nghệ 197 4.5.2 Sử dụng công nghệ dự báo ngắn hạn .198 4.5.3 Sử dụng công nghệ dự báo trung hạn mơ hình DIMOSHONG 199 4.5.4 Sử dụng cơng nghệ dự báo trung hạn HEC-HMS 201 4.6 Thử nghiệm công nghệ dự báo lũ trung hạn cho mùa lũ 2007 207 4.6.1 Dự báo thử nghiệm sử dụng mơ hình DIMOSHONG 207 4.6.2 Dự báo thử nghiệm sử dụng mơ hình HEC-HMS 210 CHƯƠNG V: CÔNG NGHỆ KẾT NỐI DỰ BÁO MƯA VỚI DỰ BÁO LŨ TRUNG HẠN 213 5.1 Xây dựng phần mềm kết nối dự báo mưa với dự báo lũ .214 5.1.1 Lập trình đọc ghi kết dự báo từ BOLAM file PPF – file số liệu vào cho DIMOSHONG 214 5.1.2 Lập trình giao diện xử lý số liệu đầu vào hiển thị kết 224 1) Chương trình Nhập liệu 224 2) Chương trình hiển thị kết chạy mơ hình DIMOSHONG 236 CHƯƠNG VI: KẾT NỐI MƠ HÌNH DỰ BÁO MƯA LŨ TRONG ĐIỀU HÀNH HỆ THỐNG HỒ CHỨA PHÒNG CHỐNG LŨ 257 6.1 Xây dựng mơ hình kết nối dự báo mưa lũ với điều hành hệ thống hồ chứa phòng lũ 257 6.1.1 Giới thiệu tóm tắt mơ hình HEC-HMS khả điều tiết hồ chứa .257 vi 6.1.2 Mô hệ thống sông Hồng phục vụ điều hành hệ thống hồ chứa HEC-HMS .261 6.2 Lựa chọn hệ thống kịch mưa dòng chảy lũ phục vụ điều hành thử nghiệm theo mơ hình dự báo lũ trung hạn 265 6.2.1 Những để thiết lập hệ thống kịch 265 6.2.2 Lựa chọn kịch mưa lũ 266 6.3 Phương án vận hành hệ thống hồ chứa phòng lũ theo thời gian dự kiến dự báo hạn trung (từ đến ngày) 268 6.3.1 Điều tiết hồ thời kỳ lũ sớm 268 6.3.2 Điều tiết hồ thời kỳ lũ vụ .269 6.3.3 Điều tiết hồ thời kỳ lũ muộn 270 6.3.4 Vận hành đảm bảo an tồn cơng trình .271 6.4 Kết vận hành thử nghiệm hệ thống hồ chứa theo kịch đề 272 Trận lũ lịch sử năm 1971 xảy thời kỳ lũ sớm từ 15 tháng đến 15 tháng 272 Trận lũ lịch sử năm 1971 xảy thời kỳ lũ vụ từ 16 tháng đến 25 tháng 275 Trận lũ lịch sử năm 1971 xảy thời kỳ lũ muộn từ 26 tháng đến 15 tháng 279 6.5 Thử nghiệm mơ hình vận hành mùa lũ 2007 283 CHƯƠNG VII: KẾT LUẬN VÀ KIẾN NGHỊ 287 7.1 Kết luận 287 a Về khoa học công nghệ 287 b Về mặt hợp tác quốc tế 290 c Về mặt đào tạo .290 d Về phục vụ sản xuất .291 7.2 Kiến nghị 291 TÀI LIỆU THAM KHẢO 293 vii MỤC LỤC HÌNH VẼ Hình 2.1: Vị trí địa lý lưu vực sơng Hồng – Thái Bình lãnh thổ Việt Nam 19 Hình 2.2: Giao diện trang đầu Chương trình quản lý sở liệu 38 Hình 2.3: CSDL quản lý đồ DEM có độ phân giải khác nước 39 Hình 2.4: CSDL quản lý toàn đồ quét tỷ lệ 1:50.000 nước 39 Hình 2.5: Sơ đồ cấu trúc CSDL WRDB 42 Hình 2.6: Sơ đồ mạng lưới trạm khí tượng thủy văn hệ thống sơng Hồng Thái Bình 43 Hình 2.7: Giao diện bảng support 48 Hình 2.8: Các yếu tố đo 49 Hình 2.9: Phân nhóm yếu tố đo 50 Hình 2.10: Phân nhóm yếu tố đo 50 Hình 2.11: Sơ đồ nhập liệu cho CSDL KTTV 57 Hình 2.12: Ví dụ nhập liệu mực nước lưu lượng trạm Sơn Tây 58 Hình 2.13: Dữ liệu mực nước lưu lượng trạm Sơn Tây bảng Import 58 Hình 2.14: Bảng thống kê theo trạm thơng số 59 Hình 2.15: Lọc liệu mực nước trạm có mã trạm 74165 60 Hình 2.16: Sắp xếp liệu từ lớn đến nhhỏ theo mã trạm theo kết 60 Hình 2.17 Lựa chọn liệu mực nước trạm Sơn Tây QBE Query 61 Hình 2.18: Đường trình mực nước trạm Mường Tè Lai Châu trục tọa độ 61 Hình 2.19: Mối tương quan mực nước trạm Mường Tè Lai Châu 62 Hình 2.20: Đường dẫn chọn vùng liệu quan tâm 62 Hình 2.21: Chọn vùng liệu quan tâm 63 Hình 2.22: Kiểm tra vùng liệu chọn 63 Hình 3.1: Biểu đồ quy trình chạy mơ hình BoLAM máy tính PC-Linux 78 Hình 3.2: Cấu trúc thư mục mơ hình BoLAM 80 viii Drofa, O., P.Malguzzi, G Fioravanti, A Buzzi, 2003: Modelling of the thermal and water exchange in the vegetated soil and interaction with the atmosphere Proceeding of the conference of "Associazione Italiana di Agrometeorologia" (Italian Agrometeorological Association) 2003, Bologna May 29-30 2003 (in Italian) Emmett, W., W., 1970 The hydraulics of overland flow on hillslopes, Geological Survey Professional Paper 662-A Georgelin, M., P Bougeault, T Black, N Brzovic, A Buzzi, J Calvo, V Cassé, M Desgagné, R ElKhatib, J F Geleyn, T Holt, S.-Y Hong, T Kato, J Katzfey, K Kurihara, B Lacroix, F Lalaurette, Y Lemaitre, J Mailhot, D Majewski, P Malguzzi, V Masson, J Mcgregor, E Minguzzi, T Paccagnella and C Wilson, 2000: The second COMPARE exercise: a model intercomparison using a case of a typical mesoscale orographic flow, the PYREX IOP3 Q J Roy Meteor Soc., 126, 991-1030 Gyakum, J.R., M Carrera, D.-L Zhang, S Miller, J Caveen, R Benoit, T Black, A Buzzi, C Chouinard, M Fantini, C Folloni, J.J Katzfei, Y.-H Kuo, F Lalaurette, S Low-Nam, J Mailhot, P Malguzzi, J.M McGregor, M Nakamura, G Tripoli and C Wilson., 1996: A regional model intercomparison using a case of explosive oceanic cyclogenesis Wea Forecasting, 11, 521-543 Horton,R.E., 1938 The interpretation and application of runoff plot experiments with reference to soil erosion problem , J Soil Sci Soc Am 56,275-370 Hydrology Earth System Sciences (HESS), 7(6), 799-811 Kain, J S and Fritsch, J M., 1990: A one-dimensional entraining/detraining plume model and its application in convective parameterization J Atmos Sci., 47, 2784–2802 Kain, J.S., 2004: The Kain-Fritsch convective parametrization: an update J App Meteorol., 43, 170-181 La Barbera., P., 1990 Problemi di preannuncio di eventi estremi: modelli geomorfo-climatici: meccanica del debris-flow, Rapporto 1988 -Linea del CNR-GNDCI, Genova, 11 pp Lacassin R, Leloup PH, Trinh PT, Tapponnier P, Unconformity of red sandstones in north Vietnam: field evidence for Indosinian orogeny in northern Indocina?, Terra Nova, 10:106-11, 1998 Lagouvardos, K., and V Kotroni, 2005: Improvement of high-resolution weather forecasts through humidity adjustment, based on satellite data Quart J Roy Meteor Soc., 131, 2695-2712 Lagouvardos, K., V Kotroni, A Koussis, C Feidas, A Buzzi, P Malguzzi, 2003: The meteorological model BOLAM at the National Observatory of Athens: Assessment of two-year operational use J Appl Meteor., 42, 1667-1678 Le Bac Huynh, “Mơ hình tính tốn dự báo lũ sơng Hồng” “Mathematical modelling and flood forecasting in Red River”, Journal of Hydro – Meteorology, - 2001 Le Van Anh, “Lũ lụt đồng Bắc Bộ thập kỷ vừa qua” “Floods in delta of North Vietnam in recent decades”, Journal of Hydro – Meteorology, – 2000 Leheman, R., 1993: On the choice of relaxation coefficients for Davies’ lateral boundary scheme for regional weather prediction models Meteorol Atmos Phys., 52, 1-14 Leloup PH, Lacassin R, Tapponier P, Schärer U, Zhong D, Xiaohan L, Liangshang Z, Shaocheng J, Trinh PT, The Ailao Shan-Red River shear zone (Yunnan, China), Tertiary trasnform boundary of Indochina, Tectonophysics, 251:3-84, 1995 Louis, J F., 1979: A parametric model of the vertical eddy fluxes in the atmosphere Boundary-Layer Meteorol., 17, 18–202 Malguzzi, P., and N Tartaglione, 1999: An economical second order advection scheme for explicit numerical weather prediction Quart J Roy Met Soc., 125, 2291-2303 Malguzzi, P., Grossi, G., Buzzi, A., Ranzi, R and Buizza, R., 2006: The 1966 ‘century’ flood in Italy: a meteorological and hydrological revisitation, J Geophys Res., 111, D24106, doi:10.1029/2006JD007111 79 Mancini, M., La modellazione della risposta idrologica: effetti della variabilità spaziale e della scala di rappresentazione del fenomeno di assorbimento, Tesi di Dottorato, Politecnico di Milano, 1990 Miglietta, M.M., and A Buzzi, 2004: A numerical study of moist stratified flow regimes over isolated topography Quart J Roy Meteor Soc., 130, 1749-1770 Ministry of Science and Technology of Vietnam (MOST), Atlas Vietnam, 1999 Mlawer, E.J., S.J Taubman, P.D Brown, M.J Iacono, and S.A Clough, 1997: Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave J Geophys Res., 102D, 16, 663-16, 682 Monin, A S and Obukhov, A., 1955: Basic laws of turbulence mixing in the surface layer of the atmosphere Trudy Geophys Inst AN SSSR, 24, 163–187 Morcrette, J.-J., 1991: Radiation and cloud radiative properties in the ECMWF operational weather forecast model J Geophys Res., 96D, 9121-9132 Nagata, M., L Leslie, H Kamahori, R Nomura, H Mino, Y Kurihara, E Rogers, R L Elsberry, B K Basu, A Buzzi, J Calvo, M Desgagne, M D'Isidoro, S.-Y Hong, J Katzfey, D Majewski, P Malguzzi, J McGregor, A Murata, J Nachamkin, M Roch, C Wilson, 2001: A Mesoscale Model Intercomparison: A Case of Explosive Development of a Tropical Cyclone (COMPARE III) J Meteorol Soc Japan, 79, 999-1033 NCHMF-National Centre for Hydro-meteorological Forecasting, “Đánh giá hình thời tiết sinh lũ lớn phục vụ dự báo cảnh báo trước khả có lũ lớn, lũ cực hạn hệ thống sơng Hồng – Thái Bình”, “Assessment of the shapes of climate causing big flood to improve the extreme flood warning and long-term forecasting in Red – Thai Binh river system”, State Project, 2001 Ngo Dinh Tuan et al, Report on Current status of the water resources of Vietnam, ADB TA 3528VIE Project, 2005 Orlandini, S and Rosso, R., 1996, Diffusion wave modelling of distributed catchment dynamics, J Hydrol Eng., ASCE, 1(3), 103-113 Pappenberger, F., Beven, K.J., Hunter, N.M., Bates, P.D., Gouweleeuw, B.T., Thielen, J., de Roo, A.P.J., 2005: cascading model uncertainty from medium range weather forecasts (10 days) through a rainfall-runoff model to flood inundation predictions within the European Flood Forecasting System (EFFS), Hydrol Earth System Sci 9(4), 381-393 Philip, J.R., 1958 The theory of infiltration, Soil Sci., 85, 333-337 Pilotti, M., and Menduni, G., 2001 Beginning of sediment transport of incoherent grains in shallow shear flows, J Hydraulic Res., IAHR, 39, 115-124 Ponce, V.M., and Theurer, F.D., 1982, Accuracy criteria in diffucion routing, J Hydraul Div., ASCE, 112(8), 747-757 Pressman D.YA., 1994 Chislennaja model' gidrotermicheskikh processov v pochve kak chast' skhemy mezomasshtabnogo prognoza (A numerical model of hydrothermal processes in soil as a part of a mesoscale weather forecast scheme) − Meteorologija i gidrologija (Meteorology and Hydrology), No.11, 62-75, in Russian Ranzi, R., B Bacchi and G Grossi, Runoff measurements and hydrological modelling for the estimation of rainfall volumes in an alpine basin, Quart Journal Royal Meteorol Soc., Vol 129 Part B, n 588, 653-672, 2003 Ranzi, R., M Bochicchio And B Bacchi, Effects on floods of recent afforestation and urbanisation in the Mella River (Italian Alps), Hydrol Earth System Sci., 6(2), 239-265, 2002 Ranzi, R., M Zappa and B Bacchi, Hydrological aspects of the Mesoscale Alpine Programme: findings from field experiments and simulations, Quart Journal Royal Meteorol Soc., 2007 (in press) Richard, E, Buzzi, A and Zaengl, G., 2007: Quantitative precipitation forecasting in mountainous regions: The push ahead obtained by the Mesoscale Alpine Programme (MAP) Quarterly J Royal Meteorol Soc, (in press) 80 Richard, E., A Buzzi and G Zängl, 2006: Quantitative precipitation forecasting in mountainous regions: The push ahead obtained by the Mesoscale Alpine Programme (MAP) Submitted to Quart J Roy Meteorol Soc Richard., E., S Cosma, R Benoit, P Binder, A Buzzi and P Kaufmann, 2003: Intercomparison of mesoscale meteorological models for precipitation forecasting Roger F, Leloup PH, Jolivet M, Lacassin R, Trinh PT, Brunel M, Seward D, Long and complex thermal hystory fo the Song Chay metamorphic dome (Northern Vietnam) by multi-system geochronology, Tectonophysics, 321:449-466, 2000 Schultz, P., 1995: An explicit cloud physics parameterization for operational numerical weather prediction Mon Wea Rev., 123, 3331–3343 Steinshouer DW, Qiang J, McCabe PJ, Ryder RT, Maps Showing Geology, Oil and Gas Fields, and Geologic Provinces of the Asia Pacific Region, USGS-CRERT Open-File Report 97-470F, on website http://pubs.usgs.gov/of/1997/ofr-97-470/OF97-470F/index.html#TOC Tran Tuat, Tran Thanh Xuan, Nguyen Duc Nhat, “Sơng ngịi Việt Nam” “Hydrogeography of rivers in Vietnam”, Science and Technology Publishing House, Hanoi, 1987 Trinh Quang Hoa, “Quy hoạch phịng chống lũ đồng sơng Hồng” “Assessment of the Planning of Flood Control for the Red River Delta”, State Project, 2001 USDA-NRCS-CED, Urban Hydrology for Small Watershed, TR-55, 1986 USGS-EROS, Global DEM GTOPO30, on website http://edc.usgs.gov/products/elevation/gtopo30/gtopo30.html USGS-EROS, Global Land Cover Characterisation, V.2.0, on website http://edcsns17.cr.usgs.gov/glcc/ Vu Minh Cat, Le Van Nghinh, Ngo Le An, “Lũ sông Hồng dự báo lũ Hà Nội” “Red River Flood and its forecasting at Hanoi Station”, Workshop on Natural environment, sustainable protection and conservation: Italy – Vietnam cooperation perspectives, 2004 Zampieri, M., P Malguzzi and A Buzzi, 2005: Sensitivity of quantitative precipitation forecasts to boundary layer parameterization: a flash flood case study in theWestern Mediterranean Natural Hazard Earth System Sci., 5, 603-612 Zonn SV, Tropical and Subtropical Soil Science, MIR, Moscow, 1986 81 Annexes A Hydrometic stations Vector Layer Name Vector Layer Type Reference System Reference Units Unit Distance ID/Value Type Number of Features : : : : : : : HYDROMETRICSTATIONS SONG HONG Text utm-48n m Integer 20 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 1 538406.773606 Hoa Binh 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 2 528750.119081 Lo@Tuyen Quang 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 3 528750.119081 Thao@Phu Tho 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 4 559329.525077 Hong@Son Tay 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 5 593932.537125 Hong@Ha Noi 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption : : : : 6 485777.848246 Thao@Yen Bai 2304165.419501 2415066.099773 2367652.040816 2339525.056689 2321845.238095 2400058.526536 82 Style Rotation Justification X Justification Y : : : : 135 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 7 393458.118225 Thao@Lao Cai 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 8 538650.0 @Na Hang 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 9 537650.0 @Bac Me 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 10 10 528650.0 2439050.0 @Chiam Hoa 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 11 11 497650.0 @Ha Giang 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 12 12 488650.0 @Vinh Tuy 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption : : : : 13 13 505650.0 @Ham Yen 2488970.771270 2471050.0 2516050.0 2521050.0 2464050 2441050.0 83 Style Rotation Justification X Justification Y : : : : 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 14 14 502650.0 @Thac Ba 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 15 15 526650.0 2389050.0 Lo@Vu Quang 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 16 16 545650.0 2348050.0 Lo@Viet Tri 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X : : : : : : : 17 17 251650.0 @Muong Te 90 Justification Y : 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 18 18 337650 2441050.0 @Lai Chau 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 19 19 377650 Ban Cung 90 0.5 Feature Number ID or Value Coordinates (X,Y) : 20 : 20 : 407650.0 2406050.0 2490050 2409050.0 2371050.0 84 Caption Style Rotation Justification X Justification Y : : : : : @Ta Bu 90 0.5 Vector Layer Name Vector Layer Type Reference System Reference Units Unit Distance ID/Value Type Number of Features : : : : : : : RAINGAUGE2002_NEW Text utm-48n m Integer 55 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 1 246111.861024 PacMa (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 2 276922.865923 MuongTe (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 3 343896.812438 TamDuong (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 4 318030.142797 SinHo (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption : : : : 5 330826.878313 PaTan (Da) Style Rotation Justification X Justification Y : : : : 90 0.5 B Raingauge stations 2503560.609169 2476864.291886 2479735.319627 2468943.745566 2538934.820948 85 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 6 310213.133802 NamGiang (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 7 250671.524285 NamPo (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 8 322619.741693 NamMuc (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 9 357520.433906 BinhLu (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 10 10 293455.589500 ThanUyen (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 11 11 401845.524521 2418372.058035 MuCangChai (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 12 12 378006.139758 BanCung (Da) 90 0.5 2461654.108360 2445874.703232 2419058.871879 2468531.428294 2428640.562858 2408942.543343 86 Feature Number : 13 ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : 13 365385.189403 MuongSai (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 14 14 384187.268058 SonLa (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 15 15 408674.662886 MaiSon (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 16 16 309106.849250 LaiChauKT (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 17 17 401761.134113 TaBu (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 18 18 461956.387635 PhuYen (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 19 19 475040.454556 BanChieng (Da) 90 0.5 2383580.652019 2357604.975870 2346749.833177 2441364.149187 2370699.230224 2349860.368061 2311099.095244 87 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 20 20 497918.393197 Km46 (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 21 21 487507.887193 Km22 (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 22 22 534685.207253 HoaBinh (Da) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 23 23 496921.492658 HaGiangTV (Lo) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 24 24 486285.982463 BacQuang (Lo) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 25 25 487980.203141 VinhTuy (Lo) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 26 26 503439.081022 HamYen (Lo) 90 0.5 2297800.095404 2294484.510224 2301893.731259 2523594.353791 2488177.765843 2462348.356480 2440206.448316 88 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 27 27 509217.895657 YenMinh (Lo) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 28 28 568347.240114 BaoLac (Lo) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 29 29 537758.743329 NaHang (Lo) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 30 30 575455.612816 ChoRa (Lo) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 31 31 527499.828507 ChiemHoa (Lo) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 32 32 443321.251134 PhoRang (Lo) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 33 33 456996.064315 BaoYen (Lo) 90 0.5 2552378.637513 2538139.860863 2469768.959716 2482821.445468 2447609.331576 2458757.002464 2434730.364285 89 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 34 34 501723.548632 ThacBa (Lo) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 35 35 520670.628940 2412548.181447 TuyenQuang (Lo) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 36 36 526745.698590 VuQuang (Lo) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 37 37 391997.850005 2488550.400129 LaoCaiTV (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 38 38 417467.782871 VanBan (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 39 39 431265.682024 BaoHa (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 40 40 448359.002853 NgoiHut (Thao) 90 0.5 2403311.279488 2384886.476064 2441898.492622 2451425.708399 2423687.298627 90 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 41 41 463830.478463 2414420.579870 NgoiThia (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 42 42 449965.068584 VanChan (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 43 43 484486.292922 2401474.003323 YenBaiTV (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 44 44 505190.650696 MinhDai (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 45 45 525944.583844 2346147.542873 ThanhSon (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 46 46 541483.460597 VietTri (Hong) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 47 47 551918.609027 SonTay (Hong) 90 0.5 2386786.610300 2340594.114532 2357247.321066 2336985.810130 91 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 48 48 581403.755836 HaNoi (Hong) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 49 49 118973.986320 China1 (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 50 50 214624.179598 China2 (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 51 51 128517.094095 China3 (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 52 52 149958.398905 China4 (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 53 53 238607.410840 China5 (Thao) 90 0.5 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 54 54 431385.909656 China6 (Thao) 90 0.5 2324191.786069 2481795.147218 2469760.407742 2771103.327426 2818239.933728 2768553.944520 2766227.189437 92 Feature Number ID or Value Coordinates (X,Y) Caption Style Rotation Justification X Justification Y : : : : : : : : 55 55 327826.469462 China7 (Thao) 90 0.5 2545613.655934 93