Ví dụ: 7h số đọc thang bốc hơi là 1,2 mm; 19h số đọc thang bốc hơi là 5,7 mm. Khi đó, lượng bốc hơi từ 7h đến 19h là: 5,7 – 1,2 = 4,5 mm. Chú ý, khi đọc trị số bốc hơi thì đọc số tại vị trí ngang mặt lõm của mực nước trong ống Pichẹ
Cách thay giấy thấm: Hàng ngày sau quan trắc 7h phải đổ thêm nước và thay
giấy thấm ống Pichẹ Sau khi đổ thêm nước, chờ cho nước thấm hết giấy rồi đọc trị số, ghi bên cạnh số đọc lúc 7h để tính lượng bốc hơi lúc 19h. Phải giữ ống bốc hơi Piche sạch sẽ, không cáu bẩn, dùng nước mưa, nước sạch để đổ vào ống (Bộ Tài nguyên Môi trường 2012)
1.2.1.2. Phương pháp thủy tiêu kế (Lysimeter)
Bản chất của phương pháp thủy tiêu kế được mô tả như sau:
Thủy tiêu kế (Lysimeter) là một thiết bị dùng để xác định giá trị bốc thoát hơi tham chiếu (ET0) của một cây trồng theo một điều kiện tưới chủ động. Bằng cách đo thể tích nước hay trọng lượng ta có thể xác định lượng bốc thoát hơi dựa vào phương trình cân bằng nước.
(Nguồn: Lê Anh Tuấn 2009) Hình 1.4. Bố trí thiết bị đo (ET0) theo phương pháp thủy tiêu kế
Thủy tiêu kế có dạng là một thùng hình trụ tròn có đường kính khoảng 30 cm và độ cao 25cm được đổ đầy đất như loại đất canh tác. Đáy thùng có chỗ để nước thoát ra nhằm đo lượng nước thấm sâụ Hầm chứa lượng nước thấm sâu có đường kính 30cm và chiều cao 42cm được bố trí ngay gần đó. Bên cạnh đó, thiết bị đo mưa bằng thùng đo mưa cũng được lắp đặt (Thiết bị đo mưa có thể chứa được 5,5 lít nước). Mặt trên của thùng, cây trồng được gieo cấy đều đặn giống như môi trường bên ngoài (Hình 1.4).
Một cách tổng quát, bằng cách đo lượng mưa rơi trong khu vực (R), lượng tưới (I) và lượng thấm sâu xuống đất (P), lượng bốc thoát hơi (ET0) sẽ được xác định theo:
ET0 = R + I – P (1.1) Một số nơi, người ta dùng cân (đặt ở dưới thủy tiêu kế) để xác định sự thay
đổi lượng nước ở thủy tiêu kế để xác định lượng bốc thoát hơi nước (Lê Anh Tuấn, 2009).
Ưu nhược điểm của phương pháp:
+ Ưu điểm: Thiết bị đo cho độ chính xác xác định ET0 là tương đối khá cao, dễ sử dụng.
+ Nhược điểm: Bố trí các điểm đo phức tạp, kinh phí xây dựng và duy trì bảo quản thiết bị đo lớn.
Đo mưa (R) Bốc thoát hơi
(ET ) Tưới
(I)
Hầm đo thấm sâu DP
1.2.2. Các mô hình sử dụng năng lượng bức xạ mặt trời (radiaton-based models)
Có nhiều mô hình khác nhau được sử dụng để ước tính lượng bốc thoát hơi nước từ bề mặt lớp phủ. Trong đó, các mô hình ước tính lượng bốc thoát hơi nước ET sử dụng năng lượng bức xạ mặt trời được các nhà khoa học đặc biệt quan tâm. Các mô hình tiêu biểu sử dụng năng lượng bức xạ mặt trời để tính lượng bốc thoát hơi nước bề mặt như sau:
1.2.2.1. Mô hình Makkink
Năm 1957 mô hình Makkink đã được đề xuất và sử dụng rộng rãi cho đến ngày nay do đã giảm bớt một số các thông số đo ngoại nghiệp mà phương pháp FAO Penman đang sử dụng. Độ chính xác của mô hình phụ thuộc vào hệ số tuyến tính u, v của phương trình, đòi hỏi phải có kết quả kiểm chứng thực địạ Mô hình Makkink đề xuất công thức tính toán lượng bốc thoát hơi nước từ năng lượng bức xạ mặt trời như sau: 𝐸𝑇𝑎 = 𝑢 ∆ ∆ + 𝛾 𝑅𝑆 𝜆 + 𝑣 (1.2)
Trong đó: ETa- là lượng bốc thoát hơi nước thực tế từ bề mặt lớp phủ (mm/ngày); Rs - bức xạ mặt trời đi tới bề mặt đất (MJ/m2/ngày); Δ - Độ dốc đường cong áp suất hơi nước bão hòa của không khí (kPa/°C), γ - Hằng số Psychrometric (kPa/°C), λ - Giá trị nhiệt ẩn của quá trình bốc thoát hơi nước (MJ/kg); u, v - hệ số tuyến tính của mô hình Makkink.
Hệ số tuyến tính của mô hình Makkink sử dụng trong việc tính lượng bốc thoát hơi nước từ bề mặt địa hình được xác định vào thời điểm đề xuất mô hình (1957) có trị số là: u = 0,61 và v = 0,12. Tuy nhiên, hệ số tuyến tính u, v phụ thuộc vào điều kiện khí hậu và yếu tố địa hình của từng khu vực nghiên cứụ Theo nghiên cứu của Hasen xác định tại Hà Lan năm 1984, hệ số u, v có giá trị là u = 0,70 và v = 0. Kết quả nghiên cứu giữa Trường Đại học Uppsala (Thụy Điển) và Trường Đại học Louisiana (Mỹ) do Xue và Singh thực hiện khảo sát vào năm 1999, hệ số u = 0,77 và v = 0,22 (Makkink 1957)
Năm 1972, Priestley - Taylor đề xuất mô hình tính lượng bốc thoát hơi nước từ năng lượng bức xạ mặt trời có dạng sau:
𝐸𝑇𝑎 = 𝑎 ∆
∆ − 𝛾 𝑅𝑛
𝜆 + 𝑏 (1.3)
Trong đó: ETa- là lượng bốc thoát hơi nước từ bề mặt lớp phủ (mm/ngày); Rn – Giá trị năng lượng bức xạ ròng mặt trời đi tới bề mặt đất (MJ/m2/ngày); Δ - Độ dốc đường cong áp suất hơi nước bão hòa của không khí (kPa/°C), γ - Hằng số Psychrometric (kPa/°C), λ - Giá trị nhiệt ẩn của quá trình bốc thoát hơi nước (MJ/kg); a, b là hệ số tuyến tính của mô hình Priestley – Taylor.
Hệ số tuyến tính a, b của mô hình Priestley - Taylor (1972) sử dụng trong việc tính lượng bốc thoát hơi nước từ bề mặt lớp phủ có trị số là: a = 0,61 và b = 0,12; kiểm định tại Thụy Sĩ năm 1984, cho kết quả a = 0,90 và b = 0; kiểm định tại châu Á (Đài Loan), 2005, a = 1,00 và b = 0.
Mô hình Priestley-Taylor là mô hình tính lượng bốc thoát hơi nước được dùng khá phổ biến trên thế giớị Ưu điểm của mô hình là cho độ chính xác tính toán lượng bốc thoát hơi nước cao, yêu cầu số liệu đầu vào ít, dễ tính toán. Độ chính xác tính lượng bốc thoát hơi nước phụ thuộc vào kết quả lựa chọn tham số a, b của mô hình. Hệ số tuyến tính a, b phụ thuộc vào điều kiện địa hình, khí hậu của từng khu vực nghiên cứu cụ thể, cần có các nghiên cứu thực nghiệm để xác định (Priestley và Taylor 1972)
1.2.2.3. Mô hình Irmak
Irmak đã sử dụng dữ liệu trong nhiều năm tại Mỹ để nghiên cứu thực nghiệm, năm 2003 đã đề xuất mô hình tính lượng bốc thoát hơi nước có dạng như sau:
ETa = 0,289Rnd + 0,023T + 0,489 (1.4)
Trong đó ETa – Lượng bốc thoát hơi nước thực tế trung bình ngày (mm/ngày), Rnd - Giá trị năng lượng bức xạ ròng trung bình ngày (MJ/m2/ngày) T - Giá trị nhiệt độ trung bình ngày (0C)
Irmak là mô hình đơn giản để tính lượng bốc thoát hơi nước thực tế trung bình ngàỵ Kết quả tính lượng bốc thoát hơi nước bề mặt lớp phủ theo mô hình Irmak có
độ chính xác khá cao, sử dụng rất ít tham số đầu vào và thường được sử dụng trong thực tế (Irmak và cộng sự 2003)
1.2.2.4. Kết quả nghiên cứu tính lượng bốc thoát hơi nước từ một số mô hình sử dụng năng lượng bức xạ mặt trời
Theo nghiên cứu của (Xu1 và Singh 2001) thực nghiệm sử dụng 4 mô hình tính lượng bốc thoát hơi nước thực tế từ năng lượng mặt trời, so sánh với lượng bốc thoát hơi nước đo trực tiếp từ thiết bị đo là chậu Pan. Kết quả và sai số tính toán được thể hiện ở bảng sau:
Bảng 1.1. Bảng thống kê sai số của một số mô hình sử dụng bức xạ mặt trời tính toán lượng bốc thoát hơi nươc ET0 so với thiết bị đo trực tiếp chậu Pan
Mùa Epan (mm) EHag (mm) Sai số ETa (%) EMak (mm) Sai số ETa (%) EPrt (mm) Sai số ETa (%) ETur (mm) Sai số ETa (%) Đông 1,04 0,53 -49,04 0,83 -20,19 0,93 -10,6 0,35 -66,34 Xuân 2,80 2,58 -7,85 2,93 4,64 2,77 -1,07 2,51 -10,35 Hè 4,37 4,43 1,37 4,51 3,20 4,60 5,26 4,42 1,14 Thu 1,71 1,39 -18,71 1,64 -4,09 1,60 -6,43 1,50 -12,28 Trung bình năm 2,48 2,23 -10,00 2,48 0,00 2,48 0,00 2,19 -11,69
Chú thích: Kí hiệu các phương pháp: Hag = Hargreaves, Mak = Makkink, Prt = Priestley -Taylor và Tur = Turn.
Qua kết quả nghiên cứu tác giả đã kết luận sử dụng các mô hình Makkink và Priestley - Taylor ước tính lượng bốc thoát hơi nước sử dụng năng lượng mặt trời cho kết quả với độ tin cậy cao nhất.
1.2.3. Các mô hình kết hợp (combined models)
Các mô hình tính toán lượng bốc thoát hơi nước sử dụng phương pháp kết hợp dữ liệu được phát triển bởi công thức Penman (1948).
1.2.3.1. Mô hình Penman
Năm 1948, Penman đã đề xuất mô hình tính toán lượng bốc hơi nước bằng phương pháp kết hợp giữa phương pháp cân bằng năng lượng bức xạ mặt trời và phương pháp chuyển đổi khối lượng có dạng như sau:
𝐸𝑇0 = ∆
∆ + 𝛾(𝑅𝑛− 𝐺) + 𝐾𝑤∆ + 𝛾𝛾 (𝑎𝑤 + 𝑏𝑤𝑢2)(𝑒𝑠− 𝑒𝑎) 𝜆
(1.5) Trong đó: ET0 - Lượng bốc thoát hơi nước (mm/ngày); λ - Giá trị nhiệt ẩn của quá trình bốc thoát hơi nước (MJ/kg); Δ - Độ dốc đường cong áp suất hơi nước bão hòa của không khí (kPa/°C); γ - Hằng số Psychrometric (kPa/°C); G - Mật độ dòng nhiệt trong đất (MJ/ m2/ngày); Rn - Giá trị năng lượng bức xạ ròng (MJ/m2/ngày); u – Tốc độ gió ở độ cao 2m (m/s); 𝑒𝑠 - Áp suất hơi nước bão hòa ở nhiệt độ trên bề mặt nước (kPa); 𝑒𝑎 - Áp suất hơi nước bão hòa ở độ cao đo trên mặt nước (kPa); 𝑎𝑤, 𝑏𝑤- hàm số của gió; 𝐾𝑤- hằng số.
Ưu điểm của mô hình là ước tính lượng bốc thoát hơi nước tham chiếu với độ chính xác caọ Tuy nhiên, mô hình sử dụng nhiều tham số đầu vào, đòi hỏi nhiều số liệu về khí tượng, các tham số phải tính toán qua nhiều bước trung gian, phức tạp (Penman 1948).
1.2.3.2. Mô hình Penman-Monteith
Năm 1965, Monteith đã phát triển mô hình của Penman 1948 thành mô hình tính lượng bốc thoát hơi nước áp dụng cho các khu vực có diện tích lớn như sau:
𝜆𝐸𝑇0 = [∆(𝑅𝑛− 𝐺)] + [𝜌𝑎𝐶𝑝(𝑒𝑟𝑠− 𝑒𝑎) 𝑎 ] ∆ + 𝛾(1 +𝑟𝑟𝑠 𝑎) (1.6)
Trong đó: ET0 - Lượng bốc thoát hơi nước (mm/ngày); λ - Giá trị nhiệt ẩn của quá trình bốc thoát hơi nước (MJ/kg); Δ - Độ dốc đường cong áp suất hơi nước bão hòa của không khí (kPa/°C); γ - Hằng số Psychrometric (kPa/°C); G - Mật độ dòng nhiệt trong đất (MJ/ m2/ngày); Rn - Giá trị năng lượng bức xạ ròng (MJ/m2/ngày); Cp - là nhiệt lượng riêng của không khí (MJ/kg/0C); ρa - là mật độ không khí trung bình ở áp suất không đổi (kg/m3); rs và ra - là sức cản bề mặt và khí động học (s/m); 𝑒𝑠 - Áp suất hơi nước bão hòa ở nhiệt độ trên bề mặt nước (kPa); 𝑒𝑎 - Áp suất hơi nước bão hòa ở độ cao đo trên mặt nước (kPa).
1.2.3.3. Mô hình FAO 56 Penman – Monteith
Hầu hết các tài liệu nghiên cứu về hiện tượng bốc thoát hơi nước bề mặt lớp phủ đều tiếp cận tới mô hình FAO 56 Penman – Monteith để tính lượng bốc thoát hơi
nước tham chiếu (ET0). Mô hình FAO 56 Penman – Monteith nằm trong nhóm mô hình kết hợp và là một trong những mô hình chính xác nhất để xác định lượng bốc thoát hơi nước trong các thời gian khác nhaụ Mô hình đã được Tổ chức Lương thực và Nông nghiệp thế giới (FAO) đề xuất năm (1990), và cụ thể thành tiêu chuẩn áp dụng trên toàn thế giớị Phương trình FAO 56 Penman – Monteith dùng để xác định giá trị bốc thoát hơi nước tham chiếu, là một hàm số phụ thuộc nhiều dữ liệu thời tiết tại chỗ và xung quanh khu vực xem xét. Các thông số này được mô tả chi tiết trong tài liệu hướng dẫn tính toán của FAO (Allen và cộng sự 1998).
Công thức tính lượng bốc thoát hơi tham chiếu ET0 theo FAO 56 Penman – Monteith như sau:
𝐸𝑇0 =0,48∆(𝑅𝑛 − 𝐺) + 𝛾 900 𝑇 + 273𝑢2(𝑒𝑠− 𝑒𝑎) ∆ + 𝛾(1 + 0,3𝑢2) (1.7) Trong đó:
ET0 - lượng bốc thoát hơi tham chiếu chung đối với cây trồng (mm/ngày); Rn - Giá trị năng lượng bức xạ ròng trên bề mặt cây trồng (MJ/m2/ngày); G - Mật độ dòng nhiệt trong đất (MJ/m2/ngày); T - Nhiệt độ trung bình ngày tại vị trí 2m từ mặt đất (°C); u2 - Tốc độ gió tại chiều cao 2 m từ mặt đất (m/s); es - Áp suất hơi nước bão hòa (kPa); ea - Áp suất hơi nước thực tế (kPa); Δ - Độ dốc đường cong áp suất hơi nước bão hòa của không khí (kPa/°C); γ - Hằng số Psychrometric (kPa/°C).
Các thông số trong công thức trên có thể tính từ số liệu do ngành khí tượng cung cấp kết hợp với công thức và bảng tra theo tài liệu của FAO, Granier (1985). Mô hình FAO 56 Penman - Monteith cho kết tính toán với độ chính xác rất cao nhưng khối lượng tính toán lớn, phức tạp và phải có đủ tài liệu ban đầụ
1.3. Các mô hình xác ước tính lượng bốc thoát hơi nước từ dữ liệu ảnh vệ tinh
Các mô hình ước tính lượng bốc thoát hơi nước ứng dụng dữ liệu ảnh vệ tinh hiện nay trên thế giới đang được nghiên cứu và phát triển theo 2 hướng chính:
1. Cân bằng năng lượng bề mặt đất EB (Energy Balance), đây là mô hình sử dụng phản xạ, bức xạ bề mặt lớp phủ trên ảnh vệ tinh trong dải phổ nhìn thấy và cận hồng ngoại của phổ điện từ, nhiệt độ bề mặt từ kênh ảnh nhiệt hồng ngoạị
2. Phản xạ dựa trên hệ số cây trồng (Kc) và phương pháp xác định lượng bốc thoát hơi nước tham chiếu ET0, trong đó hệ số cây trồng Kc liên quan đến chỉ số thực vật bắt nguồn từ giá trị phản xạ tán lá.
Trong hai phương pháp này thì phương pháp tính toán sử dụng mô hình cân bằng năng lượng bề mặt được các nhà khoa học trên thế giới đặc biệt quan tâm và ứng dụng chủ yếu trong các nghiên cứụ Các ứng dụng dựa trên nguyên lý cân bằng năng lượng bề mặt gồm các mô hình tiêu biểu như sau:
1.3.1. Mô hình cân bằng năng lượng bề mặt đất SEBAL (Surface Energy Balance Algorithms for Land)
Nguyên lý của mô hình SEBAL (Bastiaanssen và cộng sự 1998) là sử dụng phương trình cân bằng năng lượng bề mặt đất. Lượng bốc thoát hơi nước ET được tính toán từ ảnh vệ tinh và dữ liệu khí tượng. Ảnh vệ tinh cung cấp thông tin trong một khoảng thời gian tức thời, mô hình SEBAL tính lượng bốc thoát hơi nước tức thời tại thời điểm chụp ảnh. Bốc thoát hơi nước ET được tính toán cho từng điểm ảnh theo phương trình cân bằng năng lượng lượng bề mặt như sau:
𝜆𝐿𝐸 = 𝑅𝑛− 𝐺 − 𝐻 (1.8)
Trong đó: Rn - Giá trị năng lượng bức xạ ròng mặt trời (W/m2/giờ); 𝜆𝐿𝐸- Năng lượng dòng nhiệt ẩn (W/m2/giờ), dòng nhiệt gây ra hiện tượng bốc thoát hơi nước; G - Thông lượng nhiệt mà mặt đất hấp thụ (W/m2/giờ); H - Năng lượng nhiệt cảm ứng (W/m2/giờ).
(Nguồn: Allen và cộng sự 2002) Hình 1.5. Nguyên lý cân bằng năng lượng bề mặt SEBAL
Giá trị năng lượng bức xạ ròng mặt trời (Rn) đại diện cho năng lượng bức xạ thực tế tại bề mặt đất. Nó là hiệu số giữa bức xạ mặt trời đi tới bề mặt đất và lượng bức xạ phản xạ, phát xạ từ bề mặt đất trở về khí quyển. Bức xạ ròng hấp thụ bởi mặt đất được mô tả bằng biểu thức sau:
Rn = RS↓ - α RS↓ + RL↓ - RL↑ - (1-εo) RL↓ (1.9) Trong đó: RS↓ - Năng lượng tới sóng ngắn (W/m2/giờ); α RS↓ - Năng lượng phản xạ sóng ngắn (W/m2/giờ); RL↓ - Năng lượng tới sóng dài (W/m2/giờ); RL↑ - Năng lượng phát xạ sóng dài (W/m2/giờ); (1-εo)RL↓ - Năng lượng phản xạ sóng dài(W/m2/giờ).
Trong công thức (1.9) tính bức xạ ròng bề mặt đất hấp thụ, tia tới sóng ngắn RS↓ còn lưu lại trên mặt đất là một hàm số của suất phân sai bề mặt đất α. Suất phân sai bề mặt là một hệ số phản xạ được định nghĩa là tỷ số giữa năng lượng bức xạ, phản xạ đến các năng lượng bức xạ liên quan trong quang phổ mặt trờị Nó được tính toán bằng cách sử dụng thông tin hình ảnh vệ tinh về quang phổ cho mỗi kênh ảnh của vệ tinh. Tia tới sóng ngắn R S↓ được tính toán sử dụng hằng số mặt trời, góc tới năng lượng mặt trời, khoảng cách tương đối giữa Mặt trời và Trái đất, tính toán truyền
CÂN BẰNG NĂNG LƯỢNG BỐC THOÁT HƠI NƯỚC ET
ET = Rn – G – H
G (Thông lượng nhiệt mặt đất hấp thụ)
Rn (Bức xạ ròng mặt trời) ET
(Bốc thoát hơi nước) H (Thông lượng nhiệt
cảm ứng)
dẫn qua khí quyển. Tia tới sóng dài RL↓ được tính toán bằng sử dụng phương trình Stefan-Boltzmann với sự truyền dẫn qua khí quyển và nhiệt độ tham chiếu bề mặt được chọn. Tia phát xạ sóng dài RL↑ được tính toán sử dụng phương trình Stefan- Boltzmann với tính toán sự phát xạ bề mặt và nhiệt độ bề mặt. Nhiệt độ bề mặt được tính toán từ thông tin ảnh vệ tinh trên kênh bức xạ nhiệt. Độ phát xạ bề mặt là tỷ số giữa bức xạ phát xạ thực tế từ bề mặt phát xạ bởi một vật đen ở nhiệt độ bề mặt tương ứng. Trong mô hình SEBAL, phát xạ được tính toán từ hàm chỉ số thực vật. Cuối