Đánh giá tác động của biến đổi khí hậu đến dòng chảy trạm đồng trăng, sông cái nha trang, tỉnh khánh hoà

Ứng dụng mô hình SWAT đánh giá tác động của biến đổi khí hậuđến dòng chảy mùa lũ trạm Đồng Trăng sông Cái Nha Trang ...74 3.5.. Để khảo sát bài toán biến đổi dòng chảy trong tương lai ở

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ

ĐẠI HỌC QUỐC GIA HÀ NỘI

TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN

LUẬN VĂN THẠC SĨ KHOA HỌC

NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS TRẦN NGỌC ANH

Hà Nội – 2017

MỤC LỤC

DANH MỤC HÌNH 1

DANH MỤC BẢNG 6

DANH MỤC CÁC KÝ HIỆU VIẾT TẮT 7

MỞ ĐẦU 9

CHƯƠNG 1 - TỔNG QUAN 12

1.1 Đặc điểm địa lý tự nhiên 12

1.1.1 Vị trí địa lý 12

1.1.2 Địa hình, địa mạo 12

1.1.3 Địa chất, thổ nhưỡng 13

1.1.4 Thảm thực vật 15

1.1.5 Khí hậu 16

1.1.6 Thủy văn – Hải văn 17

1.1.7 Mạng lưới trạm khí tượng thủy văn trên lưu vực sông 21

1.2 Đặc điểm kinh tế- xã hội 21

1.2.1 Tình hình kinh tế 21

1.2.2 Đặc điểm xã hội 23

1.3 Tổng quan về biến đổi khí hậu 23

1.3.1 Khái niệm 23

1.3.2 Biểu hiện về BĐKH 23

1.3.3 Kịch bản biến đổi khí hậu cho khu vực nghiên cứu 29

1.4 Một số nghiên cứu tương tự và trên lưu vực sông Cái Nha Trang tỉnh Khánh Hoà 30

CHƯƠNG 2 - GIỚI THIỆU MÔ HÌNH SWAT 33

2.1 Xuất xứ của mô hình SWAT 33

2.2 Cấu trúc của mô hình 34

2.2.1 Pha đất của chu trình thuỷ văn 34

2.2.2 Pha diễn toán của chu trình thuỷ văn 35

2.3 Các quá trình thành phần trong mô hình SWAT 35

2.3.1 Quá trình hình thành dòng chảy mặt 35

2.3.2 Quá trình hình thành dòng chảy ngầm 48

2.3.3 Quá trình diễn toán dòng chảy trong sông 50

2.4 Các số liệu đầu vào và kết quả của mô hình 56

2.4.1 Các số liệu đầu vào của mô hình 56

2.4.2 Kết quả của mô hình 57

2.5 Các thông số của mô hình 57

2.5.1 Các thông số tính quá trình hình thành dòng chảy mặt 57

2.5.2 Các thông số tính toán dòng chảy ngầm 58

2.5.3 Các thông số diễn toán dòng chảy trong kênh 58

2.6 Đánh giá mô hình 58

2.7 Tiến trình mô phỏng SWAT 59

CHƯƠNG 3 - ỨNG DỤNG MÔ HÌNH SWAT ĐÁNH GIÁ TÁC ĐỘNG CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN DÒNG CHẢY LŨ TRẠM ĐỒNG TRĂNG, SÔNG CÁI NHA TRANG, TỈNH KHÁNH HÒA 61

3.1 Cơ sở dữ liệu 61

3.1.1 Số liệu mặt đệm 61

3.1.2 Số liệu khí tượng thủy văn 61

3.2 Thiết lập mô hình SWAT 62

3.3 Hiệu chỉnh, kiểm định mô hình 69

3.3.1 Đánh giá mô hình 70

3.3.2 Các thông số mô hình 70

3.3.3 Kết quả hiệu chỉnh 71

3.3.4 Kiểm định mô hình 72

2

3.4 Ứng dụng mô hình SWAT đánh giá tác động của biến đổi khí hậu

đến dòng chảy mùa lũ trạm Đồng Trăng sông Cái Nha Trang 74

3.5 Đánh giá tác động của biến đổi khí hậu đến giá trị đỉnh lũ trạm Đồng Trăng, sông Cái Nha Trang, tỉnh Khánh Hoà 78

3.5.1 Hiệu chỉnh 83

3.5.2 Kiểm định 84

3.5.3 Ứng dụng mô hình SWAT đánh giá tác động của biến đổi khí hậu đến giá trị đỉnh lũ trạm Đồng Trăng sông Cái Nha Trang 87

3.6 Thảo luận 92

KẾT LUẬN VÀ KIẾN NGHỊ 94

TÀI LIỆU THAM KHẢO 97

DANH MỤC HÌNH

Hình 1.1 Sơ đồ lưu vực sông Cái Nha Trang

Hình 1.2 Đường mực nước trận lũ lịch sử năm 2009

Hình 1.3 Đường mực nước trận lũ lớn năm 1998

Hình 1.4 Đường mực nước trận lũ lớn năm 2010

Hình 1.5.Chuẩn sai nhiệt độ trung bình toàn cầu thời kỳ 1850-2012 (so với thời kỳ 1961-1990)

Hình 1.6 Chuẩn sai nhiệt độ trung bình toàn cầu (o C) thời kỳ 1950-2015

Hình 1.7 Biến đổi của nhiệt độ trung bình năm thời kỳ 1901-2012

Hình 1.8 Biến đổi của lượng mưa năm thời kỳ 1901-2010 và thời kỳ 1951-2010

Hình 1.9 Chuẩn sai nhiệt độ (o C) trung bình năm (a) và nhiều năm (b) trên quy mô cả nước

Hình 1.10 Chuẩn sai nhiệt độ trung bình năm ( C) đối với các trạm ven biển vào hải đảo

Hình 1.11 Thay đổi nhiệt độ trung bình năm ( C) thời kỳ 1958-2014 o Hình 1.12 Thay đổi lượng mưa năm (%) thời kỳ 1958-2014

Hình 2.1 Sự khác nhau giữa phân phối độ ẩm theo chiều sâu mô phỏng theo phương trình Green và Ampt và trong thực tế

Hình 2.2 Hình dạng kênh chính hoặc kênh nhánh

Hình 2.3 Mô phỏng đoạn sông theo phương pháp Muskingum

Hình 2.4 Tiến trình mô phỏng của SWAT

Hình 3.1 Tiến trình chạy trong SWAT

Hình 3.2 Bản đồ DEM lưu vực sông Cái Nha Trang, trạm Đồng Trăng

Hình 3.3 Phân chia tiểu lưu vực khu vực nghiên cứu

Hình 3.4 Bản đồ thổ nhưỡng lưu vực sông Cái Nha Trang

Hình 3.5 Bản đồ thảm phủ lưu vực sông Cái Nha Trang

Hình 3.6 Hình ảnh kết quả mô phỏng dòng chảy mùa lũ

Hình 3.7 Dòng chảy bình quân ngày tính toán và thực đo trạm Đồng Trăng 71 Hình 3.8 Dòng chảy bình quân ngày tính toán và thực đo trạm Đồng Trăng 72 Hình 3.9 % Thay đổi lượng mưa theo mùa theo kịch bản RCP 4.5

Hình 3.10 % Thay đổi lượng mưa theo mùa theo kịch bản RCP 8.5

