Ảnh hưởng của nhiệt độ tăng do biến đổi khí hậu đến các thông số thủy văn và thành phần cân bằng nước tại lưu vực Dong Ér, Hungary

Bài viết Ảnh hưởng của nhiệt độ tăng do biến đổi khí hậu đến các thông số thủy văn và thành phần cân bằng nước tại lưu vực Dong Ér, Hungary xem xét ảnh hưởng của nhiệt độ tăng đến các thông số thủy văn và các thành phần cân bằng nước, dựa trên kết quả mô phỏng của các kịch bản. Kết quả mô phỏng năm 2018 sẽ là số liệu tham khảo để làm cơ sở so sánh độ nhạy của của việc biến đổi số liệu nhiệt độ đầu vào đến các thông số thủy văn và thành phần cân bằng nước trong bối cảnh biến đổi khí hậu.

ẢNH HƯỞNG CỦA NHIỆT ĐỘ TĂNG DO BIẾN ĐỔI KHÍ HẬU ĐẾN CÁC THÔNG SỐ THỦY VĂN VÀ THÀNH PHẦN CÂN BẰNG NƯỚC TẠI LƯU VỰC DONG-ÉR, HUNGARY Trần Quang Hợp1, Trường Đại học Tài nguyên Môi trường Hà Nội Đại học Szeged, Hungary Tóm tắt Các tượng thời tiết cực đoan biến đổi khí hậu xảy thường xuyên khu vực Trung Âu, có Hungary Trong tương lai gần, tượng hạn hán, mưa lớn gia tăng biến đổi khó lường trước, gây nên thiệt hại nghiêm trọng đến tài nguyên nước địa phương Kết nghiên cứu biến đổi khí hậu cho thấy, tượng thời tiết cực đoan khó dự báo trước, đặc biệt bối cảnh số liệu thực đo thiếu Để đánh giá ảnh hưởng biến đổi khí hậu đến môi trường hoạt động kinh tế - xã hội cần phải nghiên cứu, tìm hiểu mối quan hệ tổng hợp thông số, như: bốc hơi, dịng chảy, độ thấm, nước mặt, nước ngầm q trình thủy văn Đây thơng số hình thành nên cân nước khu vực Trong nghiên cứu này, tác giả thiết lập lại số số liệu đầu vào phát triển thêm mơ hình MIKE SHE DHI Hungary thiết lập, nhằm mục đích mơ q trình thủy văn cân nước lưu vực Dong-ér, thuộc Hungary Quá trình thực nghiên cứu sử dụng phương pháp one-at-time để phân tích độ nhạy số liệu đầu vào Cụ thể số liệu nhiệt độ đầu vào gia tăng phù hợp với kịch biến đổi khí hậu số liệu đầu vào khác giữ nguyên theo mơ hình sở năm 2018 Bài báo xem xét ảnh hưởng nhiệt độ tăng đến thông số thủy văn thành phần cân nước, dựa kết mô kịch Kết mô năm 2018 số liệu tham khảo để làm sở so sánh độ nhạy của việc biến đổi số liệu nhiệt độ đầu vào đến thông số thủy văn thành phần cân nước bối cảnh biến đổi khí hậu Từ khóa: Kênh Dong-ér; MIKE SHE; Phân tích độ nhạy; Cân nước Abstract Effects of temperature increase by climate change on hydrological parameters and water balance composition in Dong-ér basin, Hungary Extreme weather and climate changes are emerging more frequently in Central Europe, Hungary, and in the near future the increase in prolonged droughts, high-intensity precipitation events and the temporal variations of precipitation are expected, which may increase the magnitude of local water damages As a result of climate change, these extreme weather events will be more frequent, however it is difficult to predict them, as until now insufficient number of observations are available on smaller watercourses and on refined territorial water balances For the future assessment of the environmental and economic impacts of climate change, it is essential to explore the integrated relationship of evapotranspiration, runoff, infiltration, surface and subsurface waters, and other hydrological processes, which can fundamentally describe regionally the water management conditions In this research, an earlier study on the catchment area of the main canal of the Dong-ér Brook is pursued to continue the development of the MIKE SHE model in a more complex manner During the research process, one-at-time method was used to analyze the sensitivity of the input data Specifically, the input temperature data has been increased in line with climate change scenarios and other input data not change according to the base simulation of Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững 533 2018 Based on simulation results of scenarios, considering the influence of temperature rise on hydrological parameters and water balance components Using data from 2018 as reference, the sensitivity of the changes in daily mean temperature has been assessed to estimate the effects of the future climate change on hydrological parameters and water balance components Keywords: Dong-ér Brook; MIKE SHE; Sensitivity analysis; Water balance Giới thiệu Quản lý tài nguyên nước đóng vai trị quan trọng việc giảm nhẹ ảnh hưởng biến đổi khí hậu cực đoan Hiện tượng khí hậu cực đoan nhận biết hình thức khơ hạn với thời gian dài hơn, với cường độ lớn với tần suất thường xuyên Theo kết khảo sát Ladányi cộng (2014) miền Nam Hungary, tượng khô hạn, thiếu nước gây thiệt hại lớn so với tượng ngập úng Cụ thể biến đổi khí hậu, đặc biệt khơ hạn ảnh hưởng trực tiếp đến phát triển bền vững ngành nông nghiệp, qua đó, ảnh hưởng gián tiếp đến an ninh lương thực (Singh, 2014, Ladányi, 2010) Theo quan sát thực tế Hungary, 10 năm lại có 04 năm xảy tượng khô hạn nghiêm trọng Như vậy, cần phải có phương pháp tiếp cận mang tính tổng hợp để đánh giá từ nhiều góc nhìn khác mối quan hệ phức tạp trình thủy văn thành phần cân nước Để thực điều cần phải có giải pháp ứng dụng cơng cụ mơ hình tổng hợp Nhiệm vụ ghi rõ Kế hoạch Quản lý lưu vực sông Tisza ban hành Tổng cục Quản lý nước Hungary (OVF) năm 2015 Với phát triển ngành công nghệ thông tin, ứng dụng Hệ thống thông tin địa lý (GIS) mơ hình thủy văn tạo điều kiện cho nhà nghiên cứu có nhiều cơng cụ để thực mơ Trong số mơ hình mơ hình MIKE SHE lên mơ hình tổng hợp, cơng cụ hiệu cơng tác quản lý tổng hợp tài ngun nước Mơ hình cho thấy tính hiệu việc mơ mối quan hệ tương tác nước mặt nước ngầm (Graham and Butts, 2005) Nhóm tác giả Lu cộng (2006) có nhận xét tích cực liên quan đến tính hiệu mơ hình việc mơ dịng chảy tồn lưu vực Mơ hình MIKE SHE sử dụng rộng rãi nghiên cứu xác định ảnh hưởng biến đổi sử dụng đất đến cân nước vùng chưa bão hòa vùng bão hòa (Asadusjjaman and Farnaz, 2014) Một số nghiên cứu lại kết hợp MIKE SHE GIS để mơ q trình thủy văn lưu vực khác (Paparrizos Maris, 2015; Právetz cộng sự, 2015) đánh giá ảnh hưởng biến đổi sử dụng đất biến đổi khí hậu đến nước ngầm hệ sinh thái Keilholz cộng (2015) thực Mơ hình MIKE SHE áp dụng Hungary lưu vực kênh Fehértó-majsa Lưu vực nằm cách lưu vực Dong-ér phía Nam Benyhe cộng thực năm 2015 Nhóm tác giả ứng dụng ba mơ hình bao gồm BUDYKO, HEC-RAS MIKE SHE để mơ q trình thủy văn nhận định rằng, mơ hình MIKE SHE hiệu nhiều so với hai mơ hình Nhóm tác giả Leeuwen (2016) kết hợp MIKE SHE hình ảnh viễn thám để kiểm định mơ hình cho thấy tính hiệu mơ hình MIKE SHE cơng tác dự báo mức độ, phạm vi thời gian diễn ngập úng Nhóm nghiên cứu Nagy cộng (2019) xây dựng mơ hình MIKE SHE để mô ngập úng, lập đồ ngập úng cho lưu vực kênh Dong-ér Các nghiên cứu cho thấy nhược điểm mơ hình cần nhiều số liệu đầu vào với chất lượng tốt Tại lưu vực kênh Dongér, số liệu khơng đầy đủ, khơng có số liệu mực nước lưu lượng nước kênh số thống số khác, như: số LAI, độ sâu rễ, số độ dẫn thủy lực không kiểm định Nagy cộng (2019) hiệu chỉnh mơ hình dựa số liệu thực đo mực nước ngầm vào mùa xuân năm 2015 năm 2018 Kết chênh lệch kết mô số 534 Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững liệu thực đo 45 cm - kết nói tốt bối cảnh thơng số khơng kiểm định có nhiều rủi ro Kết cho thấy, mơ hình xây dựng cho lưu vực kênh Dong-ér hiệu chỉnh cách tương đối xác kế thừa để tính tốn phân tích khía cạnh khác thủy văn mối quan hệ nước mặt nước ngầm Nhu cầu hiểu rõ cân nước lưu vực có nhiều ảnh hưởng hạn hán cần thiết đóng vai trị quan trọng công tác quản lý, phân bổ, sử dụng hiệu bền vững tài nguyên nước (The EU Water Framework Directive, 2004; OVF, 2009) Hai tác giả Sipos Právecz (2014) nhận định rằng, với công cụ mơ hình cân nước cục tính tốn đánh giá tính hiệu việc sử dụng nước bền vững mùa khơ tích trữ nước mùa mưa Một lợi mơ hình MIKE SHE bao gồm cơng cụ tính tốn cân nước tổng hợp tồn diện cục cho toàn lưu vực thời điểm mô Số liệu đầu kết cân nước bao gồm dòng chảy, biến đổi lượng nước lưu trữ hệ thống biến đổi cán cân nước (DHI, 2017) Giá trị thơng số đóng vai trị quan trọng việc đánh giá cân nước lưu vực kênh Dong-ér Có nhiều yếu tố tác động trực tiếp gián tiếp đến cân nước lưu vực Đầu tiên phải kể đến lượng mưa nhiệt độ Đây hai yếu tố đóng vai trị chủ lực việc vận hành vịng tuần hồn nước có ảnh hưởng lớn đến cân nước lưu vực Câu hỏi đặt là, thông số thủy văn thành phần cân nước thay đổi số liệu đầu vào bị biến đổi ảnh hưởng biến đổi khí hậu? Theo nghiên cứu Hamby (1994) Lenhart cộng (2002) biến đổi số liệu đầu vào ảnh hưởng đến kết mơ hình phải áp dụng phương pháp phân tích độ nhạy xem thông số thủy văn hay thành phần cân nước bị ảnh hưởng nhiều Mục tiêu nghiên cứu mô điều kiện khí hậu khác áp dụng phương pháp phân tích độ nhạy để đánh giá ảnh hưởng biến đổi số liệu đầu vào đến thông số thủy văn thành phần cân nước Có nhiều phương pháp phân tích độ nhạy, nhiên, phần lớn phương pháp có độ phức tạp khác nhau, khó áp dụng cho lưu vực kênh Dong-ér Tác giả sử dụng phương pháp đơn giản phương pháp one-at-a-time Với phương pháp thơng số biến đổi, cịn lại thông số khác giữ nguyên Nhược điểm phương pháp phân tích độ nhạy phân tích cục điểm nghiên cứu áp dụng cho toàn lưu vực (Hamby, 1994) Tuy nhiên, với việc mô linh hoạt không gian thời gian mơ hình MIKE SHE nhược điểm khắc phục Để có nhìn tổng hợp tồn diện q trình thủy văn diễn lưu vực nghiên cứu, tác giả kế thừa số liệu mô hình MIKE SHE, Nagy cộng xây dựng năm 2019 cho lưu vực kênh Dong-ér Tuy nhiên, để áp dụng phương pháp one-at-a-time, tác giả điều chỉnh, thay đổi, cập nhật thông số thiết lập mơ đun tính tốn cân nước lưu vực để đánh giá tổng hợp biến đổi nhiệt độ có ảnh hưởng đến thơng số thủy văn thành phần cân nước Để thực điều cần phải (1) Xem xét tính phù hợp, xác định ưu nhược điểm mơ hình MIKE SHE cho lưu vực nghiên cứu; (2) Mô trình thủy văn biến đổi cân nước bối cảnh số liệu nhiệt độ đầu vào thay đổi; (3) Các kết mô so sánh với để xác định độ nhạy thông số thủy văn thành phần cân nước Như vậy, mục tiêu nghiên cứu đánh giá ảnh hưởng biến đổi khí hậu, đặc biệt nhiệt độ đến thông số thủy văn thành phần cân nước tương lai lưu vực Dong-ér, Hungary Vùng nghiên cứu Lưu vực kênh Dong-ér nằm khu vực sườn cát, thuộc đồng sơng Danube sơng Tisza, cách biên giới phía Nam Hungary 50 km Lưu vực có diện tích khoảng 2127 km2 Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững 535 Hình 1: Khu vực nghiên cứu Lưu vực kênh Dong-ér bao gồm nhiều kênh nhỏ khác xây dựng nhằm mục đích thu gom nước chảy tràn đổ vào dịng kênh kênh Dong-ér Kênh Dong-ér bắt nguồn từ thị trấn Kiskunhalas đổ vào sơng Tisza khu vực thị trấn Baks (Hình 1) Quá trình khảo sát thực địa cho thấy, đoạn đầu dịng kênh có lưu lượng nước mặt thấp, chí lịng kênh khơ với nhiều loại thực vật khác Như vậy, nguồn nước mặt khu vực khơng đủ số lượng, vậy, lượng nước ngầm sử dụng (Kozák, 2020) Ở đoạn gần thị trấn Kiskunhalas bắt đầu thấy có nhiều nước kênh Lưu lượng nước kênh Dong-ér vào khoảng - m3/s, giai đoạn khơng có mưa Trong đó, vào thời điểm mùa xn lưu lượng nước kênh đạt tới 20 - 30 m3/s Tuy nhiên, vào tháng mùa hè kênh lưu vực khơ hồn tồn (K&K Mérnưki Iroda Kft, 2013) Hướng gió chủ đạo lưu vực có hướng Tây Bắc, với tốc độ gió trung bình vào khoảng - m/s Ở lưu vực có đặc điểm đụn cát gió có ảnh hưởng đến q trình kiến tạo địa hình Trong đó, địa hình lại tạo nên hệ thống dịng chảy mặt Dịng kênh Dong-ér có hướng chảy từ phía Tây Nam phía Đơng Bắc, đó, kênh phụ có hướng chảy từ phía Tây Bắc đổ vào kênh Dong-ér Đông Nam Trong mùa mưa, nước ngầm xuất vùng đất trũng tạo nên tượng ngập úng cục (Sipos Právecz, 2014) Theo nghiên cứu Trung tâm Khí tượng Hungary (2018), dựa hai mơ hình khí hậu khu vực hai kịch biến đổi khí hậu nhiệt độ trung bình Hungary tăng từ - 0C cuối kỷ chạm đến mức tăng lên 0C sớm dự kiến Lưu vực kênh Dong-ér khu vực chịu ảnh hưởng mạnh mẽ tượng khô hạn (Sipos and Právecz, 2014) Cũng theo số liệu Trung tâm Khí tượng Hungary (OMSZ) số liệu từ 25 trạm quan trắc khí tượng giai đoạn từ năm 2010 đến năm 2018, lượng mưa trung bình năm khu vực nghiên cứu vào khoảng 611 mm, vào thời điểm cực đoan chạm tới lương mưa 842 mm (trong năm 2014) thấp vào khoảng 203 mm (trong năm 2000) Trong giai đoạn từ 2000 2018, nhiệt độ trung bình tháng thấp - 5,2 0C, đo vào tháng 02 năm 536 Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững 2012 tháng 01 năm 2017, nhiệt độ cao + 24,5 0C, đo tháng 08 năm 2018 Độ dày trung bình tuyết tháng mùa Đông vào khoảng 18 - 22 cm Theo tính tốn Ủy ban liên Chính phủ biến đổi khí hậu (IPCC), nhiệt độ trung bình Trái đất tăng 1.5 0C, giai đoạn từ 2030 - 2052 (IPCC, 2018) Tại lưu vực nghiên cứu, diện tích đất trồng trọt chiếm đến ~41 %, sau diện tích đồng cỏ, chiếm ~13 %, rừng rộng chiếm ~10 % tổng diện tích lưu vực, cịn lại mục đích sử dụng đất khác, mơ hình canh tác khác nhau, khu dân cư Dựa đồ đất nông nghiệp, loại đất lưu vực đa dạng Loại đất đặc trưng loại đất cát, Chernozem đất nhiễm mặn Đây loại đất có độ thấm cao Từ năm 1970 có xu hướng suy giảm lượng mưa khu vực Điều dẫn đến việc giảm mực nước ngầm khu vực đồng hai sông Danube Tisza (Fehér, 2019) Sau năm mưa nhiều, tình trạng thiếu ngầm lưu vực khôi phục Trên thực tế, đơi mực nước ngầm tăng đến mức có hại, gây tượng ngập úng cục (Szatmári and van Leeuwen, 2013) Số liệu phương pháp nghiên cứu 3.1 Các số liệu sử dụng mô hình 3.1.1 Số liệu khí tượng Dữ liệu khí tượng quan trọng chu kỳ thủy văn lượng mưa, nhiệt độ thoát nước Tại khu vực nghiên cứu, số liệu mưa quan trắc 25 trạm khí tượng giai đoạn 2010 - 2018 Cục Quản lý nước vùng hạ lưu sơng Tisza (ATIVIZIG) cung cấp Theo đó, lượng mưa trung bình tháng khu vực nghiên cứu dao động khoảng 46 mm Vào năm khơ hạn 2000, đo khoảng 29 mm Szeged (OMSZ, 2018a), số năm mưa nhiều năm mưa 2014 với lượng mưa 89 mm Lưu trữ tuyết trạng thái vật lý lượng mưa khác thời gian chảy tràn tham gia vào chu trình nước chậm Nhiệt độ yếu tố đẩy nhanh q trình chảy tràn hội tụ nước thơng qua q trình tan chảy Nhiệt độ ảnh hưởng trực tiếp (bốc nước ao, bốc đất) gián tiếp (thoát nước thảm thực vật) đến bốc thoát nước Nhiệt độ biến đổi khiến số yếu tố khí hậu quan trọng thay đổi theo, yếu tố ảnh hưởng lớn đến cân nước khu vực Dữ liệu nhiệt độ cho giai đoạn từ năm 2014 đến năm 2018 lấy từ sở liệu Đánh giá Dữ liệu Khí hậu châu Âu (European Climate Assessment and Dataset, ECAD) Dữ liệu nhiệt độ trung bình hàng ngày năm 2000 lấy từ sở liệu OMSZ (2018a) Bốc thoát nước thơng số khó xác định Khơng có số liệu thực đo số liệu tính tốn khu vực nghiên cứu Để tính lượng bốc nước, Tổ chức Lương thực Nông nghiệp Liên Hợp Quốc (FAO) khuyến nghị kết hợp phương trình FAO - 56 PenmanMonteith (Allen et al., 1998) Do thiếu liệu đầu vào, việc sử dụng phương trình FAO - 56 Penman-Monteith bị hạn chế Giá trị bốc thoát nước tiềm (ETP) ước tính dễ dàng đơn giản hóa hàm lũy thừa sau (Fiala cộng sự, 2018): = ETP ET = e0.07T ref Trong đó: ETref giá trị bốc thoát tham chiếu; T nhiệt độ Celsius Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững 537 3.1.2 Số liệu thảm phủ thực vật Trong mơ hình xác định 23 loại hình sử dụng đất khác Mỗi loại hình thảm phủ thực vật có đặc điểm khác khả tán ngăn nước, bốc nước, nước, hình thành dịng chảy khả giữ nước rễ Mơ hình sử dụng Corine Land Cover (EEA, 2018) làm tiêu chuẩn thực tế sử dụng đất quan trắc thảm phủ thực vật cấp độ châu Âu (Feranec, 2016; Aune - Lundberg Geir-Harald, 2020) Chỉ số diện tích Lá (LAI) xác định kết hợp hình ảnh MODIS (Myneni cộng sự, 2015) giá trị độ sâu rễ tham chiếu từ kết nghiên cứu cho loại thảm phủ khác (Nagy cộng sự, 2019) 3.1.3 Số liệu đất Việc xác định xác đặc trưng đất (hàm lượng nước đất khơng bão hịa bão hịa) quan trọng việc xây dựng mơ hình, Rakonczai cộng (2011) Farsang (2014) mối quan hệ tương tác đất cân nước Các thông số đất mơ hình xây dựng dựa Cơ sở liệu thủy lực đất 3D, với độ phân giải 250 m châu Âu (Tosth cộng sự, 2017) Nghiên cứu thiết lập dựa lớp đất sâu m sở liệu này, số liệu đáng tin cậy tốt Tác giả cập nhật thông số vào mơ hình: hàm lượng nước bão hịa, độ dẫn thủy lực bão hòa đường cong giữ ẩm Các thông số lớp địa chất sâu ước tính dựa 13 hồ sơ lỗ khoan, ATIVIZIG cung cấp 3.1.4 Số liệu địa hình Các điều kiện địa hình ảnh hưởng đáng kể đến dòng chảy mặt đất dòng chảy ngầm, mơ hình số độ cao (DEM) xác quan trọng cho mơ hình tối ưu Nghiên cứu dựa DEM, có độ phân giải 10 × 10 m Để chạy mơ hình MIKE SHE, cần phải chuyển đổi liệu DEM thành tệp điểm lưới không gian (.dfs2) Tệp điểm tạo ArcGIS chuyển đổi thành tệp dfs2 hộp công cụ MIKE Zero (DHI, 2017) Số liệu lưới khơng gian sử dụng làm liệu đầu vào 3.1.5 Số liệu hệ thống kênh Trong mơ hình mạng lưới nước ATIVIZIG cung cấp phân thành hai nhóm Nhóm bao gồm mạng lưới nước từ thời kỳ trước, cung cấp liệu bình đồ, liệu mặt cắt liệu thiết lập dạng điểm GIS 3D Còn lại loại kênh khác ATIVIZIG quản lý Trong trường hợp này, số liệu bình đồ cung cấp dạng Polyline shapefile, liệu mặt cắt ngang xác định theo trắc dọc Những liệu thiết lập mô đun thủy lực chuyên dụng MIKE Hydro Mơ đun cho phép lập mơ hình hóa mạng lưới sơng liên kết với mơ hình MIKE SHE để mơ tồn diện mơ hình thủy văn lưu vực nghiên cứu 3.1.6 Dữ liệu nước ngầm Nước ngầm có mực nước ổn định bị ảnh hưởng yếu tố bên ngồi so với vùng khơng bão hịa kênh Số liệu nước ngầm ATIVIZIG cung cấp từ 33 trạm đo nước ngầm, với tọa độ cụ thể từ năm 2010 đến năm 2018 Mực nước ngầm thay đổi theo địa hình đổ vào sông Tisza Theo số liệu quan trắc giai đoạn 2010 - 2018, mực nước ngầm cao thường xuất vào tháng 04 - 05, thấp vào tháng 10 - 12 3.1.7 Các lớp tính tốn Để mơ q trình dịng chảy xảy vùng bão hịa, cần phải xác định lớp tính tốn Sự phân bố khơng gian của tầng chứa nước thấp có độ sâu - 75 m Độ dày lớp 538 Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững tối thiểu áp dụng với giá trị m, nhằm ngăn chặn lớp có độ dày mỏng Để xác định điều kiện biên ngồi, mơ hình áp dụng Fixed head type, dựa tệp dfs2, phân bố theo không gian thay đổi theo thời gian, trích xuất từ mực nước ngầm Mực nước ngầm thiết lập dựa liệu cao độ nước ngầm hàng ngày, từ năm 2010 đến năm 2018 Mực nước ngầm chia thành lưới với độ phân giải 200 m × 200 m Sau đó, mơ hình nội suy từ tệp.dfs2 theo thời gian không gian đến biên đầu bước thời gian (DHI, 2017a) 3.1.8 Dữ liệu hiệu chỉnh Mơ hình MIKE SHE xây dựng cho lưu vực Dong-ér thiếu nhiều liệu, mực nước kênh lưu lượng để hiệu chỉnh kiểm định đầy đủ Mơ hình tập trung hiệu chỉnh dựa liệu nước ngầm cho mùa xuân năm 2018, Nagy et al (2019) So sánh kết mô liệu đo được, phương sai khoảng 45 cm đạt yêu cầu bối cảnh thiếu số liệu hiệu chỉnh yếu tố khơng chắn (như đồ địa hình, đặc điểm thảm thực vật đặc điểm đất) Có nhiều thơng số cần hiệu chỉnh để mơ hình hồn thiện hơn, ví dụ, độ dẫn thủy lực bão hịa để hiệu chuẩn dịng chảy khơng bão hịa độ dẫn thủy lực, suất riêng, lưu trữ cụ thể để hiệu chuẩn dòng nước ngầm Tuy nhiên, thiếu liệu mục tiêu nghiên cứu không đề cập đến 3.2 Các số liệu đầu mơ hình Kết đầu mơ hình MIKE SHE phụ thuộc vào q trình mô chọn Theo Hướng dẫn sử dụng MIKE SHE (2017), mơ hình lưu kết 03 nhóm tệp sau: (1) Tệp ASCII, danh mục tất tệp đầu (.sheres) mô phỏng; (2) Tệp đầu nhị phân chứa tất thông tin tĩnh mô (.frf); (3) Các tệp tạo theo chuỗi thời gian, đơn giản dfs0, xác định theo 02 chiều dfs2 theo chiều dfs3 Trong nghiên cứu này, kết định dạng sheres kết chuỗi thời gian (dfs0, dfs2 dfs3) cho cân nước kiểm tra Bảng Mơ cân nước thực với mô đun riêng biệt cách tính tốn liệu lưu trữ tệp sheres Cách hiểu thông thường cho rằng, tất nước chảy vào hệ thống dương (+) tất dòng chảy nước âm (-) Dự trữ hệ thống nước dương (+), trường hợp gia tăng lưu trữ Thay đổi cân nước cho kết dương (+), nghĩa là, thay đổi lưu trữ cộng với tổng dòng chảy nhỏ tổng dòng chảy vào (∆Lưu trữ + Dòng chảy < Dịng chảy vào) Các thơng số thủy văn thành phần cân nước xem xét Bảng Bảng Bảng Các thông số thủy văn đầu sử dụng cho nghiên cứu Các thông số Đơn vị Bốc nước thực mm/ngày Thoát nước thực tế mm/ngày Bốc đất thực tế mm/ngày Độ sâu nước bề mặt Dòng chảy mặt theo hướng x y Dòng thấm vào vùng chưa bão hòa m m3/s mm/ngày Ghi Sự thoát nước thực tế phụ thuộc vào lượng nước sẵn có độ sâu rễ số diện tích (LAI) Sự nước phần lớn qua lỗ khí Sự bốc thực tế đất phụ thuộc vào yếu tố đặc điểm đất, sử dụng đất, thảm thực vật Nó lượng nước bề mặt Dịng chảy trục ngang trục dọc Dòng chảy thẳng đứng xuống có giá trị dương Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững 539 Các thông số Đơn vị Sự thiếu hụt vùng chưa bão hòa mm Lượng nước trung bình vùng rễ Lượng nước trung bình vùng chưa bão hòa Dòng chảy ngầm theo hướng x y m3/s Ghi Sự thiếu vùng chưa bão hịa lượng khơng khí có vùng chưa bão hịa Do đó, thiếu hụt giảm dần có nghĩa khơng khí bị nước đẩy ngồi nhiều hơn, đất dần trở nên ẩm ướt Nó phụ thuộc vào lượng mưa, đặc điểm đất, dòng nước ngầm đặc điểm thảm thực vật Ngược lại với thiếu hụt vùng khơng bão hịa Dòng nước ngầm chảy theo trục ngang trục dọc Nguồn: DHI, 2017 Bảng Các thành phần cân nước đầu sử dụng cho nghiên cứu Tên thành phần Lượng mưa Bốc thoát nước Các dòng chảy Dòng chảy mặt biên vào Dòng chảy mặt biên Dịng chảy mặt vào sơng Dịng chảy đất biên vào Dòng chảy đất biên Dòng chảy sở vào sơng Dịng chảy sở từ sông Các kho lưu trữ nước Lượng lưu trữ tán Lượng lưu trữ tuyết Lượng lưu trữ đất bao gồm lưu trữ vùng chưa bão hòa lưu trữ vùng bão hòa Biến đổi cân nước Đơn vị mm mm mm mm mm Nguồn: DHI, 2017 3.3 Q trình mơ Trong giai đoạn phân tích độ nhạy, liệu đầu vào chia thành hai nhóm, phụ thuộc vào thay đổi theo thời gian chúng Nhóm không thay đổi khoảng thời gian dài, bao gồm: địa hình, mạng lưới nước (sơng hồ), tính chất đất, đặc điểm địa chất Nhóm có thay đổi không đáng kể không thay đổi thời gian dài khơng có tác động người Nhóm cịn lại bao gồm số liệu thay đổi nhanh chóng, thay đổi phút, giờ, ngày hay theo mùa Đó lượng mưa ngày, nhiệt độ trung bình ngày, lượng nước tham chiếu lớp phủ đất (bao gồm thảm thực vật) Sau đó, ba kịch mơ MIKE SHE xây dựng, bao gồm mô sở (BS), mơ SIM1, SIM2 (Hình 2) Mơ sở (BS) sử dụng để mô trình thủy văn tổng hợp cân nước lưu vực Dong-ér vào năm 2018 Lý năm 2018 chọn để chạy mơ hình sở là: (1) Mơ hình hiệu chỉnh Nagy cộng (2019) theo mực nước ngầm đo vào năm 2018; (2) Dữ liệu lớp phủ thực vật lấy từ Corine Land Cover năm 2018; (3) Các thông số van Genuchten khả giữ ẩm độ dẫn thủy lực tác giả 540 Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững cập nhật từ Cơ sở liệu thủy lực đất 3D năm 2017 (số liệu nhất); (4) Lượng mưa đo năm 2018 gần với mức trung bình nhiều năm Kết mô sở sử dụng làm giá trị tham chiếu để so sánh với kịch mô khác Phân tích độ nhạy nhiệt độ xác định ảnh hưởng gia tăng nhiệt độ đến thay đổi kết thông số thủy văn thành phần cân nước Dựa kết IPCC (2018) cho thấy, đến năm 2100, nhiệt độ trung bình tồn cầu tăng lên + 1,5 0C so với giai đoạn 1986 - 2005 theo kịch RCP 2.6 theo nghiên cứu mơ hình khí hậu OMSZ (2018b), đến năm 2060, nhiệt độ tăng lên + 0C theo kịch cực đoan RCP 8.5 Trên sở này, số liệu nhiệt độ mơ hình tham chiếu 2018 (BS) tăng lên + 1,5 0C (SIM1) + 2,2 0C (SIM2) Trong khi, thông số khác số liệu mưa, thảm phủ thực vật khơng thay đổi Các thơng số có mối quan hệ trực tiếp với nhiệt độ gia tăng với mức độ khác thơng số bốc (bốc nước mặt, bốc đất, bốc lưu trữ cây, bốc tuyết, thoát nước thực vật) Sau đó, kết mơ ba mơ hình phân tích đánh giá Hình 2: Sơ đồ nghiên cứu Kết Do hạn chế thời gian khả lưu trữ kết thủy văn trích xuất từ ô lưới gần khu vực thị trấn Tiszaalpár, khơng trích xuất từ tồn lưu vực Kết phân tích độ nhạy thơng số thủy văn tác động nhiệt độ tăng + 1,5 C + 2,2 0C thể Hình Sự nước thực tế tăng 41 %, ảnh hưởng + 1,5 0C tăng 67 % nhiệt độ hàng ngày tăng + 2,2 0C Giá trị thiếu hụt vùng khơng bão hịa cho thấy độ nhạy cao, cụ thể, mức tăng 37 % 57 % Nhiệt độ tăng thêm + 1,5 0C làm cho thoát nước thực tế tăng % so với bốc đất Sự thay đổi tỷ lệ nghịch lớn dòng chảy đất liền theo hướng x, giảm 21 % + 1,5 0C giảm 28 % nhiệt độ tăng + 2,2 0C (Hình 3) Vì vậy, thơng số thủy văn nhạy cảm gia tăng nhiệt độ Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững 541 Hình 3: Các thơng số thủy văn biến đổi nhiệt độ tăng lên + 1,5 0C + 2,2 0C, so với giá trị tham chiếu vào ngày 10 tháng năm 2018 (trục Y = 0) Trong số thông số thủy văn Hình 3, thâm hụt vùng khơng bão hịa, thoát nước thực tế thẩm thấu vào vùng khơng bão hịa lên yếu tố chính, có ảnh hưởng lớn đến vịng tuần hồn thủy văn khu vực Nhiệt độ tăng có ảnh hưởng lớn đến giá trị thâm hụt khơng bão hịa Ngồi ra, nhiệt độ tăng làm tăng nước thực tế làm giảm thẩm thấu cách tự nhiên Các thơng số cịn lại có giá trị ngày nhỏ 1,0 mm Trong đó, lượng mưa vào ngày 10 tháng năm 2018 không Điều cho thấy, lượng nước mơ ngày có nguồn gốc từ nguồn nước tích trữ lớp đất từ giai đoạn trước Lượng mưa chạm đến bề mặt đất, phần lớn thấm vào đất với tốc độ chậm phần lại chảy bề mặt theo hướng x (theo hướng dốc) Hình 4: Cao độ địa hình cao độ nước ngầm Cuộc khảo sát thực địa vào tháng năm 2020 cho thấy, nước bắt đầu có kênh Dongér, khu vực bên thị trấn Kiskunhalas lượng nước kênh tăng dần phía sơng 542 Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững Tisza Dựa kết này, nhận định khẳng định nguồn gốc dòng chảy lưu vực Dong-ér lượng mưa dịng chảy từ bên ngồi vào Câu hỏi đặt kênh phạm vi nghiên cứu lấy nước từ đâu? Từ thị trấn Kiskunhalas trở xuôi xuống, nguồn cung cấp nước kênh xác định từ ba nguồn Đầu tiên nguồn nước ngầm từ khu vực có độ cao chảy xuống khu vực thấp tác động trọng lực (Hình 4) Dựa Hình theo ngun tắc bình thơng nhau, nước ngầm liên tục cung cấp cho vùng nước mặt từ hướng Tây Bắc, Tây Nam phía Tây Ở khu vực có độ chênh lệch, mực nước ngầm cao độ địa hình khơng lớn, chí mực nước ngầm cao cao độ đáy kênh khu vực gần thị trấn Bugc, Tisszaalpár Baks nước ngầm liên tục bổ sung nước cho kênh gây nên ngập úng vùng đất trũng, thị trấn ven sơng Tisza Như vậy, khẳng định, nước ngầm nguồn cung cấp nước cho kênh lưu vực nghiên cứu Ngoài nguồn nước ngầm cịn hai nguồn nước, lượng nước mưa trực tiếp bổ sung cho kênh nguồn nước thải từ khu dân cư, nhiên, theo quan sát thực tế hai lượng nước không đáng kể Để xem xét mối quan hệ mối tương quan thành phần cân nước cho toàn lưu vực, tác giả sử dụng mơ-đun tính tốn cân nước mơ hình MIKE (Hình Bảng 3) Hình 5: Các thành phần cân nước biến đổi nhiệt độ tăng lên + 1,5 0C + 2,2 0C so với giá trị tham chiếu năm 2018 (trục Y = 0) Bảng Các thành phần cân nước biến đổi nhiệt độ tăng lên + 1,5 0C + 2,2 0C so với giá trị tham chiếu năm 2018 Các thông số Lượng mưa Bốc nước Dịng chảy mặt từ bên ngồi chảy vào Dịng chảy mặt chảy bên ngồi Lưu trữ tán Lưu trữ tuyết Dòng chảy mặt bổ sung cho kênh Lưu trữ nước mặt Lưu trữ vùng khơng bão hịa Lưu trữ vùng bão hòa Giá trị tham chiếu năm 2018 (+) 589 (-) 795 (-) 1042 0 (-) 14 - 98 - 39 + 1,5 0C + 2,2 0C (+) 589 (+) 589 (-) 889 (-) 924 0 (-) 1036 (-) 1032 0 0 (-) (-) 12 11 - 153 - 168 - 55 - 63 Biến đổi (%) 12 -1 0 - 22 - 14 - 56 - 41 Biến đổi (%) 16 -1 0 - 33 - 21 - 71 - 62 Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững 543 Dòng chảy vào kênh Dòng chảy từ kênh Dòng chảy ngầm từ bên ngồi vào Dịng chảy ngầm chảy bên ngồi Dịng thấm bao gồm bốc Thốt nước bao gồm bốc thoát Biến đổi cân nước (-) 12 (+) (+) 1132 (-) 101 (+) 109 (-) 1170 (-) 357 (-) 11 (-) 10 (+) (+) (+) 1137 (+) 1140 (-) 97 (-) 95 (+) 64 (+) 43 (-) 1154 (-) 1147 (-) 506 (-) 554 -8 0 -4 - 41 -1 42 - 17 -6 - 61 -2 55 Lưu ý: Ký hiệu ngoặc hướng dòng nước Không nằm ngoặc thể giá trị thông số Dựa kết Bảng 3, thành phần dịng chảy ngầm từ bên ngồi chảy vào, dịng chảy mặt chảy bên ngồi bốc nước có ảnh hưởng lớn đến cân nước khu vực nghiên cứu Cân nước lưu vực Dong-ér phụ thuộc nhiều vào ảnh hưởng dòng chảy ngầm từ bên ngồi vào Sự nước tăng nhiệt độ tăng, làm giảm đáng kể lượng lưu trữ đất (bao gồm lưu trữ vùng khơng bão hịa vùng khơng bão hịa, từ -48,5 % đến - 66,5 %) Lưu trữ nước mặt đất giảm từ - 14 % xuống - 21 %, đó, dịng chảy vào kênh nguồn nước bổ sung cho kênh giảm theo Nhiệt độ tăng khơng làm thay đổi dịng chảy từ kênh, vậy, giá trị khơng thay đổi so với năm 2018 (tất mm) Kết luận, nhiệt độ tăng làm giảm khả giữ nước kênh, đặc biệt nơi có địa hình cao mực nước ngầm Trong ba điều kiện nhiệt độ, giá trị dịng chảy mặt từ bên ngồi chảy vào khơng, ranh giới lưu vực bề mặt xác định xác Sự thay đổi lưu trữ tán lưu trữ tuyết ba điều kiện nhiệt độ không, điều giá trị nhỏ, nhỏ 0,03 mm, nên làm trịn số khơng Dịng chảy ngầm từ bên chảy vào ổn định đóng vai trị quan trọng việc cân nước lưu vực Dong-ér Dòng chảy ngầm chảy bên ngồi có xu hướng giảm từ - % xuống - % Một phần trình di chuyển thấm xuống, bổ sung cho nước ngầm phần cịn lại bốc thực vật Khi xem xét giá trị hai thông số thấm thoát bao gồm bốc thoát hơi, kết luận rằng, nhiệt độ tăng độ thấm giảm mạnh (- 41 % - 61 %), cịn độ giảm Sự thay đổi cân nước năm 2018 vào khoảng -3 57 mm, cụ thể 758 triệu m3 nước thoát khỏi hệ thống Sự gia tăng nhiệt độ + 1,5 0C + 2,2 0C làm giảm cân nước lưu vực Dong-ér, 42 % 55 %, tức là, giảm 316 triệu m3 418 triệu m3 nước so với năm 2018 Nhiệt độ tăng có tác động trực tiếp gián tiếp làm tăng lượng nước thất thoát bốc nước đó, cân nước lưu vực bị giảm đáng kể Dựa kết này, kết luận rằng, thành phần, lưu trữ đất, dòng thấm bao gồm bốc hơi, thay đổi cân nước, dịng chảy mặt bổ sung cho kênh, lưu trữ nước mặt thay đổi mức độ lớn so với lượng bốc thoát nước, tức lớn 12 % 16 % Vì vậy, thơng số nhạy cảm với biểu nhiệt độ ngày tăng biến đổi khí hậu Kết luận Các thơng số, như: thâm hụt vùng chưa bão hòa, dòng thấm vào vùng khơng bão hịa, bốc nước dịng chảy mặt, có ảnh hưởng lớn tuần hồn thủy văn có độ nhạy cao với điều kiện biến đổi khí hậu Lượng mưa nhanh chóng thấm đất lưu trữ tạm thời lớp đất phía thấm vào lớp sâu đất với thời gian chậm trễ phần lại chảy bề 544 Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững mặt theo hướng x (theo hướng dốc) Trong thời gian này, lượng lớn nước rời khỏi hệ thống thơng qua q trình nước Dựa mối quan hệ biểu biến đổi khí hậu thành phần cân nước kết luận, dịng chảy ngầm từ bên ngồi vào bốc thoát nước hai động lực hình thành điều hịa cân nước lưu vực Dong-ér Theo kết mô cho thấy, cân nước lưu vực Dong-ér xác định đáng kể dịng chảy ngầm từ bên ngồi vào cung cấp liên tục khoảng ~ 90 % lượng nước bề mặt Dịng chảy ngầm từ bên ngồi vào không cho thấy mối quan hệ tương tác với thay đổi nhiệt độ Trong đó, bốc nước phụ thuộc nhiều vào nhiệt độ Do xu hướng khí hậu nóng lên, khả giữ nước kênh lưu vực Dong-ér có xu hướng giảm dần Đặc biệt khu vực có mực nước ngầm sâu lòng kênh Các thành phần đóng vai trị quan trọng cân nước lưu vực Dong-ér theo thứ tự giảm dần sau: Dịng chảy ngầm từ bên ngồi vào; Dịng chảy mặt chảy bên ngồi; Bốc nước lượng mưa Tuy nhiên, số thành phần, lưu trữ nước đất (bao gồm lưu trữ vùng bão hịa khơng bão hịa), dịng thấm bao gồm bốc hơi, dịng chảy mặt bổ sung cho kênh, lưu trữ nước mặt bốc thoát nước nhạy cảm với biến đổi khí hậu Dựa thay đổi cân nước mơ hình, nhiệt độ tăng + 1,5 0C + 2,2 0C làm nhiều nước so với năm 2018 Điều làm cho cân nước lưu vực Dong-ér tình trạng thiếu nước ngày trầm trọng Kết nghiên cứu cho thấy, mơ hình MIKE SHE mơ đun tính tốn cân nước cơng cụ hữu ích để phân tích đánh giá ảnh hưởng biến đổi điều kiện khí hậu đến vịng tuần hồn thủy văn Tuy nhiên, mơ hình MIKE SHE có nhược điểm, như: tính độ thấm theo phương thẳng đứng, điều dẫn đến kết tính tốn khơng hồn tồn khu vực có độ dốc lớn Bên cạnh đó, việc cần nhiều liệu đầu vào với chất lượng tốt thay đổi điều kiện môi trường tương lai vấn đề cần cải thiện nên giá trị tính tốn thơng số khơng thể xác tuyệt đối cung cấp giá trị gần xu hướng biến đổi Kết nghiên cứu tài liệu tham khảo cho nhà quản lý bên liên quan khác việc quy hoạch, quản lý sử dụng hiệu nguồn tài nguyên nước tình trạng suy giảm số lượng bối cảnh biến đổi khí hậu diễn ngày phức tạp Lời cảm ơn: Tác giả xin cảm ơn Công ty DHI Hungary cấp phép sử dụng mơ hình MIKE SHE Xin cảm ơn Cục trưởng Cục Quản lý nước vùng hạ lưu sông Tisza (ATIVIZIG) Dr Kozák Péter đồng nghiệp lời khuyên hỗ trợ chuyên môn Tôi biết ơn Prof Dr Rakonczai János hỗ trợ động viên cho việc viết báo TÀI LIỆU THAM KHẢO [1] Allen, R G., Pereira, L S., Raes, D., Smith, M (1998) Crop evapotranspiration - guidelines for computing crop water requirements FAO Irrigation and drainage paper 56 Food and Agriculture Organization, Rome http://www.fao.org/3/x0490e/x0490e00.htm [2] Asadusjjaman, S., Farnaz, A (2014) Investigation of Water Balance at Catchment Scale using MIKE-SHE International Journal of Engineering and Computer Science 3, 8882 - 8887 [3] ATIVIZIG (2016) River Basin Management Plan - 20 Lower Tisza right bank planning subunit (in Hungarian) Online available at http://www.vizugy.hu/vizstrategia/documents/ C8F99A5D-2864-455E-85D6-7FB53AF0D179/VGT2 2_20_Also_Tisza_jobb_part_vegleges.pdf [4] Aune-Lundberg, L., Geir-Harald, S (2020) The content and accuracy of the CORINE Land Cover Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững 545 dataset for Norway International Journal of Applied Earth Observation and Geoinformation 96 (102266), 1-10 DOI: 10.1016/j.jag.2020.102266 [5] Benyhe, B., Právecz, T., Sipos, Gy (2015) Assessment of different hydrological modelling software on a lowland minor catchment In: Blanka, V., Ladányi, Zs (Ed.) Drought and Water Management in South Hungary and Vojvodina University of Szeged 279 - 289 [6] DHI (2017a) MIKE SHE Volume 1: User guide [7] https://manuals.mikepoweredbydhi.help/2017/Water_Resources/MIKE_SHE_Printed_V1.pdf [8] EEA (2018) Corine Land Cover 2018 European Environmental Agency https://www.eea.europa.eu/ data-and-maps/data/external/corine-land-cover-2018 [9] Farsang, A (2014) The potential pedological effects of aridification due to climate change In: Blanka, V., Ladányi, Zs (Ed.) Drought and Water Management in South Hungary and Vojvodina University of Szeged 117 - 126 [10] Fehér, Z Zs (2019) Large scale geostatistical modelling of the shallow groundwater time series on the Southern Great Hungarian Plain Two approaches for spatiotemporal stochastic simulation of a non-complete monitoring dataset PhD dissertation University of Szeged DOI: https://doi.org/10.14232/ phd.10122 [11] Feranec, J., Soukup, T., Hazeu, G., Jaffrain, G (2016) (Ed.) European landscape dynamics Corine land cover data CRC-Press, Boca Raton - 14 DOI: https://doi.org/10.1201/9781315372860 [12] Fiala, K., Barta, K., Benyhe, B., Fehérváry, I., Lábdy, J., Sipos, Gy., Győrffy L (2018) Operational drought and water scarcity monitoring system (In Hungarian) Hungarian Journal of Hydrology http:// publicatio.bibl.u-szeged.hu/17598/1/Fiala_et_al2018HidrologiaiKozlony.pdf [13] Graham, D N., Butts, M (2005) Flexible, integrated watershed modelling with MIKE SHE In: Singh, V.P., Frevert, D.K (Ed) In Watershed Models CRC Press 245-272 DOI: 10.1201/9781420037432 ch10 [14] Hamby, D M (1994) A review of techniques for parameter sensitivity analysis of environmental models Environmental Monitoring and Assessment 32 135 - 154 DOI: https://doi.org/10.1007/BF00547132 [15] IPCC (2018) Global Warming of 1.5 0C Thematic Reports https://www.ipcc.ch/sr15/ [16] K&K Mérnöki Iroda Kft (2013) Harmonized activities related to extreme water management events - especially flood, inland inundation and drought (in Hungarian) HUSRB/1203/121/145/01, Ref No.: T-51/2013 [17] Keilholz, P., Disse, M., Halik Ü (2015) Effects of Land Use and Climate Change on Groundwater and Ecosystems at the Middle Reaches of the Tarim River Using the MIKE SHE Integrated Hydrological Model MDPI, Water 3040-3056 DOI: 10.3390/w7063040 [18] Kozák, P (2020) Changes in surface runoff on the south-eastern slope of the Danube-Tisza Interfluve Sand Ridge in the context of climate change In: Farsang, A., Ladányi, Zs., Mucsi, L (Ed.) Climate change challenges - From global to local (in Hungarian) GeoLitera, 109 - 115 [19] Ladányi, Zs (2010) Climate change impact in a sample area of Danube-Tisza Interfluve In: Kiss, T (Ed.): Natural geographical processes and forms Natural Geography Studies of the 9th National Conference of Geographical Doctoral Students, 93 - 98 (in Hungarian) http://acta.bibl.u-szeged.hu/68110/1/2009_ termeszetfoldrajzi_folyamatok_es_formak.pdf [20] Ladányi, Zs., Blanka, V., Armenski, T., Stankov, U., Sipos, Gy (2014) Interviewing agricultural stakeholders and landowners In: Blanka, V., Ladányi, Zs (Ed.) Drought and Water Management in South Hungary and Vojvodina University of Szeged 339 - 359 546 Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững [21] Lu, J., Sun, G., McNulty, G S., Comerford, B N (2006) Evaluation and application of the MIKE SHE model for a cypress­pine flatwoods watershed in North Central Florida Interagency Conference on Research in the Watersheds 2, 63 - 74 [22] Myneni, R., Knyazikhin, Y., Park, T (2015) MCD15A2H MODIS/Terra+Aqua Leaf Area Index/ FPAR 8-day L4 Global 500m SIN Grid V006 NASA EOSDIS Land Processes DAAC DOI: https://doi org/10.5067/MODIS/MCD15A2H.006 [23] Nagy, Zs.,  Pálfi, G., Priváczkiné Hajdú, Zs., Benyhe, B (2019) Operation of canal systems and multi-purpose water management - Dong-ér catchment (in Hungarian) In: Ladányi, Zs., Blanka, V (Ed.) Monitoring, risks and management of drought and inland excess water in South Hungary and Vojvodina University of Szeged 83 - 96 [24] OMSZ (2018a) Observed changes in Hungary (in Hungarian) https://www.met.hu/eghajlat/ eghajlatvaltozas/megfigyelt_valtozasok/Magyarorszag/ [25] OMSZ (2018b) To the margin of the IPCC Thematic Report assessing a 1.5 degree global temperature rise (in Hungarian) https://www.met.hu/ismerettar/erdekessegek_tanulmanyok/index php?id=2334&hir=Az_IPCC_1,5_fokos_globalis_homerseklet-emelkedest_ertekelo_Tematikus_ Jelentesenek_margojara [26] OVF (2015) Tisza River Basin Management Plan by General Directorate of Water Management (in Hungarian) [27] OVF (2016) Effects of climate change, hydrometeorological extremes (in Hungarian) http://www.ovf hu/hu/korabbi-hirek-2/a-klimavaltozas-hatasai-hidrometeorologiai-szelsosegek [28] Paparrizos, S., Maris, F (2015) Hydrological simulation of Sperchios River basin in Central Greece using the MIKE SHE model and geographic information systems Applied Water Science 7, 591 - 599 DOI 10.1007/s13201-015-0271-5 [29] Právetz, T., Sipos, G., Benyhe, B., Blanka, V (2015) Modelling runoff on a small lowland catchment, Hungarian Great Plains Journal of Environmental Geography (1 - 2), 49 - 58 DOI: 10.1515/ jengeo-2015-0006 [30] Rakonczai, J., Farsang, A., Mezősi, G., Gál, N (2011) Theoretical background of inland excess water formation (in Hungarian) Hungarian Geographical Review 135 (4), 339 - 350 http://publicatio.bibl.uszeged.hu/5631/1/1959489.pdf [31] Singh, A (2014) Conjunctive use of water resources for sustainable irrigated agriculture Journal of Hydrology 519, 1688 - 1697 http://linkinghub.elsevier.com/retrieve/pii/S0022169414007318 [32] Sipos, Gy., Právecz T (2014) Identification of water retention areas on the Dong-ér catchment using GIS In: Blanka, V., Ladányi, Zs (Ed.) Drought and Water Management in South Hungary and Vojvodina University of Szeged 157 - 167 [33] Szatmári, J., van Leeuwen, B (2013) (Ed.) Inland Excess Water - Belvíz - Suvišne Unutrašnje Vode, Szeged University of Szeged Novi Sad, University of Novi Sad DOI: 10.13140/2.1.5143.3920 [34] The EU Water Framework Directive (2004) - DOI:10.2779/50903 [35] Tóth, B., Weynants, M., Pásztor, L., Hengl, T (2017) 3D Soil Hydraulic Database of Europe at 250 m resolution Hydrological Processes DOI: 10.1002/hyp.11203 https://eusoilhydrogrids.rissac.hu/ [36] Van Leeuwen, B., Právetz, T., Liptay Z Á., Tobak, Z (2016) Physically based hydrological modelling of inland excess water Carpathian Journal of Earth and Environmental Sciences 11 (2), 497 - 510 http:// publicatio.bibl.u-szeged.hu/17155/ Ngày chấp nhận đăng: 10/11/2021 Người phản biện: PGS.TS Phạm Quý Nhân Nghiên cứu chuyển giao, ứng dụng khoa học công nghệ sử dụng hợp lý tài nguyên, bảo vệ môi trường phát triển bền vững 547 ... cứu đánh giá ảnh hưởng biến đổi khí hậu, đặc biệt nhiệt độ đến thông số thủy văn thành phần cân nước tương lai lưu vực Dong- ér, Hungary Vùng nghiên cứu Lưu vực kênh Dong- ér nằm khu vực sườn cát,... cho lưu vực nghiên cứu; (2) Mô trình thủy văn biến đổi cân nước bối cảnh số liệu nhiệt độ đầu vào thay đổi; (3) Các kết mô so sánh với để xác định độ nhạy thông số thủy văn thành phần cân nước. .. giả điều chỉnh, thay đổi, cập nhật thông số thiết lập mơ đun tính tốn cân nước lưu vực để đánh giá tổng hợp biến đổi nhiệt độ có ảnh hưởng đến thơng số thủy văn thành phần cân nước Để thực điều