Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.

Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc hơi do bức xạ mặt trời tại lưu vực sông Cầu.

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ TÀI NGUYÊN VÀ MÔI TRƯỜNG VIỆN KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ o0o NGUYỄN VĂN HÙNG NGHIÊN CỨU ỨNG DỤNG TƯ LIỆU VIỄN THÁM XÁC ĐỊNH LƯỢNG NƯỚC BỐC HƠI DO BỨC XẠ MẶT TRỜI TẠI LƯU VỰC SÔNG CẦU LUẬN ÁN TIẾN SĨ KỸ THUẬT HÀ NỘI, NĂM 2017 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ TÀI NGUYÊN VÀ MÔI TRƯỜNG VIỆN KHOA HỌC ĐO ĐẠC VÀ BẢN ĐỒ o0o NGUYỄN VĂN HÙNG NGHIÊN CỨU ỨNG DỤNG TƯ LIỆU VIỄN THÁM XÁC ĐỊNH LƯỢNG NƯỚC BỐC HƠI DO BỨC XẠ MẶT TRỜI TẠI LƯU VỰC SÔNG CẦU LUẬN ÁN TIẾN SĨ KỸ THUẬT NGÀNH: KỸ THUẬT TRẮC ĐỊA – BẢN ĐỒ MÃ SỐ: 62.52.05.03 NGƯỜI HƯỚNG DẪN KHOA HỌC: TSKH LƯƠNG CHÍNH KẾ TS ĐÀO NGỌC LONG HÀ NỘI, NĂM 2017 i LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu riêng thân tơi Tồn trình nghiên cứu tiến hành cách khoa học, số liệu, kết trình bày luận án xác, trung thực chưa cơng bố cơng trình khác Tác giả luận án Nguyễn Văn Hùng ii LỜI CẢM ƠN Trong trình thực luận án“Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc xạ Mặt Trời Lưu vực Sông Cầu”, Tôi nhận nhiều giúp đỡ, tạo điều kiện tập thể lãnh đạo, nhà khoa học, cán bộ, chuyên viên Cục Viễn thám quốc gia; tập thể Ban Lãnh đạo Viện Khoa học Đo đạc Bản đồ Tơi xin bày tỏ lịng cảm ơn chân thành giúp đỡ Tơi xin bày tỏ lòng biết ơn sâu sắc tới TSKH Lương Chính Kế, TS Đào Ngọc Long thầy giáo trực tiếp hướng dẫn bảo cho Tơi hồn thành luận án Tôi xin chân thành cảm ơn bạn bè, đồng nghiệp Tôi công tác Cục Viễn thám quốc gia, Viện Khoa học Đo đạc Bản Đồ, Trường Đại học Mỏ địa chất Hà Nội, Đại học Khoa học Tự nhiên Hà Nội gia đình động viên, khích lệ, tạo điều kiện giúp đỡ Tơi suốt q trình thực hoàn thành luận án Tác giả luận án Nguyễn Văn Hùng iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC CÁC KÍ HIỆU, CÁC CHỮ VIẾT TẮT vii DANH MỤC CÁC BẢNG BIỂU ix DANH MỤC CÁC HÌNH VẼ VÀ SƠ ĐỒ xii MỞ ĐẦU 1 Tính cấp thiết luận án Mục tiêu nghiên cứu Nội dung nghiên cứu Phương pháp nghiên cứu Phạm vi nghiên cứu Ý nghĩa khoa học thực tiễn luận án 6.1 Ý nghĩa khoa học: 6.2 Ý nghĩa thực tiễn: Luận điểm bảo vệ luận án Những điểm luận án CHƯƠNG 1: TỔNG QUAN VỀ BỐC THOÁT HƠI NƯỚC DO BỨC XẠ MẶT TRỜI VÀ PHƯƠNG PHÁP NGHIÊN CỨU 1.1 Các khái niệm bốc thoát nước yếu tố ảnh hưởng đến bốc thoát nước 1.1.1 Các khái niệm bốc thoát nước 1.1.2 Các yếu tố ảnh hưởng đến bốc thoát nước 1.2 Tổng quan phương pháp truyền thống xác định lượng bốc thoát nước 10 1.2.1 Nhóm phương pháp trực tiếp xác định lượng bốc nước ngồi thực địa .10 1.2.2 Nhóm phương pháp thực nghiệm sử dụng số liệu quan trắc khí tượng 11 1.3 Tổng quan nghiên cứu ứng dụng tư liệu viễn thám xác định lượng bốc thoát nước lượng xạ Mặt Trời 13 1.3.1 Tổng quan tình hình nghiên cứu giới 13 1.3.2 Tổng quan kết nghiên cứu nước 17 1.4 Đánh giá kết nghiên cứu nước giới 19 1.5 Những vấn đề phát triển luận án 23 CHƯƠNG 2: CƠ SỞ KHOA HỌC XÁC ĐỊNH LƯỢNG BỐC THOÁT HƠI NƯỚC THỰC TẾ BỀ MẶT ĐẤT DỰA VÀO NĂNG LƯỢNG BỨC XẠ MẶT TRỜI CHIẾT XUẤT TỪ TƯ LIỆU ẢNH VIỄN THÁM .24 2.1 Đặc tính phản xạ phổ đối tượng tự nhiên 24 2.1.1 Đặc tính phản xạ phổ thực vật 24 iv 2.1.2 Khả phản xạ phổ thổ nhưỡng 26 2.1.3 Khả phản xạ phổ nước 28 2.1.4 Sự phản xạ phổ đối tượng đô thị 30 2.2 Cơ sở khoa học ứng dụng tư liệu viễn thám xác định lượng bốc thoát nước từ lượng xạ Mặt Trời 31 2.2.1 Phương pháp cân lượng bề mặt đất 31 2.2.2 Vai trò tư liệu viễn thám việc xác định tham số để tính lượng bốc thoát nước 32 2.2.2.1 Xác định tham số tính xạ ròng hấp thụ bề mặt đất 32 2.2.2.2 Xác định tham số tính Thơng nhiệt bề mặt đất .33 2.2.2.3 Xác định tham số tính lượng nhiệt cảm ứng H 33 2.2.2.4 Xác định tham số tính lượng nhiệt ẩn LE 33 2.3 Phương pháp viễn thám xác định lượng bốc thoát nước thực tế bề mặt đất 35 2.3.1 Chiết xuất lượng xạ ròng trung bình ngày hấp thụ mặt đất .35 2.3.1.1 Hiệu chỉnh xạ ảnh vệ tinh 35 2.3.1.2 Tính tốn suất phân sai bề mặt đất α .41 2.3.1.3 Bức xạ sóng ngắn (RS↓ ) 41 2.3.1.4 Bức xạ sóng dài (RL↑ ) 42 2.3.1.5 Bức xạ tới sóng dài RL↓ 44 2.3.1.6 Bức xạ ròng hấp thụ bề mặt đất trung bình Rni 45 2.3.1.7 Bức xạ rịng trung bình ngày hấp thụ bề mặt đất Rnd .45 2.3.1.8 Phương pháp FAO-56 PM tính xạ rịng trung bình ngày 45 2.3.2 Xác định lượng bốc thoát nước thực tế bề mặt đất theo ngày 46 2.3.2.1 Xác định lượng bốc thoát nước thực tế bề mặt đất theo Phương pháp số cân lượng xạ bề mặt giản lược S-SEBI .46 2.3.2.2 Tính tốn bốc nước thực tế bề mặt đất theo phương pháp Priestley-Taylor 48 2.4 Đề xuất ứng dụng phương pháp xác định lượng bốc thoát nước thực tế bề mặt đất từ lượng xạ Mặt Trời chiết xuất từ ảnh viễn thám .49 2.4.1 Lựa chọn phương pháp S-SEBI .49 v 2.4.2 Đề xuất giải pháp nâng cao độ xác phương pháp S-SEBI 51 2.5 Kết luận Chương 51 CHƯƠNG 3: THỰC NGHIỆM XÁC ĐỊNH LƯỢNG BỐC THOÁT HƠI NƯỚC THỰC TẾ BỀ MẶT ĐẤT CỦA LƯU VỰC SÔNG CẦU BẰNG SỬ DỤNG TƯ LIỆU ẢNH LANDSAT-7 52 3.1 Đặc điểm vị trí địa lý tự nhiên .52 3.1.1 Vị trí địa lý .52 3.1.2 Đặc điểm địa hình 52 3.1.3 Đặc điểm khí hậu 54 3.1.4 Thổ nhưỡng 58 3.1.5 Thảm phủ thực vật 59 3.1.6 Mạng lưới sông suối 59 3.2 Thực nghiệm xác định lượng bốc thoát nước thực tế bề mặt đất Lưu vực Sông Cầu sử dụng tư liệu ảnh LANDSAT-7 theo phương pháp l ựa chọn S-SEBI 59 3.2.1 Dữ liệu đầu vào 59 3.2.2 Đề xuất quy trình tính tốn lượng bốc thoát nước thực tế bề mặt đất theo phương pháp S-SEBI 63 3.2.3 Chiết xuất lượng xạ rịng trung bình ngày hấp thụ bề mặt đất Rnd từ ảnh LANDSAT-7 .65 3.2.3.1 Tiền xử lý ảnh viễn thám 65 3.2.3.2 Tạo ảnh suất phân sai bề mặt đất α 65 3.2.3.3 Tạo ảnh số thực vật NDVI 67 3.2.3.4 Tạo ảnh phát xạ bề mặt εo 68 3.2.3.5 Tính ảnh nhiệt độ bề mặt Ts nhiệt độ khơng khí Ta 68 3.2.3.6 Tạo ảnh xạ ròng hấp thụ bề mặt đất trung bình ngày 73 3.2.4 Xác định lượng bốc thoát nước thực tế bề mặt đất ETa Lưu vực Sơng Cầu dựa vào lượng xạ rịng trung bình ngày chiết xuất từ ảnh viễn thám theo phương pháp SSEBI 78 3.2.4.1 Tỷ phần bốc thoát nước tức thời ETFi .78 3.2.4.2 Bốc thoát nước thực tế bề mặt đất trung bình ngày ETa (mm/ngày) 79 3.2.5 Thành lập phương trình xác định hệ số “liên hệ c” .82 3.2.6 Khảo sát xác định tham số a sử dụng phương pháp Priestley-Taylor phù hợp với điều kiện địa hình khí hậu Lưu vực Sông Cầu Việt Nam 87 3.2.7 Khảo sát bốc thoát nước tham chiếu ETo theo phương pháp FAO-56 PM từ lượng xạ ròng chiết xuất từ phương pháp khác 90 vi 3.2.8 Khảo sát mối quan hệ bốc thoát nước tham chiếu ETo_Rnd_VT với tham số Albedo, NDVI, Ts_DEM Rnd_ VT tính tốn từ ảnh viễn thám 93 3.2.9 Tính ảnh bốc nước tham chiếu ETo_Rnd_VT dựa vào số thực vật NDVI 94 3.2.10 Xác định hệ số trồng Kc 97 3.2.11 Khảo sát bốc thoát nước thực tế bề mặt đất ETa_VT3 theo phương pháp S-SEBI với tham số ETo_Rnd_VT , Rnd_VT , Ts_DEM, NDVI, Albedo chiết xuất từ ảnh viễn thám 103 3.2.12 Khảo sát mối quan hệ bốc thoát nước thực tế bề mặt đất ETa_VT3 theo phương pháp S-SEBI với tham số khí tượng đo đạc ngồi thực địa 104 3.2.13 So sánh đánh giá bốc thoát nước thực tế bề mặt đất tính theo phương pháp S-SEBI từ xạ rịng trung bình ngày tính tốn theo phương pháp khác 107 KẾT LUẬN VÀ KIẾN NGHỊ 111 DANH MỤC CÁC CƠNG TRÌNH KHOA HỌC ĐÃ CÔNG BỐ 113 TÀI LIỆU THAM KHẢO 114 PHỤ LỤC - Phụ lục 1.1: Tính ETo theo FAO-56 PM sử dụng Rnd_FAO1 (trong xạ Mặt Trời tới Rs xác định theo Ăngstrom – công thức 2.32b), ngày 04/11/2000 - Phụ lục 1.2: Tính ETo theo FAO-56 PM sử dụng Rnd_FAO2 (trong xạ Mặt Trời tới Rs xác định theo công thức thực nghiệm Việt Nam – công thức 2.32c), ngày 04/11/2000 - Phụ lục 1.3: Tính ETo theo FAO-56 PM sử dụng xạ rịng trung bình ngày chiết xuất từ ảnh LANDSAT Rnd_VT – công thức 2.31), ngày 04/11/2000 - Phụ lục 2.1: Tính ETo theo FAO-56 PM sử dụng Rnd_FAO1 (trong xạ Mặt Trời tới Rs xác định theo Ăngstrom – công thức 2.32b), ngày 23/11/2001 - Phụ lục 2.2: Tính ETo theo FAO-56 PM sử dụng Rnd_FAO2 (trong xạ Mặt Trời tới Rs xác định theo công thức thực nghiệm Việt Nam – công thức 2.32c), ngày 23/11/2001 - Phụ lục 2.3: Tính ETo theo FAO-56 PM sử dụng xạ rịng trung bình ngày chiết xuất từ ảnh LANDSAT Rnd_VT – công thức 2.31), ngày 23/11/2001 - - vii DANH MỤC CÁC KÍ HIỆU, CÁC CHỮ VIẾT TẮT Kí hiệu Tiếng Anh Nghĩa tiếng Việt E T ET ETo Evaporation Transpiration Evapotranspiration Reference evapotranspiration Bốc nước Thoát nước Bốc thoát nước Bốc thoát nước tham chiếu Bốc thoát nước tham chiếu xác định theo phương pháp FAO-56 PM sử dụng xạ ròng chiết xuất từ ảnh viễn thám Bốc thoát nước tiềm Bốc thoát nước thực tế bề mặt đất Bốc thoát nước thực tế bề mặt đất theo phương pháp SSEBI (khi chưa hiệu chỉnh tham số) Bốc thoát nước thực tế bề mặt đất theo phương pháp SSEBI (sau hiệu chỉnh nhiệt độ Ts DEM) Bốc thoát nước thực tế bề mặt đất theo phương pháp SSEBI (sau hiệu chỉnh tham số ) Bốc thoát nước thực tế bề mặt đất tính theo phương pháp Priestley-Taylor Phương pháp số cân lượng xạ bề mặt giản lược Phương pháp Chỉ số cân lượng bề mặt Phương pháp cân lượng bề mặt cho đất Phương pháp hệ thống cân lượng bề mặt Phương pháp lập Bản đồ bốc thoát nước kiểm định nội hóa Mơ hình mã hai nguồn Tỷ phần bốc thoát nước tức thời Nhiệt độ bề mặt đất ETo_Rnd_VT PET ETa Potential evapotranspiration Actual evapotranspiration ETa_VT1 ETa_VT2 ETa_VT3 ETa_PT Actual evapotranspirationPriestley Taylor S-SEBI Simplified Surface Energy Balance Index SEBI Surface Energy Balance Index Surface Energy Balance Algorithms for Land Surface Energy Balance System Mapping Evapotranspiration with Internalized Calibration SEBAL SEBS METRIC TSM ETFi Two-sources models Evapotranspiration Fraction Ts Surface Temperature viii Nhiệt độ khơng khí gần mặt đất Bức xạ rịng Bức xạ sóng ngắn Difference Chỉ số thực vật NDVI Ta Near Surface Temperature Rn Rs NDVI Net Radiation Short radiation Normalized Vegetation Index Albedo surface Surface Emissivity Crop coefficient Land + Satellite Gross Domestic Product Food and Agriculture Organization of the United56 Penman – Monteith α 𝜀𝑜 Kc LANDSAT GDP FAO - 56 PM (RMSE) SSTP (MAE) SSTĐTB Root mean square error Mean absolute error Suất phân sai bề mặt đất Phát xạ bề mặt Hệ số trồng Ảnh vệ tinh NASA, MỸ Tổng sản phẩm quốc nội Phương pháp xác định bốc thoát nước tham chiếu tiêu chuẩn Tổ chức Lương thực Nông nghiệp Liên Hiệp Quốc đề xuất năm 1990 Sai số trung phương Sai số tuyệt đối trung bình 110 Sử dụng mơ hình số độ cao (DEM) để hiệu chỉnh nhiệt độ bề mặt Ts chiết xuất từ kênh ảnh nhiệt LANDSAT, kết hợp với việc xác định hệ số “liên hệ c” độ xác xác định lượng bốc thoát nước thực tế bề mặt đất nâng cao lên 5% (minh chứng mục 3.2.3, mục 3.2.4, mục 3.2.5, bảng 3.14a, bảng 3.14b, bảng 3.15a bảng 3.15b) - Minh chứng điểm luận án: Đề xuất xác định hệ số “liên hệ c” nhằm nâng cao độ xác tính tốn lượng bốc nước thực tế bề mặt đất điều kiện địa hình khí hậu Lưu Vực Sông Cầu Việt Nam (minh chứng mục 3.2.5, công thức (3.1), bảng 3.15a bảng 3.15b); Hệ số a phương pháp Priestley-Taylor khảo nghiệm điều kiện địa hình khí hậu Lưu vực Sông Cầu Việt Nam nằm khoảng 0,951,05 cho phép tính tốn bốc nước thực tế bề mặt đất với sai số tuyệt đối trung bình khơng lớn 20% so với kết đo ngoại nghiệp (minh chứng mục 3.2.6, bảng 3.18 bảng 3.19); Dựa vào bốc thoát nước thực tế bề mặt đất tính theo phương pháp viễn thám bốc nước tham chiếu tính theo phương pháp FAO-56 PM giúp cho việc xác định hệ số trồng hỗn hợp Kc quy mô diện rộng nhanh (minh chứng mục 3.2.10); Đề xuất xây dựng áp dụng thành công quy trình tính tốn lượng bốc nước thực tế bề mặt đất theo mơ hình S-SEBI (Mơ hình số cân lượng xạ bề mặt giản lược) Lưu vực Sông Cầu Việt Nam (minh chứng mục 3.2.2, hình 3.2 mục 3.2.4) 111 KẾT LUẬN VÀ KIẾN NGHỊ A KẾT LUẬN Kết nghiên cứu luận án “Nghiên cứu ứng dụng tư liệu viễn thám xác định lượng nước bốc xạ Mặt Trời Lưu vực Sông Cầu” chứng minh đầy đủ luận điểm tính mà luận án đặt Thông qua kết nghiên cứu thực nghiệm Luận án, nghiên cứu sinh rút số kết luận sau: Có thể thay tính tốn xạ rịng theo phương pháp truyền thống (phương pháp FAO-56 PM) phương pháp viễn thám S-SEBI, sai số tuyệt đối trung bình xác định lượng bốc nước thực tế bề mặt đất ETa_VT3 Lưu vực Sông Cầu (cho hai trường hợp sau hiệu chỉnh nhiệt độ Ts hệ số “liên hệ c”) nhỏ 20% Sử dụng mơ hình số độ cao (DEM) để hiệu chỉnh nhiệt độ bề mặt Ts chiết xuất từ kênh ảnh nhiệt LANDSAT, kết hợp với việc xác định hệ số “liên hệ c” nâng cao độ xác xác định lượng bốc thoát nước thực tế bề mặt đất điều kiện địa hình khí hậu Lưu vực Sông Cầu lên 5% Sử dụng phương pháp Priestley-Taylor cho điều kiện địa hình khí hậu Lưu vực Sơng Cầu với ứng dụng ảnh LANDSAT, hệ số a phương pháp nhận tối ưu khoảng 0,95-1,05, sai số tuyệt đối trung bình tính tốn bốc nước thực tế bề mặt đất khơng lớn 20% Dựa vào bốc nước thực tế bề mặt đất tính theo phương pháp viễn thám bốc thoát nước tham chiếu tính theo phương pháp FAO-56 PM giúp cho việc xác định hệ số trồng hỗn hợp Kc quy mô diện rộng nhanh Nghiên cứu sinh xây dựng ứng dụng thành cơng quy trình tính toán lượng bốc thoát nước thực tế bề mặt đất theo phương pháp S-SEBI cho điều kiện địa hình khí hậu Lưu vực Sơng Cầu Việt Nam Có thể áp dụng quy trình tính tốn bốc thoát nước thực tế bề mặt đất cho Lưu vực Sông Việt Nam (với điều kiện xác định tham số phù hợp với điều kiện khu vực) Bốc thoát nước thực tế bề mặt đất ETa_VT3 xác định theo phương pháp S-SEBI có mối tương quan mạnh với nhiệt độ Ts_TĐ độ ẩm khơng khí RH đo ngồi thực địa, nhiệt độ Ts_VT tính tốn từ ảnh viễn thám Đây dẫn quan trọng cho cán khoa học đo đạc nhiệt độ Ts_TĐ độ ẩm khơng 112 khí RH ngồi thực địa tính tốn nhiệt độ Ts_VT từ ảnh viễn thám phục vụ tính tốn bốc nước thực tế bề mặt đất theo phương pháp viễn thám Bốc nước tham chiếu ETo_Rnd_VT có mối tương quan mạnh với số thực vật NDVI Do đó, xây dựng hàm hồi quy ETo_Rnd_VT số thực vật NDVI phục vụ cơng tác tính tốn bốc nước tham chiếu ETo_Rnd_VT cho tồn Lưu vực Sông Cầu B KIẾN NGHỊ Cần có thêm nghiên cứu ứng dụng phương pháp S-SEBI tính tốn lượng bốc nước thực tế bề mặt đất ETa_VT3 cho vùng có điều kiện địa hình khí hậu khác để kiểm chứng lại Mơ hình Quy trình nhằm ứng dụng phương pháp viễn thám quy mô lãnh thổ Việt Nam Bốc nước tham chiếu ETo có mối tương quan mạnh với số thực vật NDVI Trong tương lai, nghiên cứu sinh tiêp tục nghiên cứu tính tốn hệ số trồng Kc thơng qua số thực vật NDVI điều kiện Việt Nam 113 DANH MỤC CÁC CƠNG TRÌNH KHOA HỌC ĐÃ CƠNG BỐ Minh Le, Ke Luong Chinh, Tuong Tran Ngoc, Hung Nguyen Van, Son Le Minh, (2013), Essessing the accuracy of land surface evapotranspiration estimated by Makkink’s model based on solar radiation extracted from modis data, Proceedings of the 9th International Conference on Geoinformation for Disaster Management, Hanoi, Vietnam, 9-11 Dec 2013 TS Nguyễn Xuân Lâm, ThS Nguyễn Văn Hùng, CN Vũ Hữu Liêm, Ứng dụng ảnh vệ tinh Modis khảo sát mối quan hệ nhiệt đô bề mặt lượng xạ Mặt Trời khu vực phía bắc Việt Nam, Tạp chí Khoa học Đo đạc Bản đồ, Số 18, tháng 12/2013; trang 29-37 Tran Ngoc Tuong, Nguyen Van Hung, Pham Van Manh, Luong Chinh Ke, (2014), Effect of daily net radiation model extracted from MODIS on daily reference evapotranspiration, Proceedings of the International Symposium on GeoInformatics for Spatial-Infrastructure Development in Earth and Allied Sciences (GIS-IDEAS), Da Nang, Dec 6-9, 2014; page, 397-405 ThS Nguyễn Văn Hùng, Trần Ngọc Tưởng, Lương Chính Kế, (2014), “Xác định hệ số tuyến tính phương pháp Makkink cho tính tốn lượng bốc nước khu vực phía bắc Việt Nam sử dụng tư liệu viễn thám”, Kỷ yếu Hội thảo Khoa học Công nghệ, Viện khoa học Đo đạc Bản đồ , Hà Nội, tháng 7, 2014; trang 223-228 ThS Nguyễn Văn Hùng, ThS Trần Ngọc Tưởng, ThS Dương Đình Trường, CN Nguyễn Ngọc Anh, TSKH Lương Chính Kế, (2016), “Phương pháp ước tính bốc nước thực tế tích hợp số thực vật NDVI ảnh Modis với số liệu khí tượng cho địa hình miền Bắc Việt Nam”, Kỷ yếu Hội thảo khoa học Đo đạc Bản đồ với ứng phó biến đổi khí hậu, Viện khoa học Đo đạc Bản đồ, Hà Nội Luong Chinh Ke, Nguyen Le Dang, Nguyen Van Hung, Tran Ngoc Tuong, (2016), “Proposed empirical model of crop coefficient Kc from MODIS-NDVI for estimating actual land surface”, Geoinformatics for Natural Resources, Hazards and Sustainability - GIS IDEAS 2016, Ha Noi 12-16, November ThS Nguyễn Văn Hùng, (2016), “Xác định lượng bốc thoát nước thực tế Lưu vực Sông Cầu ứng dụng tư liệu ảnh LANDSAT”, Tạp chí Khoa học Đo đạc Bản đồ, tháng 12/2016 114 TÀI LIỆU THAM KHẢO A Tài liệu tiếng Việt Lê Vân Anh, Trần Anh Tuấn (2014), Nghiên cứu nhiệt độ bề mặt sử dụng phương pháp tính tốn độ phát xạ từ số thực vật, Tạp chí khoa học trái đất, NXB Viện KH&CN Việt Nam, tập 36, số 2, tr 184-92 Nguyễn Văn Bình cộng (2016), “Nghiên cứu, phát triển hệ thống quản lý, giám sát vùng biên giới cảng biển sử dụng ảnh viễn thám VNREDSat-1 tương đương phục vụ công tác đảm bảo an ninh quốc gia” Đề tài cấp Nhà nước, Bộ Khoa học Công nghệ ThS Nguyễn Văn Hùng, Trần Ngọc Tưởng, Lương Chính Kế (2014), Xác định hệ số tuyến tính phương pháp Makkink cho tính tốn lượng bốc nước khu vực phía bắc Việt Nam sử dụng tư liệu viễn thám, Kỷ yếu Hội thảo Khoa học Công nghệ, Viện khoa học Đo đạc Bản đồ , Hà Nội, trang 223-228 Lương Chính Kế, cộng (2011), “Ứng dụng công nghệ viễn thám để nghiên cứu khả phát giám sát số thành phần ô nhiễm mơi trường khơng khí khu vực thị khu công nghiệp”, Đề tài cấp năm 2009-2011 Lương Chính Kế, cộng (2012), “Nghiên cứu Xây dựng quy trình cơng nghệ hiệu chỉnh chuẩn hoá liệu ảnh vệ tinh quang học đa thời gian, đa đầu thu, đa độ phân giải nhằm nâng cao chất lượng ảnh”, Đề tài cấp 2010-2012 Lương Chính Kế cộng (2013), “Đánh giá độ xác lượng bốc nước bề mặt địa hình theo phương pháp Makkink dựa vào số liệu xạ Mặt Trời chiết xuất từ ảnh vệ tinh Modis”, Kỷ yếu Hội thảo Khoa học Công nghệ, trường ĐH Tài nguyên Môi trường Hà Nội, lần thứ nhất, 2013 Lương Chính Kế, Trần Ngọc Tưởng, Phạm Văn Mạnh (2014), “Khảo sát bốc thoát nước tiềm bề mặt địa hình theo phương pháp PriestleyTaylor xạ ròng chiết xuất từ ảnh vệ tinh Modis”, Kỷ yếu Hội thảo Khoa học Công nghệ, Viện khoa học Đo đạc Bản đồ, Hà Nội Tiêu chuẩn Quốc gia TCVN 8641-2011, “cơng trình thủy lợi kỹ thuật tưới tiêu cho lương thực thực phẩm” Phạm Vọng Thành, Nguyễn Trường Xuân, (2003), Công nghệ viễn thám Bài giảng dành cho học viên cao học Trường Đại học Mỏ Địa chất, Hà Nội 10 Lê Anh Tuấn (2009), Giáo trình “Hệ thống tưới tiêu”, tài liệu giảng dạy trường đại học Cần Thơ 11 Trần Thị Vân, Hoàng Thái Lan, Lê Văn Trung (2011), Nghiên cứu thay đổi nhiệt độ bề mặt thị tác động hóa q trình thị hóa Thành Phố Hồ Chí Minh phương pháp viễn thám, Tạp chí Các khoa học Trái Đất, NXB Viện KH&CN Việt Nam, tập 33, số 3, tr 347-359 12 Viện Quy hoạch Thủy lợi (2007), báo cáo tổng hợp dự án “Rà soát, cập nhật bổ sung quy hoạch thủy lợi Lưu vực Sông Cầu – sông Thương” B Tài liệu tiếng Anh 13 Allen, R.G., L Pereira, D Raes, and M Smith (1998), Crop Evapotranspiration Food and Agriculture Organization of the United Nations, Rome, Italy, ISBN 92-5-104219-5 290p 115 14 Allen, R G., Tasumi, M., Morse, A., and Trezza, R (2005), “A LANDSAT based energy balance and evapotranspiration model in western US water rights regulation and planning.” J of Irrig., and Drain Sys., 19, 251–268 15 Allen, R., Tasumi, M., Morse, A., Trezza, R., Wright, J., Bastiaanssen, W., Kramber, W., Lorite, I., Robison, C (2007a), Satellite-Based Energy Balance for Mapping Evapotranspiration with Internalized Calibration (METRIC)— Applications Journal of Irrigation and Drainage Engineering 133(4): 395406 16 Atmospheric correction mudule: QUAC and FLAASH user’s guide 17 Ångström, Anders (1920), Applications of heat radiation measurements to the problems of the evaporation from lakes and the heat convection at their surfaces Geografiska Annaler, H 3, 16 pp Stockholm 18 Allen, R G., Tasumi, M., and Trezza, R., (2007), “Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC) odel.” J Irrig Drain Eng., 133_4_, 380–394 19 Alexandris, S, Stricevic, R and Petkovic, S (2008) Comparative analysis of reference evapotranspiration from the surface of rainfed grass in central Serbia, calculated by six empirical methods against the Penman-Monteith formula European Water 21(22): 17-28 20 Bastiaanssen, W G M., Menenti, M., Feddes, R A., and Holtslag, A A M (1998a), “A remote sensing surface energy balance algorithm for land (SEBAL) 1: Formulation.” J Hydrol., 212–213, 198–212 21 Bastiaanssen WGM, Pelgrum H, Wang J, Ma Y, Moreno JF, Roerink GJ (1998b), “A remote sensing surface energy balance algorithm for land (SEBAL) Validation” J Hydrol (Amst) 212/213:213–229 22 Bastiaanssen, W.G.M Allen, R G., Tasumi, M., and Trezza, R (2002), SEBAL “Surface energy balance algorithm for land” Version1 23 Bastiaanssen, W.G.M (1995), Regionalization of surface flux densities and moisture indicators in composite terrain Ph.D Thesis, Wageningen Agricultural University, Wageningen, The Netherlands 24 Bastiaanssen, W.G.M., (2000), SEBAL-based sensible and latent heat fluxes in the irrigated Gediz Basin, Turkey Journal of Hydrology, 229:87-100 25 Bastiaanssen, W.G.M., E.J.M Noordman, H Pelgrum, G Davids, B.P Thoreson and R.G Allen., (2005), SEBAL model with remotely sensed data to improve water resources management under actual field conditions Journal of irrigation and drainage engineering, 131 (1):85-93 26 Er Anju Bala , Dr Ram Karan Singh , Dr V Gayathri (2013), Review on different Surface Energy Balance algorithms for Estimation of Evapotranspiration through Remote Sensing, International Journal of Emerging Technology and Advanced Engineering, ISSN, 2250-2459 27 Chen J, F; Yen H-F., Lee C-H., Lo W-C., (2005), Optimal comparison of empirical equations for estimating potencial evapotranspiration in Taiwan XXXI IAHR CONGRESS; September 11~16, 2005, Seoul, Korea 28 C.-Y Xu1* and V P Singh (2000), “Evaluation and generalization of radiation-based methods for calculating evaporation”, Hydrol Process., Vol 14, 339-349 116 29 Choudhury, B.J.; Menenti, M (1993), Parameterization of Land Surface Evaporation by means of Location Dependent Potential Evaporation and Surface Temperature Range; Department for Environment, Food and Rural Affairs (Defra): London, UK; Volume 212, pp 561–568 30 Droogers, P., Allen, R.G (2002), Estimating reference evapotranspiration under inaccurate data conditions Irrigation and Drainage Systems 16: 33–45 31 Dalton, (1802), Experiments and observations to determine whether the quantity of rain and dew is equal to the quantity of water carried off by rivers and raised by evaporation; with an inquiry into the origin of springs (Cited by Brutsaert, 1982, p 31) 32 Drexler, J.Z., Snyder, R.L., Spano, D., Paw U, K.T (2004), A review of models and micrometeorological methods used to estimate wetland evapotranspiration Hydrological Processes 18(11): 2071-2101 33 FAO Irrigation and Drainage Paper No 56, (2006): Guidelines for computing crop water requirement, book 34 Irmak, A., Ratcliffe, I., Ranade, P., Hubbard, K., Singh, K.R (2011), Estimation of Land Surface Evapotranspiration with a Satellite Remote Sensing Procedure Great Plains Research 21: 73-88 35 Gowda, P.H., Chávez, J.L., Colaizzi, P.D., Evett, S.R., Howell, T.A., Tolk, J.A (2007), Remote sensing based energy balance algorithms for mapping ET: Current status and future challenges American Society of Agricultural and Biological Engineers 50(5): 1639-1644 36 Gómez, M., Olioso, A., Sobrino, J.A., Jacob, F (2005), Retrieval of evapotranspiration over the Alpilles/ReSeDA experimental site using airborne POLDER sensor and a thermal camera Remote Sensing of Environment 96: 399 – 408 37 Gibson, L.A., Măunch, Z., Engelbrecht, J 2011 Particular uncertainties encountered in using a pre-packaged SEBS model to derive evapotranspiration in a heterogeneous study area in South Africa Hydrology and Earth System Sciences 15: 295–310 38 Hossein Tabari, Mark E Grismer & Slavisa TrajkovicPaulson, (2005) “Comparative analysis of 31 reference evapotranspiration methods under humid conditions” Irrigation Science, ISSN 0342-7188 39 Jackson R D et al., (1983), Estimation of the daily evapotranspiration from one-time-of-day measurements Agric Water Manage 7(3), pp 351-362 40 Jun, X., Bingfang, W., Nana, Y., Yuan, Z., Shufu, L (2010), Estimation and Validation of Land Surface Evaporation Using Remote Sensing and Meteorological Data in North China Selected Topics in Applied Earth Observations and Remote Sensing, IEEE Journal of 3(3): 337-344 41 Jetse D Kalma, Tim R McVicar E Matthew F McCabe (2008), Estimating Land Surface Evaporation: A Review of Methods Using Remotely Sensed Surface Temperature Data, Surv Geophys 29:421–469 42 Lott, R.B., Hunt, R.J (2001), Estimating evapotranspiration in natural and constructed wetlands Wetlands 21(4): 614–628 43 Minh Le, Ke Luong Chinh, Tuong Tran Ngoc, Hung Nguyen Van, Son Le Minh (2013), Essessing the accuracy of land surface evapotranspiration estimated by Makkink’s model based on solar radiation extracted from modis 117 data, Proceedings of the 9th International Conference on Geoinformation for Disaster Management, Hanoi, Vietnam, 9-11 Dec 44 Menenti, M., Jia, L., Su, H (2003), On SEBI-SEBS validation in France, Italy, Spain, USA and China International Commission on Irrigation and Drainage, New Delhi, India 45 Michael J Johnson (1995), Summer Evapotranspiration Rates, by BowenRatio and Eddy-Correlation Methods, in Boulder Flat and in Maggie Creek Area, Eureka County, Nevada, 1991-1992 Water-Resources Investigations Report 95-4079 Carson City, Nevada 46 Makkink GF (1957), Testing the Penman formula by means of lysimeters Journal of the Institution of Water Engineers 11: 277–288 47 Mutiga, J.K., Su, Z., Woldai, T (2010), Using satellite remote sensing to assess evapotranspiration: Case study of the upper Ewaso Ng’iro North Basin, Kenya International Journal of Applied Earth Observation and Geoinformation 12, Supplement 1(0): S100-S108 48 Norman, J.M.; Becker, F (1995), Terminology in thermal infrared remote sensing of natural surfaces Agric For Meteorol 77, 153–166 49 Nagler, P.L., Cleverly, J., Glenn, E., Lampkin, D., Huete, A., Wan, Z 2005 Predicting riparian evapotranspiration from MODIS vegetation indices and meteorological data Remote Sensing of Environment 94(1): 17-30 50 Ruhoff, A.L., Paz, A.R., Collischonn, W., Aragao, L.E.O.C., Rocha, H.R., Malhi, Y.S (2012), A MODIS-based energy balance to estimate evapotranspiration for clear-sky days in Brazilian Tropical Savannas Remote Sensing 4: 703-725 51 Roerink, G.J.; Su, Z.; Menenti, M (2000), S-SEBI A simple remote sensing algorithm to estimate the surface energy balance Phys Chem Earth B 25, 147–157 52 Priestley CHB, Taylor RJ (1972), On the assessment of surface heat fluxes and evaporation using large-scale parameters Monthly WeatherReview 100: 81–92 53 Penman HL (1948), Natural evaporation from open water, bare and grass Proceedings of the Royal Society of London, Series A: Mathemat ical, Physical and Engineering Sciences 193: 120–145 54 Su, Z (2009), Surface Energy Balance System for estimation of turbulent heat fluxes and evapotranspiration 55 Santos, C.; Lorite, I.J.; Tasumi, M.; Allen, R.G.; Fereres, E (2007), Integrating satellite-based evapotranspiration with simulation models for irrigation management at the scheme level Irrig Sci., 26, 277–288 56 Su, H., M.F McCabe, E F Wood, Z Su and J Prueger., (2005), Modeling evapotranspiration during SMACEX: comparing two approaches for local and regional scale prediction Journal of Hydrometeorology, 6(6): 910-922 57 Su, Z., Timmermans, W., Gieske, A., Jia, L., Elbers, J A., Olioso, A., Timmermans, J., Van Der Velde, R., Jin, X., Van Der Kwast, H., Nerry, F., Sabol, D., Sobrino, J A., Moreno, J and Bianchi, R., (2008), Quantification of land-atmosphere exchanges of water, energy and carbon dioxide in space and time over the heterogeneous Barrax site Vol 29, 5215–5235 118 58 Sobrino, J.A.; Gómez, M.; Jiménez-Muñoz, J.C.; Oliosob, A.; Chehbouni, G (2005), A simple algorithm to estimate evapotranspiration from DAIS data: Application to the DAISEX campaigns J Hydrol, 315, 117–125 59 C Shoko, (2013), Use of Remote sensing for estimating total evaporation, Submitted in partial fulfilment of the requirements for the degree of MSc Bioresource Systems, Centre for Water Resources Research University of KwaZulu-Natal Pietermaritzburg 60 Sobrino, J.A.; Gómez, M.; Jiménez-Moz, J.C.; Olioso, (2007), A Application of a simple algorithm to estimate the daily evapotranspiration from NOAA-AVHRR images for the Iberian Peninsula Remote Sens Environ 110, 139–148 61 Swinbank, W.C (1951), The measurement of vertical transfer of heat and water vapor by eddies in the lower atmosphere: Journal of Meteorology, v 8, p 135-145 62 Su, Z, (2002), The Surface Energy Balance System (SEBS) for estimation of turbulent heat fluxes, Hydrology and Earth System Sciences, 6(1), 85-99 63 Tran Ngoc Tuong, Nguyen Van Hung, Luong Chinh Ke, Pham Van Manh, (2014), “Effects of daily net radiation model extracted from modis on daily reference evapotranspiration” The International Conference on GeoInformatics for Spatial-Infrastructure Development in Earth & Allied Sciences (GIS-IDEAS), DaNang 64 Tasumi, M.; Trezza, R.; Allen, R.G.; Wright, J.L (2005), Operational aspects of satellite-based energy balance models for irrigated crops in the semi-arid U.S Irrig Drain Sys 19, 355–376.Y Oku, H Ishikawa, 65 Van de Griend, A A and Owe, M., (1993), On the relationship between thermal emissivity and the normalized difference vegetation index for natural surfaces, Int J Remote Sens, 14, 1119–1131 66 Valor, E and Caselles, V., (1996), Mapping land surface emissivity from NDVI: application to european, african and south American areas , Remote Sens Environ, 57, 167–184, doi: 10.1016/0034-4257(96)00039-9 67 Yuei-An Liou * and Sanjib Kumar Kar, (2014), Evapotranspiration Estimation with Remote Sensing and Various Surface Energy Balance Algorithms—A Review Energies, 7, 2821-2849; doi:10.3390/en7052821 68 Zhang, R (2016), A review of remote sensing based actual evapotranspiration estimation, https://www.researchgate.net/publication/304708984 -1- PHỤ LỤC Phụ lục 1.1: Tính ETo theo FAO-56 PM sử dụng Rnd_FAO1 (trong xạ Mặt Trời tới Rs xác định theo Ăngstrom – công thức 2.32b), ngày 04/11/2000 Rn-G (MJ/m2/d) Cao độ (m) Uz (m/s) U2 (m/s) TT B (oC) Tmin (oC) Tmax (oC) eT (kPa) 9.90331302 2.9097276 20.55 13.8 27.3 1.57591 9.92247496 2.9097276 20.65 13.3 28 1.52542 10.2712924 17 3.6371595 21.25 14.8 27.7 1.68131 9.1158565 685 4.3645914 17.2 12.2 22.2 1.41931 10.2337331 36 2.9097276 20.9 13.2 28.6 1.51549 9.35356349 138 2.9097276 19.95 11.1 28.8 1.31973 es-ea (kPa) P (kPa) Gama γ (kPa/ oC) 0.9065 101.3 0.89769 eT max(kPa) es (kPa) RH eo(T min) ea (kPa) (%) (kPa) 3.6238213 2.599864569 70 2.419 1.693 3.7749796 2.650198368 72 2.434 1.753 3.7095384 2.695422474 76 2.525 1.919 2.6728283 2.046067573 54 1.962 1.06 3.9088829 2.712187044 73 2.472 1.804 3.9544265 2.637080598 79 2.331 1.842 Lamda λ (MJ kg-1) tu_delta MauSo_delta 900/(T+273) (o C ) Tu_ETo Mau_ETo ETo (mm/ngày) 0.067 Delta Δ (kPa/ o C) 0.149 2.436 9913.4 66486.62 3.06591722 1.1470419 0.283066 4.05 101.2 0.067 0.15 2.444 9974.7 66538.2 3.064873148 1.14573799 0.283808 4.04 0.77608 101.1 0.067 0.155 2.446 10349 66848.1 3.058623619 1.22924271 0.305189 4.03 0.98636 93.46 0.062 0.124 2.444 8042 64770.25 3.101309442 1.29159326 0.278544 4.64 0.90779 100.9 0.067 0.152 2.452 10129 66667.24 3.062266077 1.17700945 0.285385 4.12 0.79555 99.68 0.066 0.144 2.444 9552.7 66177.56 3.072196621 1.02227756 0.276214 3.70 -2Phụ lục 1.2: Tính ETo theo FAO-56 PM sử dụng Rnd_FAO2 (trong xạ Mặt Trời tới Rs xác định theo công thức thực nghiệm Việt Nam – công thức 2.32c), ngày 04/11/2000 Rn-G (MJ/m2/d) Cao độ (m) Uz (m/s) Tmax (oC) eT (kPa) eT max (kPa) es (kPa) RH (%) eo(T min) (kPa) ea (kPa) U2 (m/s) 12.635262 2.9097276 20.55 13.8 27.3 1.57591 3.6238213 2.599864569 70 2.419 1.693 12.8646837 2.9097276 20.65 13.3 28 1.52542 3.7749796 2.650198368 72 2.434 1.753 13.0186433 17 3.6371595 21.25 14.8 27.7 1.68131 3.7095384 2.695422474 76 2.525 1.919 11.1970743 685 4.3645914 17.2 12.2 22.2 1.41931 2.6728283 2.046067573 54 1.962 1.06 12.3181531 36 2.9097276 20.9 13.2 28.6 1.51549 3.9088829 2.712187044 73 2.472 1.804 14.3961104 138 2.9097276 19.95 11.1 28.8 1.31973 3.9544265 2.637080598 79 2.331 1.842 es-ea (kPa) P (kPa) Gama γ (kPa/ oC) 0.9065 101.3 0.89769 TT B (oC) Tmin (oC) Lamda λ (MJ kg-1) tu_delta MauSo_delta 900/(T+273) (o C ) Tu_ETo Mau_ETo ETo (mm/ngày) 0.067 Delta Δ (kPa/ o C) 0.149 2.436 9913.4 66486.62 3.06591722 1.3132384 0.283066 4.64 101.2 0.067 0.15 2.444 9974.7 66538.2 3.064873148 1.32569215 0.283808 4.67 0.77608 101.1 0.067 0.155 2.446 10349 66848.1 3.058623619 1.40278148 0.305189 4.60 0.98636 93.46 0.062 0.124 2.444 8042 64770.25 3.101309442 1.39702401 0.278544 5.02 0.90779 100.9 0.067 0.152 2.452 10129 66667.24 3.062266077 1.30622445 0.285385 4.58 0.79555 99.68 0.066 0.144 2.444 9552.7 66177.56 3.072196621 1.31925509 0.276214 4.78 -3- Phụ lục 1.3: Tính ETo theo FAO-56 PM sử dụng xạ rịng trung bình ngày chiết xuất từ ảnh LANDSAT Rnd_VT – công thức 2.31), ngày 04/11/2000 Rn-G (MJ/m2/d) Cao độ (m) Uz (m/s) Tmax (oC) eT (kPa) eT max (kPa) es (kPa) RH (%) eo(T min) (kPa) ea (kPa) U2 (m/s) 11.8704 2.9097276 20.55 13.8 27.3 1.57591 3.6238213 2.599864569 70 2.419 1.693 10.2957 2.9097276 20.65 13.3 28 1.52542 3.7749796 2.650198368 72 2.434 1.753 11.5005 17 3.6371595 21.25 14.8 27.7 1.68131 3.7095384 2.695422474 76 2.525 1.919 13.921 685 4.3645914 17.2 12.2 22.2 1.41931 2.6728283 2.046067573 54 1.962 1.06 11.885 36 2.9097276 20.9 13.2 28.6 1.51549 3.9088829 2.712187044 73 2.472 1.804 12.018 138 2.9097276 19.95 11.1 28.8 1.31973 3.9544265 2.637080598 79 2.331 1.842 P (kPa) Gama γ (kPa/ oC) 0.90650 101.3 0.89769 TT B (oC) Tmin (oC) Lamda λ (MJ kg-1) tu_delta MauSo_delta 900/(T+273) (o C ) Tu_ETo Mau_ETo ETo (mm/ngày) 0.067 Delta Δ (kPa/ o C) 0.149 2.436 9913.4 66486.62 3.06591722 1.26670847 0.283066 4.47 101.2 0.067 0.15 2.444 9974.7 66538.2 3.064873148 1.16856554 0.283808 4.12 0.77608 101.1 0.067 0.155 2.446 10349 66848.1 3.058623619 1.30688666 0.305189 4.28 0.98636 93.46 0.062 0.124 2.444 8042 64770.25 3.101309442 1.53501825 0.278544 5.51 0.90779 100.9 0.067 0.152 2.452 10129 66667.24 3.062266077 1.27939151 0.285385 4.48 0.79555 99.68 0.066 0.144 2.444 9552.7 66177.56 3.072196621 1.17919193 0.276214 4.27 es-ea (kPa) -4Phụ lục 2.1: Tính ETo theo FAO-56 PM sử dụng Rnd_FAO1 (trong xạ Mặt Trời tới Rs xác định theo Ăngstrom – công thức 2.32b), ngày 23/11/2001 e Tmin (kPa) eTmax(kPa) es (kPa) RH (%) eo(Tmin) (kPa) ea (kPa) 25.1 1.49777 3.186695 2.34223 72 2.2042 1.6864 12.4 25.7 1.43998 3.302286 2.37113 73 2.2042 1.7309 19.1 13.1 25.1 1.50759 3.186695 2.34716 76 2.2111 1.7838 5.09202326 14.05 10 18.1 1.22796 2.077002 1.65248 59 1.6037 0.9749 2.90972758 18.35 11.5 25.2 1.35698 3.205712 2.28134 72 2.1098 1.6425 2.18229568 16.05 8.5 23.6 1.10985 2.913023 2.01143 80 1.8241 1.6091 Rn-G (MJ/m2/d) Nắng thực (giờ) Cao độ (m) Uz (m/s) U2 (m/s) TTB ( oC) Tmin (oC) Tmax( oC) 9.4713347 9.3 2.90972758 19.05 13 9.4652334 9.2 2.18229568 19.05 9.3352376 8.8 17 2.18229568 8.7365243 9.6 685 8.9846555 8.6 36 8.5603445 7.9 138 es-ea (kPa) P (kPa) Gama γ (kPa/ oC) Lamda λ (MJ kg-1) tu_delta MauSo_delta 900/(T+273) (oC) Tu_ETo Mau_ETo ETo (mm/ngày) 0.0673 Delta Δ (kPa/ o C) 0.1375 0.6558 101.2645 2.4361 9033.0456 65715.3225 3.081664 0.927185 0.271419 3.42 0.6402 01.2173 0.0673 0.1375 2.4443 9033.0456 65715.3225 3.081664 0.820633 0.254709 3.22 0.5633 101.0992 0.0672 0.1378 2.4462 9061.2502 65740.9600 3.081137 0.779625 0.254948 3.06 0.6775 93.4607 0.0622 0.1040 2.4436 6572.3693 63176.8225 3.135342 1.043094 0.273785 3.81 0.6387 100.8752 0.0671 0.1323 2.4517 8646.1689 65356.9225 3.089068 0.870101 0.265738 3.27 0.4022 99.6794 0.0663 0.1165 2.4443 7475.1692 64186.2225 3.113648 0.587948 0.231931 2.54 -5Phụ lục 2.2: Tính ETo theo FAO-56 PM sử dụng Rnd_FAO2 (trong xạ Mặt Trời tới Rs xác định theo công thức thực nghiệm Việt Nam – công thức 2.32c), ngày 23/11/2001 Nắng Cao độ Rn-G (MJ/m2/d) thực (giờ) (m) Uz (m/s) U2 (m/s) TTB (oC) Tmin (oC) Tmax(oC) eTmax e Tmin (kPa) (kPa) es (kPa) RH (%) eo(Tmin) (kPa) ea (kPa) 13.2164 9.3 2.90972758 19.05 13 25.1 1.4978 3.1867 2.3422 72.000 2.2043 1.6864 13.4372 9.2 2.18229568 19.05 12.4 25.7 1.4400 3.3023 2.3711 73.000 2.2043 1.7309 14.0929 8.8 17 2.18229568 19.1 13.1 25.1 1.5076 3.1867 2.3471 76.000 2.2111 1.7838 11.4986 9.6 685 5.09202326 14.05 10 18.1 1.2280 2.0770 1.6525 59.000 1.6038 0.9750 13.9524 8.6 36 2.90972758 18.35 11.5 25.2 1.3570 3.2057 2.2813 72.000 2.1099 1.6426 14.9748 7.9 138 2.18229568 16.05 8.5 23.6 1.1099 2.9130 2.0114 80.000 1.8241 1.6092 es-ea (kPa) P (kPa) Gama γ (kPa/ oC) Lamda λ (MJ kg-1) tu_delta MauSo_delta 900/(T+273) (o C ) Tu_ETo Mau_ETo ETo (mm/ngày) 0.0673 Delta Δ (kPa/ o C) 0.1375 0.6558 101.2645 2.4361 9033.0456 65715.3225 3.0817 1.1372 0.2714 4.19 0.6402 101.2173 0.0673 0.1375 2.4443 9033.0456 65715.3225 3.0817 1.0434 0.2547 4.10 0.5633 101.0992 0.0672 0.1378 2.4462 9061.2502 65740.9600 3.0811 1.0472 0.2549 4.11 0.6775 93.4607 0.0622 0.1040 2.4436 6572.3693 63176.8225 3.1353 1.1603 0.2738 4.24 0.6388 100.8752 0.0671 0.1323 2.4517 8646.1689 65356.9225 3.0891 1.1382 0.2657 4.28 0.4023 99.6794 0.0663 0.1165 2.4443 7475.1692 64186.2225 3.1136 0.8927 0.2319 3.85 -6- Phụ lục 2.3: Tính ETo theo FAO-56 PM sử dụng xạ rịng trung bình ngày chiết xuất từ ảnh LANDSAT Rnd_VT – công thức 2.31), ngày 23/11/2001 Nắng Cao độ Rn-G (MJ/m2/d) thực (giờ) (m) Uz (m/s) U2 (m/s) TTB (oC) Tmin (oC) Tmax(oC) eTmax e Tmin (kPa) (kPa) es (kPa) RH (%) eo(Tmin) (kPa) ea (kPa) 10.429 9.3 2.90972758 19.05 13 25.1 1.4978 3.1867 2.3422 72.00 2.2043 1.6864 9.889 9.2 2.18229568 19.05 12.4 25.7 1.4400 3.3023 2.3711 73.00 2.2043 1.7309 10.639 8.8 17 2.18229568 19.1 13.1 25.1 1.5076 3.1867 2.3471 76.00 2.2111 1.7838 12.290 9.6 685 5.09202326 14.05 10 18.1 1.2280 2.0770 1.6525 59.00 1.6038 0.9750 11.019 8.6 36 2.90972758 18.35 11.5 25.2 1.3570 3.2057 2.2813 72.00 2.1099 1.6426 10.744 7.9 138 2.18229568 16.05 8.5 23.6 1.1099 2.9130 2.0114 80.00 1.8241 1.6092 es-ea (kPa) P (kPa) Gama γ (kPa/ oC) Lamda λ (MJ kg-1) tu_delta MauSo_delta 900/(T+273) (o C ) Tu_ETo Mau_ETo ETo (mm/ngày) 0.0673 Delta Δ (kPa/ o C) 0.1375 0.6558 101.2645 2.4361 9033.0456 65715.3225 3.081664 0.980882 0.271419 3.61 0.6402 101.2173 0.0673 0.1375 2.4443 9033.0456 65715.3225 3.081664 0.844399 0.254709 3.32 0.5633 101.0992 0.0672 0.1378 2.4462 9061.2502 65740.9600 3.081137 0.852937 0.254948 3.35 0.6775 93.4607 0.0622 0.1040 2.4436 6572.3693 63176.8225 3.135342 1.193921 0.273785 4.36 0.6388 100.8752 0.0671 0.1323 2.4517 8646.1689 65356.9225 3.089068 0.979883 0.265738 3.69 0.4023 99.6794 0.0663 0.1165 2.4443 7475.1692 64186.2225 3.113648 0.691687 0.231931 2.98