Hình 3.11 Thay đổi dòng chảy trung bình các tháng mùa lũ vào giữa thế kỷ, cuối thế kỷ so với giai đoạn nền theo kịch bản RCP 4.5

Hình 3.12 Thay đổi dòng chảy trung bình các tháng mùa lũ vào giữa thế kỷ, cuối thế kỷ so với giai đoạn nền theo kịch bản RCP 8.5

Hình 3.13 Dữ liệu mưa đầu vào thời đoạn tính toán ngắn hơn (theo giờ, phút, giây) trong mô hình SWAT

4

Hình 3.14 Cài đặt bước thời gian tính toán thời đoạn ngắn hơn (theo giờ,

phút, giây) trong mô hình SWAT 80

Hình 3.15 File cio được sửa file đọc kết quả 81

Hình 3.16 File.cio được sửa theo mưa đầu vào thời đoạn ngắn hơn 82

Hình 3.17 File bsn được sửa để đọc kết quả theo thời đoạn ngắn hơn 83

Hình 3.18 Đường quá trình lưu lượng tính toán và thực đo tại trạm Đồng Trăng của trận lũ từ 1 - 4/XI năm 2009 84

Hình 3.19 Đường quá trình lưu lượng tính toán và thực đo tại trạm Đồng Trăng của trận lũ từ 11-14/XII/2016 85

Hình 3.20 % Thay đổi lượng mưa 5 ngày lớn nhất theo các kịch bản 87

Hình 3.21 So sánh mức độ phù hợp giữa đường tần suất lý luận và đường thực nghiệm của lượng mưa 5 ngày lớn nhất giai đoạn 1986-2005 89

Hình 3.22 Thay đổi dòng chảy cực trị giữa thế kỷ và cuối thế kỷ so với thời đoạn nền 1986-2005 theo kịch bản RCP 4.5 91

Hình 3.23 Thay đổi dòng chảy cực trị giữa thế kỷ và cuối thế kỷ so với thời đoạn nền 1986-2005 theo kịch bản RCP 8.5 92

DANH MỤC BẢNG

Bảng 1.1 Diện tích các loại đất đai tỉnh Khánh Hòa 14

Bảng 1.2 Cán cân nước các lưu vực[7] 18

Bảng 1.3 Các trạm đo KTTV trên lưu vực sông 21

Bảng 1.4 Thay đổi lượng mưa (%) trong 57 năm qua (1958-2014) 29

Bảng 1.5 % Thay đổi lượng mưa 5 ngày lớn nhất theo các kịch bản 29

Bảng 1.6 % thay đổi lượng mưa mùa theo các kịch bản 30

Bảng 3.1 Diện tích các tiểu lưu vực trên lưu vực sông Cái Nha Trang tính đến trạm Đồng Trăng 63

Bảng 3.2 Phân loại đất lưu vực sông Cái Nha Trang 66

Bảng 3.3 Phân loại các loại thảm phủ lưu vực sông Cái Nha Trang 67

Bảng 3.4 Số liệu mưa đầu vào mô hình SWAT 68

Bảng 3.5 Mức độ mô phỏng tương ứng với chỉ số Nash 70

Bảng 3.6 Kết quả bộ thông số khi hiệu chỉnh mô hình SWAT 73

Bảng 3.7 % Thay đổi lượng mưa theo mùa theo các kịch bản 74

Bảng 3.8 Kết quả hiệu chỉnh bộ thông số cho lưu vực sông Cái Nha Trang 85 Bảng 3.9 % Thay đổi lượng mưa 5 ngày lớn nhất theo các kịch bản 87

Bảng 3.10 Bảng số liệu lượng mưa 1 ngày, 3 ngày, 5 ngày lớn nhất 89

6

DANH MỤC CÁC KÝ HIỆU VIẾT TẮT

Digital Elevation Model (Mô hình số độ cao)Danish Hydraulic Institute (Viện nghiên cứu thủy lựcĐan Mạch)

Environmental Systems Research Institute (Viện nghiêncứu hệ thống môi trường)

Global Climate Model (Mô hình khí hậu toàn cầu)Hydrologic Engineering Center’s Hydrologic ModellingSystem (Hệ thống mô hình hóa thủy văn trung tâm kỹthuật thủy văn)

Soil Conservation ServiceSoil and Water Assessment Tool (Mô hình đánh giá môitrường lưu vực sông)

KH KTTV& MT Khoa học Khí tượng Thủy văn và Môi trường

WHO World Health Organization

LỜI CÁM ƠN

Luận văn được hoàn thành tại Khoa Khí tượng Thủy văn và Hải dương học,trường Đại học Khoa học Tự nhiên, Đại học Quốc Gia Hà Nội Trong quá trình thựchiện, tác giả xin bày tỏ lòng biết ơn sâu sắc đến các thầy cô trong khoa Khí tượngThủy văn và Hải dương học về sự hỗ trợ chuyên môn và kĩ thuật Tác giả cũng xinchân thành cám ơn TS Nguyễn Xuân Hiển - Giám đốc trung tâm Nghiên cứu Khítượng Thủy văn biển - Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu,cùng các anh chị em trong trung tâm đã tạo điều kiện và giúp đỡ Đặc biệt, xin bày

tỏ sự cám ơn chân thành nhất đến thầy hướng dẫn: PGS.TS Trần Ngọc Anh đãhướng dẫn tận tình và định hướng nghiên cứu để tác giả hoàn thành luận văn này

Học viên Nguyễn Thị Phương

8

MỞ ĐẦU

Sông Cái Nha Trang (còn có tên gọi là sông Phú Lộc, sông Cù) là sông lớnnhất tỉnh Khánh Hòa có chiều dài khoảng 79km, tuy ngắn nhưng giữ vai trò vô cùngquan trọng trong sự phát triển kinh tế - xã hội tỉnh Khánh Hòa Nhìn chung, dòngchảy trên lưu vực sông Cái Nha Trang chia thành 2 mùa rõ rệt:

Mùa lũ kéo dài 4 tháng, bắt đầu từ tháng IX và kết thúc vào tháng XII Lũcác sông thuộc lưu vực sông Cái Nha Trang tương tự các sông khác ở khu vực miềnTrung, khi có nước tập trung nhanh, lên xuống đột ngột, đường quá trình có dạngrăng cưa Nguyên nhân gây lũ như đã nêu chủ yếu do mưa có cường độ lớn, lũ docác hình thế thời tiết đơn độc gây ra thường là lũ thời kỳ đầu và cuối mùa, đườngquá trình lũ có đỉnh nhọn (tháng IX, tháng XII); lũ do tổ hợp nhiều hình thế thời tiếtphức tạp thường là thời kỳ giữa mùa lũ (tháng X – XI), đường quá trình lũ kéo dài

và có nhiều ngọn kế tiếp nhau Lượng dòng chảy 4 tháng mùa lũ chiếm từ 65% –66% lượng dòng chảy cả năm Tổng lượng dòng chảy (W) phụ thuộc vào độ sâudòng chảy (Y) và diện tích lưu vực F, có giá trị lớn nhất đối với sông Cái Nha Trang

2319 triệu m3

Mùa kiệt bắt đầu từ tháng I, kết thúc vào tháng VIII, kéo dài tới 8 tháng.Dòng chảy cạn chủ yếu là phần nước còn lại của mùa lũ năm trước sau đó giảmnhanh chóng theo đường nước rút và thường đạt trị số thấp nhất vào thời kỳ cuốitháng III và tháng IV Lượng dòng chảy tháng III, IV chỉ đạt 2,8% – 3% tổng lượngdòng chảy năm Sang tháng V, VI một số con lũ tiểu mãn xuất hiện làm lượng dòngchảy mùa cạn đã tăng lên một ít, nhưng còn ở mức thấp Tháng VII, VIII dòng chảytrên các khu vực tiếp tục giảm chậm, nhiều năm dòng chảy nhỏ nhất năm xảy ratrong thời kỳ này Tính chung trong toàn tỉnh lượng dòng chảy 8 tháng mùa cạn chỉchiếm khoảng 34% – 35% tổng lượng dòng chảy năm, trong đó nhu cầu dùng nướctrong thời kỳ này lại rất lớn, nếu như không có biện pháp tích trữ và sử dụng nướchợp lý sẽ bị thiếu nước nghiêm trọng

Dòng chảy mùa lũ trên các sông tỉnh Khánh Hòa có xu hướng giảm trongnhững năm gần đây, số trận lũ ngày càng ít, những năm xuất hiện ít lũ và mực nước

đỉnh lũ thấp ngày càng nhiều Xen kẽ những năm có đỉnh lũ thấp lại có những năm

có đỉnh lũ cao, giữa năm 2002 và 2004 là hai năm có mực nước đỉnh lũ năm rất thấpthì xen vào đó là năm 2003 có mực nước đỉnh lũ rất cao Năm 2006 và 2012 là năm

có đỉnh lũ năm rất thấp thì năm 2009 lại xuất hiện lũ lịch sử Như vậy dòng chảy lũcàng về sau càng không ổn định, có những năm mực nước đỉnh lũ năm rất thấp, xen

kẽ vào đó có những năm xuất hiện mực nước đỉnh lũ năm rất cao, tần suất xuất hiện

sự bất ổn định dòng chảy lũ ngày càng nhiều Đây có thể là một trong những biểuhiện của tác động biến đổi khí hậu (BĐKH) đến dòng chảy lũ, những trận mưa lớntập trung trong thời gian ngắn xuất hiện càng nhiều nên trận lũ lớn xuất hiện nhiềuhơn, số trận lũ ngày càng ít đi, mùa lũ có xu hướng ngắn đi, trong đó đặc trưng củaBĐKH là tăng tính cực đoan của dòng chảy trong đó có dòng chảy lũ thể hiện ngàycàng rõ[7]

Hơn thế nữa, trong những năm gần đây, tác động của thiên tai lũ lụt vàBĐKH đã gây ra những hậu quả nghiêm trọng, ảnh hưởng đến đời sống kinh tế và

xã hội của người dân trên lưu vực, đặc biệt là khu vực Thành phố Nha Trang phíadưới hạ lưu Để chỉ ra được sự tác động đó như thế nào thì cần phải hiểu được sựbiến đổi khí hậu tác động đến dòng chảy và các hiện tượng thủy văn theo mùa vàcon lũ như thế nào Kiến thức và các thông tin này là cần thiết cho việc quy hoạch

sử dụng đất để bảo vệ nguồn nước và quản lý hiệu quả lưu vực và cũng rất quantrọng đối với môi trường sinh thái và hoạt động kinh tế xã hội của lưu vực

Để khảo sát bài toán biến đổi dòng chảy trong tương lai ở lưu vực sông CáiNha Trang trong bối cảnh của biến đổi khí hậu, việc sử dụng các mô hình diễn toánmưa – dòng chảy là phù hợp do các kịch bản BĐKH mới chỉ cung cấp các thông tin

về sự biến đổi của các yếu tố khí tượng như nhiệt độ, lượng mưa, Một số mô hìnhmưa - dòng chảy được ứng dụng nhiều như: SWAT, HEC-HMS, MIKE-SHE, SAC-SMA, NASIM, HBV v.v… trong đó, mô hình SWAT (Arnold và cộng sự, 2002) đãđược chứng minh là một công cụ hiệu quả để đánh giá tài nguyên nước và

ô nhiễm với phạm vi lớn và các điều kiện môi trường trên toàn cầu Đồng thời, môhình SWAT còn được xây dựng để đánh giá tác động của việc sử dụng đất, của xóimòn và việc sử dụng hoá chất trong nông nghiệp trên một hệ thống lưu vực sông

10

Bên cạnh đó, mô hình SWAT cũng đã được nhiều tác giả ứng dụng để khảo sát sựbiến đổi của dòng chảy dưới các kịch bản về khí hậu và mặt đệm khác nhau và đãchứng tỏ được những ưu điểm Chính vì vậy luận văn đã lựa chọn sử dụng mô hình

SWAT để thực hiện mục tiêu của đề tài, và do đó, “Đánh giá tác động của biến đổi

khí hậu đến dòng chảy trạm Đồng Trăng, sông Cái Nha Trang, Tỉnh Khánh Hoà” là một đề tài có tính khoa học và thực tiễn, sẽ góp phần lớn đáp ứng yêu cầu

phòng chống lũ lụt cho người dân và chính quyền địa phương nhằm giảm thiểu đếnmức thấp nhất thiệt hại do lũ lụt gây ra Việc đánh giá tác động của biến đổi khí hậuđến dòng chảy trạm Đồng Trăng sông Cái Nha Trang tỉnh Khánh Hòa còn nhằm tạotiền đề cho việc xây dựng các giải pháp giảm thiểu những tác hại do biến đổi khíhậu, giúp nhà quản lý, các nhà hoạch định chính sách xác định chiến lược phát triểnkinh tế bền vững và đảm bảo an sinh xã hội Luận văn có bố cục gồm 3 chươngcùng với phần mở đầu, kết luận, phụ lục và tài liệu tham khảo:

Chương 1: Tổng quan

Chương 2: Mô hình SWAT

Chương 3: Ứng dụng mô hình SWAT để đánh giá tác động của biến đổi khíhậu đến dòng chảy lũ trạm Đồng Trăng, sông Cái Nha Trang, tỉnh Khánh Hòa

Kết luận và kiến nghị

CHƯƠNG 1 - TỔNG QUAN 1.1 Đặc điểm địa lý tự nhiên

1.1.1 Vị trí địa lý

Sông Cái Nha Trang là con sông lớn nhất tỉnh Khánh Hòa có vị trí địa lý nằmtrong khoảng 120°03' - 120°37' vĩ độ Bắc, 108°41' - 109°12' kinh độ Đông Sôngbắt nguồn từ đỉnh núi phía Tây Bắc xã Khánh Thượng cao 1477.5m, chảy theohướng Tây Bắc - Đông Nam Khi đến thôn Trang xã Khánh Thượng sông đổi sanghướng Tây - Đông và nhập với sông Chò tại thôn 1 xã Diên Đồng Sông tiếp tụcchảy theo hướng Tây - Đông và nhập với sông Suối Dầu tại Cầu Hà Dừa - thị trấnDiên Khánh, sông đổ ra vịnh Nha Trang tại cầu Trần Phú Diện tích lưu vực2000km2, chiều dài sông chính 79km, độ rộng bình quân lưu vực là 25.3km với hệ

số uốn khúc 1.4, hệ số hình dạng là 0.3, độ dốc sông 3.70/00, mật độ lưới sông0.8 km/km2

Lưu vực sông Cái Nha Trang bao trùm toàn bộ Thành phố Nha Trang, huyệnDiên Khánh, Khánh Vĩnh và một phần diện tích của huyện Cam Lâm tỉnh KhánhHòa và MaĐrăk tỉnh ĐăkLăk Phía bắc giáp lưu vực sông Dinh Ninh Hòa, phíaNam giáp lưu vực sông Cái Phan Rang, phía Tây giáp lưu vực sông ĐăkRông, phíaĐông giáp Biển Đông

1.1.2 Địa hình, địa mạo

Địa hình lưu vực sông Cái Nha Trang có xu hướng dốc dần từ Tây sangĐông, đỉnh núi cao nhất ở phía tây là đỉnh Hòn Giao Bắc có độ cao 2038.2m, khuvực trung du có độ cao phổ biến từ 10 đến 25m, khu vực đồng bằng có độ cao phổbiến từ 2 đến 10m Độ dốc trung bình vùng núi là 16%, độ dốc trung bình lưu vực2,8% Lưu vực sông Cái Nha Trang nằm ở sườn Đông của dãy Trường Sơn, địa hình

bị chia cắt nhiều bởi những ngọn núi, khu vực ven biển có những dãy núi đâmngang ra biển Khu vực đồng bằng và trung du xen kẽ những đỉnh núi nhỏ có độ cao

từ 40 đến 50m, cá biệt có đỉnh núi Chín Khúc có độ cao 592.6m, các quả đồi nhỏ có

độ cao 20 đến 30m Địa hình cao nhất và dốc nhất là khu vực phía Tây nam của lưu

12

vực, khu vực này là thượng nguồn của các nhánh sông đổ vào sông Thác Ngựa lànhánh sông chính của sông Cái Nha Trang.

Hình 1.1 Sơ đồ lưu vực sông Cái Nha Trang 1.1.3 Địa chất, thổ nhưỡng

- Đất mùn vàng trên núi cao 900-1000m

- Đất thung lũng có thành phần cơ giới nhẹ đến trung bình và đất phù sa phân bố dọc các sông suối trong tỉnh

- Đất pha cát thành phần cơ giới nhẹ và thô, kết cấu rời rạc, phân bố phần lớn vùng ven biển phía đông

Bảng 1.1 Diện tích các loại đất đai tỉnh Khánh Hòa

VII.Nhóm đất mùn vàng đỏ trên núi

14

1 Đất mùn vàng đỏ trên đá Macma acid

sự phân hóa về thảm phủ thực vật như sau:

- Khu vực đỉnh núi ở độ cao trên dưới 2000m: Thảm thực vật thân gỗ chỉ caokhoảng 7 - 10m, tán xen kẽ Sự phân tầng của thảm thực vật không thật sự rõ rệt Ởđây chủ yếu là sự góp mặt của các loài thực vật có nhiều thân trên một gốc Trênthân cây phủ lớp rêu mỏng Một số loài thực vật bì sinh phát triển trên cành, nhánhcây Điểm đáng lưu ý là các sườn của đỉnh núi rất dốc, nhưng vẫn được bao phủ bởithảm thực vật rừng khá dày với sức sống tốt Dưới tán rừng, lớp lá rụng 3 - 4cm vớiquá trình phân giải chậm Cây có độ cao đạt tới 10m, đường kính 50 - 100cm, tán đan xen Cây có sự phân cành sớm, ở độ cao khoảng 2 - 3m

- Khu vực đỉnh - sườn núi ở độ cao 1400 - 1700m: Rừng á nhiệt đới thườngxanh cây lá rộng và cây lá rộng, lá kim núi trung bình Rừng có diện tích lớn với cây

đa trội và cấu trúc thảm phức tạp Rừng phân ra nhiều tầng nhưng có 2 tầng chủđạo, các tầng trung gian không liên tục, xen kẽ nhau Ngoài ra ở độ cao này còn córừng lá kim núi trung bình (rừng thông 3 lá tự nhiên)

- Khu vực có độ cao dưới 1400m: Rừng chủ yếu là rừng trồng, loại cây lárộng thường là cây keo, cao su, cà phê, hồ tiêu, ca cao và cây lá kim là cây thông

- Khu vực trung du có độ cao dưới 100m thường là các trảng cây cỏ bụi, cáccây lá rộng với tán cây có đường kính từ 2 đến 4m, thân cây có đường kính từ 10đến 30cm, độ cao từ 2 đến 5m Các cây được trồng chủ yếu là mía và cây keo từ 2 đến 5 năm tuổi

Theo niên giám thống kê tỉnh Khánh Hòa, khu vực trung du và đồng bằnglưu vực sông Cái Nha Trang là đất nông nghiệp với diện tích 1262.2km2, đất phinông nghiệp với diện tích 131.1km2, đất chưa sử dụng 363.9km2.

a Mưa

Do địa hình phức tạp nên lượng mưa giữa các khu vực có sự chênh lệch nhaukhá lớn Lượng mưa trung bình năm trên lưu vực sông Cái Nha Trang thì lớn hơnkhoảng từ 1300 - 1931mm Tổng lượng mưa mùa mưa khoảng 900 - 1059mm,chiếm khoảng 67 - 75% lượng mưa năm; tổng lượng mưa mùa khô khoảng 300 -450mm Những trận mưa, lũ lớn chủ yếu tập trung vào tháng X và XI

b. Nhiệt độ không khí

Nhiệt độ trung bình năm ở vùng đồng bằng và ven biển dao động từ 26,3°Cđến 26,9°C, lên đến độ cao 400m nhiệt độ giảm xuống khoảng 23,0 -24,0°C Nhìnchung trong năm ít xảy ra các đợt nắng nóng kéo dài nhiều ngày gây ảnh hưởng lớnđến sản xuất và đời sống

c Độ ẩm tương đối

Độ ẩm khá thấp, trung bình năm vào khoảng 80% tại Nha Trang Hàng nămchỉ có 2 - 3 tháng đầu mùa đông (tháng X, XI, XII) là khá ẩm với độ ẩm trung bình85% Còn trong nửa cuối mùa đông độ ẩm giảm xuống 80 - 83% Tháng ẩm nhất làtháng XI có độ ẩm khoảng 85 - 87%

d Bốc hơi

16

Tổng lượng bốc hơi trung bình nhiều năm dao động từ 1000 - 1100mm/nămtức là bằng 2/3 lượng mưa Trong 3 tháng (từ tháng VI - VIII), mỗi tháng lượng bốchơi đạt tới 120 - 150mm, vượt quá lượng mưa các tháng này Thời kỳ bốc hơi ít nhất

là các tháng mùa mưa, từ tháng X - XII, lượng bốc hơi chỉ khoảng 60 - 80mm

bờ biển

Khi xét ảnh hưởng của lũ và triều ta thấy khi gặp triều lên mực nước có thểdâng cao hơn trường hợp không có triều từ 20 - 30cm, khi triều xuống mực nước cóthể chênh lệch so với mực nước không ảnh hưởng triều từ 9 - 10cm

Dòng chảy năm

Dòng chảy trên các sông chủ yếu do mưa cung cấp, nên sự phân bố của dòngchảy tương tự sự phân bố của mưa Độ sâu dòng chảy trên lưu vực sông Cái NhaTrang là 1159mm

Bảng 1.2 Cán cân nước các lưu vực[7]

STT

1 2 3

Lượng dòng chảy 4 tháng mùa lũ chiếm từ 65% - 66% lượng dòng chảy cảnăm, lượng dòng chảy 8 tháng mùa cạn chỉ chiếm từ 34% - 35% lượng dòng chảy

cả năm Tổng lượng dòng chảy (W) sông Cái Nha Trang 2319km3

Dòng chảy lũ giai đoạn những năm gần đây

a) Chế độ dòng chảy mùa lũ

Mùa lũ trên các sông suối tỉnh Khánh Hòa bắt đầu từ tháng IX đến tháng XII.Nguyên nhân chủ yếu do vùng thấp, áp thấp, bão hoặc cũng có thể do không khílạnh kết hợp với vùng thấp Hàng năm trên các sông suối tỉnh Khánh Hòa xuất hiện

từ 1 đến 4 trận lũ, trung bình có 2,3 trận lũ trong năm, cá biệt một số năm khôngxuất hiện lũ Số lượng các trận lũ ở tỉnh Khánh Hòa không nhiều so với các tỉnhkhác ở khu vực Nam Trung Bộ do đặc điểm khí hậu ở đây ôn hòa hơn, không khắcnghiệt như các vùng khác thuộc khu vực Nam Trung Bộ Các trận lũ lớn chủ yếuxảy ra từ tháng X đến tháng XI, đó là thời kỳ hoạt động của xoáy thuận nhiệt đới,bão ảnh hưởng mạnh nhất đối với khu vực Nam Trung Bộ, mặt khác đây cũng làthời kỳ rãnh thấp xích đạo đi qua khu vực tỉnh Khánh Hòa

Tốc độ dòng chảy lớn nhất trên sông Cái Nha Trang đạt 2.96m/s xuất hiện vàonăm 1999 Nhìn chung, tốc độ dòng chảy của các sông tỉnh Khánh Hòa không lớn lắm.Tốc độ dòng chảy chịu ảnh hưởng bởi khả năng tập trung nước sườn dốc lưu vực, độdốc và mặt cắt ngang sông Trong những năm gần đây diện tích rừng bị thu hẹp nhiều

và mức độ che phủ rừng giảm do chặt phá rừng đầu nguồn làm mức độ tập trung dòngchảy sườn dốc nhanh Độ dốc địa hình và lòng sông giảm trong những năm gần đây do

18

bào mòn bề mặt lưu vực mạnh, lòng sông độ dốc giảm do xói mòn đáy sông vùngthượng lưu vì vùng núi là khu vực có dòng chảy rối và dòng chảy xiết Hiện tượng khaithác cát cũng làm tăng đáng kể diện tích mặt cắt ngang sông Do giảm độ dốc bề mặt,lòng sông và mở rộng mặt cắt ngang đã là nguyên nhân tốc độ dòng chảy trong sôngnhững năm gần đây giảm so với những năm trước Ngoài ra, tốc độ dòng chảy trongsông còn phụ thuộc chế độ mưa, những biến đổi lượng mưa là nguyên nhân chính thayđổi chế độ dòng chảy sông ngòi trong đó có vận tốc dòng chảy trong sông.

Biên độ lũ trung bình tại trạm Đồng Trăng là 6.2m, lớn nhất là 9.11m trongtrận lũ tháng XI năm 2009 Cường suất lũ trung bình tại trạm Đồng Trăng là0.75m/h, lớn nhất là 2.01m/h xuất hiện trong trận lũ tháng XI năm 1978 Biên độ vàcường suất lũ tại trạm thủy văn phụ thuộc vào đặc điểm mưa của từng đợt mưa trêncác lưu vực sông Đối với tỉnh Khánh Hòa, các lưu vực sông thuộc loại vừa và nhỏ,trong khi đó các trận mưa lũ lớn xảy ra trên diện rộng do các hình thế thời tiết lớngây ra nên đặc điểm mưa của các trận mưa lớn trên các lưu vực sông có sự khácnhau là do tác động của địa hình Ngoài ra, địa hình, thảm phủ thực vật, đặc điểmhình dạng sông cũng tác động đáng kể đến biên độ và cường suất lũ Nhìn chung,cường suất và biên độ lũ vùng thượng lưu lớn hơn vùng hạ lưu, lưu vực sông hìnhdạng cành cây như sông Cái Nha Trang có biên độ và cường suất lũ cao

là trận lũ có cường suất lũ lên nhanh

19

1350 1150 950 750

550 350

Hình 1.2 Đường mực nước trận lũ lịch sử năm 2009

Trận lũ lớn năm 1998 với mực nước đỉnh lũ đạt 13.12m vào hồi 18h ngày 11

tháng XII, cường suất lũ lên cao nhất đạt 0.50m/h, biên độ lũ đạt 7.37m Lưu lượng

đỉnh lũ đạt 2,037m3/s, thời gian xuất hiện từ ngày 9 đến ngày 12 tháng XII Đường

quá trình lũ lên từ từ nhưng xuống nhanh

1350

cm1250

1150 1050 950 850 750 650 550

Hình 1.3 Đường mực nước trận lũ lớn năm 1998

Trận lũ lớn năm 2010 với mực nước đỉnh lũ đạt 12.66m xuất hiện vào lúc

15h ngày 05 tháng XI, cường suất lũ lớn nhất đạt 0.79m/h, biên độ lũ đạt 8,14m

Lưu lượng đỉnh lũ đạt 1,890m3/s, thời gian của trận lũ kéo dài từ ngày 29 tháng X

đến ngày 13 tháng XI Trận lũ này có đường quá trình diễn biến rất phức tạp với lũ

20

Qua số liệu thống kê ở trên có thể thấy rằng các trận lũ lớn thừờng xảy ra vào

tháng XI kéo dài 4-5 ngày Càng về sau, lũ xuất hiện càng phức tạp hơn gây bất lợi

nghiêm trọng cho người dân phía dưới hạ lưu

1350cm1150

1.1.7 Mạng lưới trạm khí tượng thủy văn trên lưu vực sông

Hiện trên địa bàn tỉnh Khánh Hoà, có tất cả 4 trạm đo cơ bản (Hình 1.1)

Trong đó, có 1 trạm khí tượng Nha Trang, 1 trạm mưa Khánh Vĩnh, 2 trạm thuỷ văn

là Diên An đo mực nước và Đồng Trăng đo lưu lượng, mực nước, mưa

Bảng 1.3 Các trạm đo KTTV trên lưu vực sông

21

38279 tấn Trên lưu vực có 15 doanh nghiệp nhà nước (chiếm 78,9% so với toàntỉnh), 488 doanh nghiệp tư nhân (chiếm 75,9%), 11 doanh nghiệp có vốn đầu tưnước ngoài (chiếm 39,3%) Có 2551 cơ sở sản xuất, chiếm 39,3% và 7586 côngnhân, chiếm 37,9% Doanh nghiệp thương mại là 20 (chiếm 100%), doanh nghiệpvận tải là 3 (chiếm 75%) Trọng điểm của sản xuất kinh doanh thuộc địa phận thànhphố Nha Trang, trong đó có ngành du lịch rất phát triển.

* Một số khu công nghiệp và vùng kinh tế trọng điểm

Trên lưu vực sông Cái Nha Trang bao gồm huyện Khánh Vĩnh, huyện DiênKhánh và Thành phố Nha Trang Trong đó về kinh tế TP Nha Trang là trung tâmchính trị, văn hóa, kinh tế của tỉnh Khánh Hòa Sự phát triển của Nha Trang luônđược đánh giá là có tính chất đầu tàu, làm động lực phát triển chung cho cả tỉnh.Hiện nay, cơ cấu kinh tế của TP Nha Trang đã hình thành khá rõ nét theo hướngdịch vụ, du lịch - công nghiệp - nông nghiệp Hướng tới, TP Nha Trang sẽ tiếp tụcđầu tư phát triển mạnh lĩnh vực dịch vụ, du lịch nhằm xây dựng TP Nha Trangthành một trung tâm du lịch có tầm cỡ quốc gia, quốc tế

- Khu công nghiệp Suối Dầu: đã có 21 doanh nghiệp trong và ngoài nướcđang hoạt động; khu công nghiệp Ninh Thủy đang xây dựng; ngoài ra các khu côngnghiệp Nam và Bắc như: Cam Ranh, Vạn Ninh chuẩn bị xây dựng sẽ tạo điều kiệnthuận lợi cho các nhà đầu tư

- Cụm công nghiệp Bình Tân (Nha Trang): bao gồm các ngành sản xuấtchính là chế biến thủy sản, dệt, nhuộm, thuốc lá, song mây Trong cụm có nhà máyDệt Tân Tiến, nhà máy thuốc lá Khánh Hòa, xí nghiệp chế biến song mây xuất khẩuNha Trang, công ty sản xuất hàng thủ công xuất khẩu Nha Trang

- Cụm công nghiệp Suối Hiệp (Diên Khánh): bao gồm các ngành sản xuấtchính là đường Diên Khánh, nhà máy bia Rồng Vàng, phân xưởng sản xuất cồncông nghiệp, phân xưởng sản xuất giấy Hoa Hồng, nước ngọt Suối Tiên và cụm khucông nghiệp Suối Dầu

- Cụm công nghiệp Hòn Khô (Nha Trang): bao gồm các ngành khai thác vàchế biến thủy sản, sản xuất vật liệu xây dựng

1.2.2 Đặc điểm xã hội

* Dân số

Theo số liệu niên giám thống kê tỉnh Khánh Hoà năm 2009 dân số trên lưuvực sông Cái Nha Trang là 552060 người (trong đó nam là 271602 người chiếm49,2%, nữ là 280458 người chiếm 50,8%) Dân số ở thành thị là 317496 ngườichiếm 57,5%, nông thôn là 234564 người chiếm 42,5%; nơi có mật độ dân số caonhất là Thành phố Nha Trang với 1540 người/km2

* Giáo dục

Trên lưu vực sông Cái Nha Trang có 82 trường tiểu học, chiếm 25,6% so vớitoàn tỉnh; 37 trường trung học cơ sở, chiếm 37%; 14 trường trung học phổ thông,chiếm 50% Giáo viên tiểu học là 2010, chiếm 42,4% so với toàn tỉnh; 1756 giáoviên trung học cơ sở, chiếm 41,8%; 873 giáo viên trung học phổ thông, chiếm47,9% Trên lưu vực có 5 trường đại học, 3 trường cao đẳng, tập trung ở thành phốNha Trang

* Y tế

Tổng số giường bệnh các huyện, thành phố trên lưu vực sông Cái Nha Trang

là 1941, chiếm 64,5% so với toàn tỉnh, có 410 bác sĩ (chiếm 68,8%), 399 y sĩ(chiếm 60,1%), 600 y tá (chiếm 64,2%), 20 dược sĩ cao cấp (chiếm 76,9%), 124dược sĩ trung cấp (chiếm 76%), 39 dược tá (chiếm 52%) [3]

1.3 Tổng quan về biến đổi khí hậu

1.3.1 Khái niệm

Biến đổi khí hậu, theo cách sử dụng của IPCC, chỉ sự biến đổi trong trạng thái khí hậu nhận biết được thông qua những thay đổi về giá trị trung bình hoặc tính chất của nó diễn ra trong một thời đoạn dài hàng thập kỷ hoặc hơn thế Nó chỉ ra bất

cứ thay đổi nào của khí hậu theo thời gian cho dù là do biến đổi tự nhiên hay do tác động của con người [1]

1.3.2 Biểu hiện về BĐKH

a) Quy mô toàn cầu

Nhiệt độ

23

Nhiệt độ trung bình toàn cầu tăng khoảng 0,89oC (dao động từ 0,69 đến1,08oC) trong thời kỳ 1901-2012 Nhiệt độ trung bình toàn cầu có chiều hướng tăngnhanh đáng kể từ giữa thế kỷ 20 với mức tăng khoảng 0,12oC/thập kỷ trong thời kỳ1951-2012.

Giáng thủy trung bình toàn cầu kể từ năm 1901 có xu thế tăng ở vùng lục địa

vĩ độ trung bình thuộc bắc bán cầu

Số ngày và số đêm lạnh có xu thế giảm, số ngày và số đêm nóng cùng vớihiện tượng nắng nóng có xu thế tăng rõ rệt trên quy mô toàn cầu từ khoảng năm

1950 Mưa lớn có xu thế tăng trên nhiều khu vực, nhưng lại giảm ở một số ít khuvực

Hình 1.5.Chuẩn sai nhiệt độ trung

bình toàn cầu thời kỳ 1850-2012

(so với thời kỳ 1961-1990)

Hình 1.7 Biến đổi của nhiệt độ trung bình năm thời kỳ 1901-2012

(Nguồn: IPCC, 2013)

Chú thích: Các ô lưới được thể hiện (được tô màu) nếu đảm bảo điều kiện: có đủ tối thiểu 70% số liệu trong thời kỳ 1901-2012; trong đó, tối thiểu giai đoạn đầu chỉ được thiếu 20% số liệu và giai đoạn cuối thiếu 10% số liệu Những ô lưới màu trắng (không được tô màu) là những ô không đảm bảo điều kiện tính toán Những ô được đánh dấu + là những ô lưới có xu thế biến đổi ở mức ý nghĩa 10% (hay mức tin cậy 90%) trở lên.

25

Hình 1.8 Biến đổi của lƣợng mƣa năm thời kỳ 1901-2010 và thời kỳ 1951-2010

(Nguồn: IPCC, 2013)

b) Biến đổi khí hậu ở Việt Nam

Việt Nam được đánh giá là một trong những quốc gia bị ảnh hưởng nặng nềnhất của biến đổi khí hậu Trong những năm qua, dưới tác động của biến đổi khíhậu, tần suất và cường độ các thiên tai ngày càng gia tăng, gây nhiều tổn thất to lớn

về người, tài sản, cơ sở hạ tầng, về kinh tế, văn hoá, xã hội, tác động xấu đến môitrường Tác động của biến đổi khí hậu đối với nước ta là rất nghiêm trọng, là nguy

cơ hiện hữu cho mục tiêu xóa đói giảm nghèo, cho việc thực hiện các mục tiêu thiênniên kỷ và sự phát triển bền vững của đất nước Việt Nam đã rất nỗ lực ứng phó vớibiến đổi khí hậu, thể hiện qua các chính sách và các chương trình quốc gia

Năm 2009, Bộ Tài nguyên và Môi trường đã xây dựng và công bố kịch bảnbiến đổi khí hậu và nước biển dâng cho Việt Nam và được cập nhật qua các năm

2011, năm 2016

Với phiên bản cập nhật năm 2016, mô hình khí hậu toàn cầu và khu vực lànhững công cụ chính được sử dụng để đánh giá xu thế và mức độ biến đổi của khíhậu tương lai, đặc biệt là các cực đoan khí hậu Các mô hình được sử dụng trongtính toán xây dựng kịch bản biến đổi khí hậu cho Việt Nam, gồm: (i) Mô hìnhAGCM/MRI của Viện Nghiên cứu Khí tượng Nhật Bản, (ii) Mô hình PRECIS củaTrung tâm Hadley - Vương quốc Anh, (iii) Mô hình CCAM của Cơ quan Nghiên

cứu Khoa học và Công nghiệp Liên bang Úc (CSIRO), (iv) Mô hình RegCM củaTrung tâm quốc tế về Vật lý lý thuyết của Ý (ICTP), (v) Mô hình clWRF của Mỹ.

Nhiệt độ

Nhiệt độ có xu thế tăng ở hầu hết các trạm quan trắc, tăng nhanh trong nhữngthập kỷ gần đây Trung bình cả nước, nhiệt độ trung bình năm thời kỳ 1958-2014tăng khoảng 0,62oC, riêng giai đoạn (1985-2014) nhiệt độ tăng khoảng 0,42oC Tốc

độ tăng trung bình mỗi thập kỷ khoảng 0,10oC, thấp hơn giá trị trung bình toàn cầu(0,12oC/thập kỷ, IPCC 2013)

Nhiệt độ tại các trạm ven biển và hải đảo có xu thế tăng ít hơn so với cáctrạm ở sâu trong đất liền Có sự khác nhau về mức tăng nhiệt độ giữa các vùng vàcác mùa trong năm Nhiệt độ tăng cao nhất vào mùa đông, thấp nhất vào mùa xuân.Trong 7 vùng khí hậu, khu vực Tây Nguyên có mức tăng nhiệt độ lớn nhất, khu vựcNam Trung Bộ có mức tăng thấp nhất[1]

Hình 1.9 Chuẩn sai nhiệt độ ( o C) trung bình năm (a) và nhiều năm (b) trên

quy mô cả nước

27

Hình 1.10 Chuẩn sai nhiệt độ trung bình năm ( o C) đối với các trạm ven biển

Đối với các khu vực phía Bắc, lượng mưa chủ yếu giảm rõ nhất vào cáctháng mùa thu và tăng nhẹ vào các tháng mùa xuân Đối với các khu vực phía Nam,

lượng mưa các mùa ở các vùng khí hậu đều có xu thế tăng; tăng nhiều nhất vào cáctháng mùa đông (từ 35,3% ÷ 80,5%/57 năm) và mùa xuân (từ 9,2% ÷ 37,6%/57năm)[1].

Bảng 1.4 Thay đổi lƣợng mƣa (%) trong 57 năm qua (1958-2014)

1.3.3 Kịch bản biến đổi khí hậu cho khu vực nghiên cứu

Trong khuôn khổ của luận văn, tác giả sử dụng kịch bản RCP4.5 và RCP 8.5

về sự thay đổi % lượng mưa do biến đổi khí hậu cho lượng mưa 5 ngày trung bình

lớn nhất (Bảng 1.5) và mưa trung bình theo mùa (Bảng 1.6) cho lưu vực nghiên cứu

được chiết xuất từ báo cáo “Kịch bản biến đổi khí hậu, nước biển dâng cho ViệtNam”(năm 2016) [1] để tính toán lưu lượng dòng chảy 5 ngày lớn nhất, dòng chảytrung bình tháng mùa lũ vào giữa thế kỷ, cuối thế kỷ Từ đó đánh giá tác động củabiến đổi khí hậu đến dòng chảy lũ trạm Đồng Trăng, sông Cái Nha Trang tỉnhKhánh Hòa

Bảng 1.5 % Thay đổi lƣợng mƣa 5 ngày lớn nhất theo các kịch bản

Yếu tố

5 ngàylớnnhất

29

Bảng 1.6 % thay đổi lƣợng mƣa mùa theo các kịch bản

Bùi Văn Chanh (2013)[2] đã ứng dụng mô hình thủy lực một và hai chiều kếthợp HDM (Hydro Dynamic Model) xây dựng bản đồ ngập lụt hạ lưu sông Cái Nha

3%, 5% và 10% cho vùng hạ lưu sông Cái Nha Trang Mô phỏng diễn biến ngập lụttheo thời gian, tạo bản đồ ngập lụt và tích hợp chương trình sử dụng Kết quả tínhtoán dự tính vùng ngập lụt hạ lưu sông Cái Nha Trang bao gồm hầu hết vùng đồngbằng của Huyện Diên Khánh và ngoại thành Nha Trang Ngoài ra nhánh sông Chòcòn gây ngập lụt lớn tại thôn Khánh Xuân xã Diên Lâm Vùng nội thành Nha Tranghầu như không bị ảnh hưởng ngập lụt của lũ sông Cái mà chỉ ảnh hưởng ngập úng

do mưa lớn gây ra Vùng ngập sâu nhất là xã Diên Thạnh, Diên Toàn và Diên Lạccủa Huyện Diên Khánh Diện ngập và độ sâu thay đổi theo các kịch bản với tần suất

lũ 1%, 3%, 5% và 10%

Nguyễn Ý Như (2012)[8] đã đánh giá dòng chảy năm tỉnh Khánh Hòa trongbối cảnh biến đổi khí hậu Mô hình MIKE – NAM được áp dụng đối với tỉnh KhánhHòa để đánh giá tác động của biến đổi khí hậu đối với dòng chảy Biến đổi khí hậuđược dự tính theo Bộ tài nguyên và Môi trường (kịch bản B1, B2 và A2) cho giai

30

đoạn 2050 – 2100 Theo kịch bản lượng mưa trung bình năm tăng 0.7 đến 1.7% vànhiệt độ tăng 0.4 – 1oC Kịch bản biến đổi khí hậu được sử dụng như tác nhân bênngoài đối với mô hình mưa dòng chảy NAM, được cài đặt và thực hiện cho 18 tiểulưu vực Dòng chảy trung bình năm trên lưu vực dao động trong khoảng 1- 7m3/s(31 – 45%) khi so sánh kịch bản với giai đoạn nền Những biến đổi lớn hơn đượcnhận thấy ở giá trị cực trị Dòng chảy theo mùa được dự đoán thay đổi với dòngchảy cao hơn đáng kể trong mùa mưa, giảm đáng kể trong dòng chảy mùa kiệt.

Đặng Đình Đức (2012)[4] đã khôi phục số liệu dòng chảy tỉnh Khánh Hòabằng mô hình NAM Lí do bởi các lưu vực sông thuộc tỉnh Khánh Hòa có nguồn sốliệu dòng chảy đo đạc rất hạn chế Toàn tỉnh chỉ có duy nhất trạm thủy văn ĐồngTrăng có số liệu dòng chảy đo đạc (từ năm 1983 tới nay), trong khi đó số trạm đomưa trong tỉnh tương đối đầy đủ và đồng bộ từ năm 1977 tới nay Bài toán đặt ra làcần phải khôi phục số liệu dòng chảy từ tài liệu đo mưa Tác giả đã sử dụng phươngpháp mô hình toán, cụ thể là ứng dụng mô hình NAM để khôi phục số liệu dòngchảy các lưu vực sông thuộc tỉnh Khánh Hòa Bài báo đã tính toán lưu lượng trungbình tháng từ 1977-2010 cho 18 tiểu lưu vực thuộc tỉnh Khánh Hòa Qua số liệutính toán được từ mô hình NAM đã đưa ra được một bức tranh tương đối đầy đủ vềtài nguyên nước mặt của tỉnh, sự phân bố theo không gian và thời gian trong năm.Các số liệu dòng chảy tính toán cùng với số liệu dòng chảy thực đo sẽ là một nguồntài liệu quý giá, làm đầu vào cho các bài toán khác, đặc biệt là bài toán cân bằngnước, quy hoạch, sử dụng hợp lý tài nguyên nước của tỉnh

Lại Thị Hương (2012)[13] đã đánh giá tác động của biến đổi khí hậu đến tàinguyên nước mặt tỉnh Khánh Hoà của thời kỳ 2000-2009 so giai đoạn 1980-1999với việc sử dụng phương pháp thống kê Kết quả của nghiên cứu đưa ra kết luận:Biến đổi khí hậu tác động đến tài nguyên nước thông qua việc làm thay đổi lượngmưa và phân bố mưa giữa các vùng trên địa bàn Nhiệt độ tăng sẽ làm lượng nướcbốc hơi nhiều hơn kết hợp với địa hình đặc thù tỉnh Khánh Hoà dẫn đến lượng mưanhiều hơn Một hậu quả của BĐKH là thay đổi về thời gian mùa mưa Trong khoảng

5 năm trở lại đây, mùa mưa lũ trên các sông trong tỉnh đến sớm và cũng kết

thúc sớm hơn so với TBNN khoảng gần 01 tháng, thường bắt đầu từ tháng IX và kếtthúc vào giữa tháng XI.

Như vậy, các nghiên cứu thực hiện trên lưu vực sông Cái Nha Trang mặc dù

đã sử dụng mô hình thông số tập trung để đánh giá tác động biến đổi khí hậu đếndòng chảy mặt trên địa bàn tỉnh Khánh Hoà, tuy nhiên vấn đề cấp thiết đặt ra ở đây

là cần cập nhật kịch bản biến đổi khí hậu bản mới nhất hiện nay cũng như sử dụng

mô hình thông số phân bố trong việc đánh giá tác động của biến đổi khí hậu đếndòng chảy mới mang ý nghĩa cho thực tiễn Trong những năm gần đây, Việt Nam đã

có một số công trình nghiên cứu đánh giá tác động của BĐKH đến lũ lụt cũng nhưchứng minh hiệu quả của SWAT trên nhiều lưu vực[6, 7, 8, 9, 10, 11, 20] Chính vìvậy, để giải quyết bài toán này mô hình SWAT được lựa chọn là công cụ chính Hơnthế nữa, mô hình đánh giá được thay đổi thuỷ văn trong lưu vực sông Cái NhaTrang dưới tác động của BĐKH và phản ánh đuợc những thay đổi các đặc trưngmùa lũ như: dòng chảy mùa lũ, dòng chảy 5 ngày lớn nhất Thời gian gần đây, nó đã

mô phỏng được khá tốt cho dòng chảy ngày và giờ [15, 18, 19, 21, 22] là điều kiệntốt cho việc tính toán dòng chảy mùa lũ và dòng chảy 5 ngày lớn nhất dưới tác độngcủa biến đổi khí hậu theo các kịch bản

32

CHƯƠNG 2 - GIỚI THIỆU MÔ HÌNH SWAT 2.1 Xuất xứ của mô hình SWAT

Mô hình “Công cụ đánh giá đất và nước“ SWAT (Soil and Water AssementTools) là một mô hình vật lý được xây dựng từ những năm 90 do tiến sỹ Jeff Arnoldthuộc trung tâm nghiên cứu đất nông nghiệp USDA- Agricultural Research Service(ARS) xây dựng nên Mô hình này được xây dựng để mô phỏng ảnh hưởng của việcquản lý sử dụng đất đến nguồn nước, bùn cát và hàm lượng chất hữu cơ trong đấttrên hệ thống lưu vực sông trong một khoảng thời gian nào đó Tiền thân của môhình SWAT là mô hình SWRRB ((Simulator for Water Resources in Rural Basins)(Williams et al., 1985; Arnold et al., 1990)) và mô hình ROTO ((Routing Outputs toOutlet) (Arnold et al., 1995))

Các phiên bản của mô hình SWAT là:

+ SWAT94.2: Mô phỏng đường lưu lượng đơn vị

+ SWAT96.2: Trong phiên bản này được cập nhật thêm phần quản lý về hàmlượng chất hữu cơ trong đất, trong đó nghiên cứu ảnh hưởng của sự thay đổi khí hậutới sự phát triển của cây trồng Phương trình chất lượng nước từ mô hình QUAL2Eđược sử dụng đến

+ SWAT98.1: Trong phiên bản này thêm phần diễn toán dòng chảy do tuyết tan, chất lượng nước trong sông

+ SWAT99.2: Phiên bản này cập nhật thêm diễn toán chất lượng nước cho hồchứa, phần thuỷ văn đô thị được cập nhật từ mô hình SWMM

+ SWAT2000: Cập nhật thêm phương trình thấm của Green & Ampt, cậpnhật thêm các yếu tố khí tượng thời tiết như bức xạ mặt trời, tố độ gió , cho phépgiá trị bốc thoát hơi tiềm năng của lưu vực được lấy từ giá trị đầu vào của mô hìnhhoặc được tính toán theo phương trình Đặc biệt trong phiên bản này có sử dụngmôi trường ARCVIEW là môi trường giao diện

+SWAT2012: Đặc biệt trong phiên bản này có sử dụng môi trường ARCGIS

là môi trường giao diện Trong luận văn này đã sử dụng phiên bản SWAT 2012 đểtính toán

2.2 Cấu trúc của mô hình

Mô hình SWAT mô phỏng quá trình thuỷ văn diễn ra trên lưu vực

Mô hình chia dòng chảy thành ba pha: pha mặt đất, pha dưới mặt đất (sát mặt, ngầm) và pha trong sông

 Pha mặt đất diễn tả các thành phần dòng chảy mặt, phần xói mòn

 Pha sát mặt diễn tả các thành phần dòng chảy sát mặt, dòng chảyngầm

 Pha trong sông diễn tả diễn toán lượng dòng chảy tới mặt cắt cửa ra của lưu vực

2.2.1 Pha đất của chu trình thuỷ văn

Chu trình thuỷ văn được mô tả trong mô hình SWAT dựa trên phương trìnhcân bằng nước như sau:

t

SW t = SW 0 +∑(R day − Q surf − E a − W seep − Q gw ) (2.1)

trong đó:

SWt là tổng lượng nước tại cuối thời đoạn tính toán (mm)

SWo là tổng lượng nước ban đầu tại ngày thứ i (mm)

T là thời gian (ngày)

Rday là tổng lượng mưa tại ngày thứ i (mm)

Qsurf là tổng lượng nước mặt của ngày thứ i (mm)

Ea là lượng bốc thoát hơi tại ngày thứ i (mm)

Wseep là lượng nước đi vào tầng ngầm tại ngày thứ i (mm)

Qgw là lượng nước hồi quy tại ngày thứ i (mm)

1 Các yếu tố khí hậu

Các yếu tố khí hậu của lưu vực cung cấp số liệu đầu vào của mô hình đểkiểm soát cân bằng nước và xác định mối liên quan giữa các thành phần khác nhautrong chu trình thuỷ văn Các biến khí hậu được sử dụng trong mô hình SWAT làmưa, nhiệt độ tối cao, nhiệt độ tối thấp, độ ẩm, vận tốc gió, bức xạ Các số liệu nàyđược lấy ở những trạm đo khí tượng và trạm đo mưa

34