Trang 1 ỨNG DỤNG MƠ HÌNH VẬN ĐỘNG KHÔNG ỔN ĐỊNH MÔ PHỎNG HỆ THỐNG NƯỚC DƯỚI ĐẤT LƯU VỰC SÔNG RIO COBRE, JAMAICAĐỗ Thị Thùy Dung1, Uton Henry2, Yangxiao Zhou3, Bùi Du Dương41Trường Đại họ
ỨNG DỤNG MƠ HÌNH VẬN ĐỘNG KHƠNG ỔN ĐỊNH MÔ PHỎNG HỆ THỐNG NƯỚC DƯỚI ĐẤT LƯU VỰC SÔNG RIO COBRE, JAMAICA Đỗ Thị Thùy Dung1, Uton Henry2, Yangxiao Zhou3, Bùi Du Dương4 1Trường Đại học Tài nguyên Môi trường Hà Nội 2Water Resources Authority Jamaica 3IHE Delft Institute for Water Education 4Trung tâm Quy hoạch Điều tra tài nguyên nước Quốc gia Tóm tắt Nghiên cứu mơ dịng chảy nước đất theo thời gian phạm vi lưu vực Rio Cobre, Jamaica Mơ hình nhằm nghiên cứu biến động lịch sử tầng chứa nước với tác động khai thác nước đất tượng tự nhiên Mơ hình nước đất 3D sở phần mềm GMS giai đoạn (2000-2020) xây dựng hiệu chỉnh theo mực nước đất quan trắc lưu vực Rio Cobre Lượng nước sơng tính tốn dựa dòng chảy kiệt đo đạc Ở lưu vực Thượng Rio Cobre, tầng chứa nước cung cấp lượng lớn cho hệ thống sông suối chiếm 33,7 % tổng lượng nước tự nhiên Mực nước đất thập kỷ qua thay đổi tùy theo nguồn bổ cập tự nhiên cho nước đất, khơng có xu hướng giảm liên tục Lượng bổ cập cho nước đất từ mưa nguồn cung cấp chính, sau nguồn nước đất sơng chiếm tỷ lệ nhỏ tổng lượng xả thải (~26 %) Phần tích chứa giảm chủ yếu vào năm khô hạn Nghiên cứu cho thấy tầng chứa nước không bị khai thác mức, nhiên, tác động biến đổi khí hậu, nhu cầu sử dụng ngày tăng, nguồn nước sẵn có lâu dài khơng cịn trì đối mặt với nguy cạn kiệt Từ khóa: Nước đất; Mơ hình dịng chảy khơng ổn định; Lưu vực sông Rio Cobre, Jamaica Abstract Using a transient groundwater flow model to simulate groundwater dynamics in Rio Cobre basin, Jamaica A transient groundwater flow model was conducted in the Rio Cobre basin, Jamaica The primary objective was to investigate the historical changes in the basin’s aquifer, considering both anthropogenic impacts, such as groundwater exploitation, and natural events The 3D groundwater simulation tool (GMS) was utilized to accurately reproduce historical changes in the aquifer resulting from human activities and climate variations The transient groundwater model for the period 2000-2020 was constructed and calibrated using measured groundwater levels from the Rio Cobre basin The computed groundwater levels exhibited variations consistent with observed heads, while the computed discharges to rivers were found to correspond to measured low flows Notably, the Upper Rio Cobre basin demonstrated a substantial discharge of water into multiple springs, contributing to approximately 33.7 % of the total natural discharge Analysis of the groundwater levels over the past two decades revealed variations in response to natural recharge, with no discernible trend of continuous decline Precipitation infiltration emerged as the dominant recharge mechanism, followed by river leakage It is worth mentioning that abstraction accounted for a relatively small proportion (~26 %) of the total discharge The results indicate that the aquifer is not overdeveloped, however, local depletion occurs and was caused by intensive urban water supply in the lower alluvial aquifer The potential for further groundwater development from the limestone aquifer in the upper basin is large and can be considered With proper planning and management, groundwater can meet increased demand in the basin Keywords: Groundwater; Transient flow model; Rio Cobre basin, Jamaica Hội thảo Khoa học Quốc gia 2023 167 Đặt vấn đề Việc sử dụng mức nước đất để hỗ trợ phát triển nông nghiệp nhu cầu nước sinh hoạt báo cáo nhiều nơi giới làm cạn kiệt nguồn dự trữ nước đất vòng nhiều thập kỷ Khi phân bố không gian thời gian lượng mưa bị ảnh hưởng nóng lên tồn cầu, tốc độ thời gian thấm bị ảnh hưởng Sử dụng phát triển đất ngăn chặn khuếch đại phản ứng nước đất nóng lên tồn cầu Mở rộng sở hạ tầng nhân tạo dẫn đến tăng nhu cầu, trồng rừng dẫn đến tăng độ che phủ [1] Điều quan trọng phải quan sát dao động khoảng thời gian dài mực nước đất tầng chứa nước biến đổi khí hậu gây Cùng với xu hướng chung giới, vùng Caribe nói chung, Jamaica nói riêng phải đối mặt với vấn đề tương tự Nhiều nghiên cứu thảo luận khí hậu khu vực Caribe thay đổi [2, 3] Những thay đổi dẫn đến nhiều biến động chu trình thủy văn Jamaica, cụ thể, nhiệt độ báo cáo tăng lên lượng mưa bị suy giảm Do đó, số thiên tai cực đoan bão hạn hán kéo dài ghi nhận từ năm 1990 Những kiện thời tiết cho có tác động bất lợi đến tài nguyên nước, đặc biệt tài nguyên nước đất quốc đảo Tương tự nước phát triển khác, kinh tế Jamaica phụ thuộc chủ yếu vào tài nguyên thiên nhiên Nền kinh tế đóng góp du lịch thủy lợi hoạt động khai thác mỏ, thủy lợi bao gồm trồng mía, cà chua, dứa, bí ngơ tiêu thụ 75 % tổng lượng nước đất khai thác Nhận thức vấn đề này, điều quan trọng phải hiểu hệ thống nước đất khu vực lưu vực sông Rio Cobre phản ứng định lượng với thay đổi khí hậu hoạt động người để áp dụng quản lý thích ứng nhằm giảm thiểu vấn đề cải thiện điều kiện sống người cư trú khu vực dễ bị ảnh hưởng vấn đề liên quan tài nguyên nước đất Jamaica Nghiên cứu nhằm mục đích đánh giá lượng nước đất, ảnh hưởng khai thác khí hậu cách xây dựng mơ hình dịng chảy nước đất Mơ hình mơ diễn biến lịch sử tài nguyên nước đất theo tăng giảm khai thác Giới thiệu chung khu vực nghiên cứu Lưu vực Rio Cobre, nằm Jamaica có diện tích cắt ngang qua năm giáo xứ chủ yếu giáo xứ St Catherine Lưu vực bao gồm tiểu lưu vực Rio Cobre (phía Bắc); Một vùng nơng thơn đến vùng nội địa bán cơng nghiệp phần phía Nam (Hạ Rio Cobre); Đồng ven biển (Hình 1) Hình 1: Sơ đồ khu vực nghiên cứu 168 Hội thảo Khoa học Quốc gia 2023 Sông Rio Cobre kênh dẫn cho nước mặt nước đất từ lưu vực xuống lưu vực lưu vực Rio Cobre Ở phần lưu vực, hệ thống nước đất lưu vực biệt lập, khơng có dịng chảy vào khỏi tiểu lưu vực Nước đất thoát từ tầng chứa nước đá vơi phía kết thúc tạo thành nhiều suối đổ sông Rio Cobre Tại Bog Walk, nơi phân chia hai tiểu lưu vực, nước mặt từ thượng nguồn sông Rio Cobre chảy vào tiểu lưu vực phía Ở độ cao khoảng 76 m so với mực nước biển trung bình, Rio Cobre (sơng) chảy từ tiểu lưu vực thượng nguồn Rio Cobre Tài nguyên nước lưu vực sông Rio Cobre Jamaica bị tình trạng căng thẳng nghiêm trọng nhu cầu vượt khả cung ứng 94 % [4] Nước mặt nước đất đóng vai trị hai nguồn nước uống khu vực nghiên cứu theo Ủy ban Tài nguyên nước Quốc gia (NWC) Các giếng NWC khoan vào tầng ngậm nước cung cấp nước đất cho khu vực 90 % nhu cầu nước đảo đáp ứng nước đất, chủ yếu lưu trữ tầng ngậm nước đá vôi Việc khai thác nước đất thực phạm vi hai tầng chứa nước sau: Tầng chứa nước đá vôi thuộc bồn trũng Rio Cobre: Từ năm 1930, tầng chứa nước đá vôi Thượng Rio Cobre thức khai thác nguồn cung cấp nước đất cho hoạt động trồng mía Trong thời kỳ hạn hán 1975-1976, 291.000 m3 nước đất khai thác ngày năm 1976 dẫn đến xâm nhập mặn vào sâu đất liền với chiều dài 1,5 km, sản lượng khai thác lưu vực đánh giá 278.000 m3/ngày Ước tính thiết lập khơng xét đến áp thấp bơm thấp mực nước biển 37 giếng cấp phép quan Tài nguyên nước Jamaica (WRA) với công suất khai thác 120.893 m3/ngày, chiếm 43 % sản lượng khai thác tầng chứa nước Ngồi ra, 62 giếng có cơng suất bơm 396.871 m3/ngày phép khai thác tùy theo hoạt động bơm công suất cấp phép trước Tầng chứa nước phù sa Rio Cobre: Nước đất tầng chứa nước phù sa phía Đơng Nam Giáo xứ Catherine sử dụng để liên kết với hệ thống kênh tưới tiêu cho việc canh tác khu vực tưới chuối Bernard Lodge, nơi trồng vườn mía Tầng chứa nước có sản lượng 43.708 m3/ngày Cơ quan quản lý Tài nguyên nước cấp phép cho 21 giếng để khai thác 48.101 m3/ngày tầng chứa nước phù sa Rio Cobre (Đông) chiếm 4.993 m3/ngày (12 %) tổng sản lượng khai thác tầng chứa nước Một số giếng bơm có cơng suất cho phép 19.100 m3/ngày đặt phía Tây Bernard Lodge 52 giếng bơm khác với tốc độ bơm 126.026 m3/ngày cấp phép để đáp ứng hoạt động khai thác nước trước đây, WRA không chịu trách nhiệm cấp phép cho giếng bơm Phương pháp nghiên cứu liệu sử dụng Phần mềm GMS giao diện kết nối mơ hình nước đất MODFLOW, MODPATH, MT3D, PEST Nó cung cấp khả phân tích liệu thu chỗ, xây dựng mơ hình, hiệu chỉnh, đánh giá trình bày trực quan cho giai đoạn mô nước đất Cả hai mô hình sai phân hữu hạn phần tử hữu hạn, chẳng hạn MODFLOW 2000, ART3D, UTCHEM, FEMWATER SEEP2D GMS hỗ trợ dạng 2D 3D Cấu trúc module phần mềm cho phép người lập mơ hình chọn tham gia module tích hợp cụ thể, cho phép họ chọn kỹ lập mơ hình nước đất cần thiết Hội thảo Khoa học Quốc gia 2023 169 3.1 Xây dựng mơ hình dịng chảy nước đất cho khu vực nghiên cứu Mơ hình khái niệm thể thành phần thủy văn, vị trí địa lí lớp chứa nước, điều kiện biên xây dựng nghiên cứu trước “Đánh giá tài nguyên nước đất lưu vực sông Rio Cobre, Jamaica - Phương pháp tiếp cận phân tích hệ thống dịng chảy nước đất khu vực” [8] Ranh giới lưu vực xác định dựa phân chia nước đất thay phân chia thành phần thủy văn có điều chỉnh nhỏ diện tích nước (+113,3 km2) thực Q trình phân định khai thác độ cao bề mặt thực phần mềm QGIS phiên 3.16.1 - Hannover với hệ tọa độ Jamaica JAD - 2001 (Hình 2) Hình 2: Ranh giới khu vực nghiên cứu Lưới phân sai Lưu vực rời rạc hóa thành lưới tính tốn để tích phân hệ phương trình áp dụng mơ hình Modflow nhằm mơ tả q trình động thái nước đất Từ điều kiện số liệu địa hình tầng chứa nước, khu vực nghiên cứu chia thành mạng lưới (cells) với kích thước km × km gồm 72 hàng 120 cột với ô không hoạt động (màu trắng) Mơ hình hóa với mạng lưới với kích thước 200 × 200 m, bao gồm 96 hàng × 100 cột Độ phân giải thể tốt hệ thống dịng chảy khu vực với mơ dịng nước đầy đủ Có 39.188 nút giao 28.800 vng với 11.377 ô hoạt động 17.423 ô không hoạt động miền mơ hình (Hình 3) Hình 3: Lưới phân sai lưu vực nghiên cứu 170 Hội thảo Khoa học Quốc gia 2023 Các lớp mơ hình Mơ hình xây dựng gồm tầng dựa độ cao, đặc điểm địa chất thủy văn thạch học Mỗi lớp có độ cao mặt - mặt cụ thể lớp mơ hình thứ lớp bao gồm hệ thống chứa nước cao nguyên miền núi, lớp có độ cao dao động từ 446-976 m so với mực nước biển trung bình (a.m.s.l) Độ cao lấy từ đồ độ cao kỹ thuật số QGIS xử lý, độ cao giả định 100 m a.m.s.l Tầng chứa nước tạo thành từ đá vơi, phù sa bazan hình thành Lớp thứ hai thành tạo chủ yếu đá vôi đá bazan có độ cao từ 16-466 m a.m.s.l, cao độ đáy thấp m a.m.s.l Lớp lớp phủ thấp hệ thống tầng chứa nước xem xét mơ hình Nó bị chi phối hình thành đá vơi phù sa Vùng phù sa đại diện cho đồng St Catherine đóng vai trị quan trọng việc chi phối dòng nước đất Cao độ tầng dao động từ 0-109 m a.m.s.l cao độ mặt -200 m a.m.s.l cao độ mặt Điều kiện biên điều kiện ban đầu Như trình bày Hình 4, có hai điều kiện biên đặt cho khu vực nghiên cứu ranh giới khơng có dịng chảy bên phía Bắc (the lateral no-flow boundary), phía Tây, phía Đơng ranh giới kiểm sốt dịng chảy bên ngồi phía Nam (the external head-controlled boundary) lưu vực, biển Caribe Tầng chứa nước phù sa lớp mơ hình thứ ba hoạt động dịng chảy chắn, định nghĩa mơ hình khái niệm cách sử dụng HFB (Horizontal Flow Barrier) Hình 4: Điều kiện biên ranh giới dòng chảy ngang Điều kiện ban đầu Các mực nước đất không đổi mô từ nghiên cứu “Đánh giá tài nguyên nước đất lưu vực sông Rio Cobre, Jamaica - Phương pháp tiếp cận phân tích hệ thống dịng chảy nước đất khu vực” sử dụng làm mực nước ban đầu mơ hình tạm thời [8] Lượng bổ cập Lượng bổ cập nước mưa lưu vực xác định cách nhân lượng mưa hàng ngày khu vực với hệ số thấm nước mưa [8, 9] Biên sông (hệ thống sông Rio Cobre) Sông Rio Cobre nằm tầng chứa nước thứ hai ba Các thuộc tính sơng bao gồm độ cao đáy độ dẫn thủy lực, hạn chế liệu thu thập nên giả định đưa cột áp độ cao, độ dẫn thủy lực xem xét đặc điểm địa lý thổ nhưỡng khu vực sông Hội thảo Khoa học Quốc gia 2023 171 Điểm lộ Điểm lộ sử dụng để đại diện cho dòng chảy suối lưu vực MODFLOW-2000 [10] tính tốn lưu lượng từ tầng chứa nước cách lấy hiệu số mực nước đất cao trình điểm lộ Các điểm lộ mơ độ cao độ dẫn cho ô vị trí điểm lộ Hệ số nhả nước Một thông số địa chất thủy văn định nghĩa tầng chứa nước hệ số nhả nước đàn hồi (Specific Storage - Ss) hệ số nhả nước trọng lực (Specific Yield - Sy) Ước tính hệ số lưu trữ ban đầu xác định dựa nghiên cứu trước [11] Hệ số thấm Hệ số thấm (K) tính tốn dựa hình thành địa chất ba lớp chứa nước (m/d) 3.2 Dữ liệu chạy mơ hình dịng chảy khơng ổn định Từ kết thu từ mơ hình dịng chảy ổn định, mơ hình dịng chảy khơng ổn định xây dựng với dự liệu đầu vào thay đổi theo thời gian Mơ hình chạy hiệu chỉnh dựa liệu thu thập từ lỗ khoan quan trắc tương ứng với bước thời gian (time step) bước thời đoạn (stress period) thiết lập, cho kết đường đẳng mực nước biến đổi theo bước thời gian mô Bước thời đoạn (stress period) bước thời gian (time step) Liên quan đến việc thu thập liệu, liệu thay đổi theo thời gian bao gồm lượng mưa, mực nước đất quan trắc lưu lượng sơng liệu có sẵn giai đoạn 2000- 2020 với tần xuất đo đạc hàng tháng Do đó, bước thời đoạn với tần xuất hàng tháng chọn để mô Trong 20 năm, 252 bước thời đoạn chạy với bước thời gian, bước ngày 01 tháng 01 năm 2000 ngày 01 tháng 12 năm 2020 Giếng khoan khai thác nước 110 giếng bơm mô tầng chứa nước thứ có giếng bơm, tầng chứa nước thứ hai có 23 giếng tầng chứa nước thứ ba có 84 giếng bơm Các giếng khai thác tự phát bơm bất hợp pháp coi yếu tố khơng chắn mơ hình Hình 5: Vị trí phân bổ giếng khoan quan trắc tầng chứa nước Kết tính tốn Đường đẳng mực nước đất mơ từ mơ hình dịng chảy phụ thuộc vào thời gian trình bày so sánh với mực nước giếng quan trắc thể kết 172 Hội thảo Khoa học Quốc gia 2023 trình hiệu chỉnh đồng thời kết khác mơ hình bao gồm phương trình cân nước, dịng chảy tích lũy thay đổi lượng nước đất lưu trữ tầng chứa nước phân tích nhằm đưa hiểu biết cụ thể hệ thống nước đất lưu vực sông Rio Cobre Jamaica 4.1 So sánh mực nước đất quan trắc mô Hình biểu thị cột nước đất quan trắc tính tốn tầng chứa nước lưu vực sông Rio Cobre Hệ số nhả nước trọng lực đàn hồi điều chỉnh để giảm thiểu khác biệt cột nước đất quan trắc cột nước đất mơ từ mơ hình Những tương đồng hình thái xu chứng tỏ độ xác mơ hình hiệu chuẩn Khoảng 5-15 m khác biệt cột nước đất quan sát mô Các cột áp nước đất từ ba giếng quan trắc có vị trí dọc theo sơng Rio Cobre lưu vực phía cho thấy tương quan cột áp nước đất mô quan trắc Ngoài dao động theo mùa (mùa mưa mùa khô) tăng giảm theo năm tương ứng với lượng mưa, biểu đồ thể cột nước có áp khơng cho thấy suy giảm đáng kể suốt 20 năm quan trắc mơ Hình 6: So sánh cốt cao mực nước đất quan trắc thực tế tính tốn giếng quan trắc 4.2 Dòng chảy sở (measured baseflow) dòng cấp nước đất (computed baseflow) Hình 7: So sánh dòng chảy sở dòng cấp nước đất hai trạm Bog Walk Spanish Town Hội thảo Khoa học Quốc gia 2023 173 Hình mơ tả so sánh dịng chảy mùa kiệt sơng dòng chảy sở nước đất cung cấp hai trạm quan trắc Bog Walk (New) Spanish Town Trạm Bog Walk (New) nằm tầng chứa nước thứ hai nơi mà ba nhánh thượng nguồn sông Rio Cobre gặp Trạm Spanish Town thuộc tầng chứa nước thứ nằm nơi thượng lưu sông Rio Cobre đổ vào lưu vực hạ lưu Mặc dù thực tế dòng chảy kiệt quan sát hai trạm hiển thị đỉnh cực trị khác suốt thời gian mô phỏng, chúng tuân theo xu hướng thay đổi định theo mùa theo năm 4.3 Tính tốn cân nước Các thành phần dòng chảy vào mơ hình hiệu chỉnh trình bày Hình 8, năm khô hạn 2013-2016, cho thấy suy giảm lượng mưa bổ sung bắt đầu dòng chảy từ biển vào Tuy nhiên, mức độ xâm nhập không đáng kể Biểu đồ biểu thị thành phần dòng chảy cho thấy sau đạt đỉnh vào năm nhiều mưa 2002 2005, chúng có xu hướng suy giảm cuối đạt mức thấp vào năm khô hạn 2013 đến 2016 Hình 8: Thành phần dịng chảy vào dịng chảy mơ từ mơ hình 4.4 Tích lũy thay đổi sản lượng lưu trữ Hình cho thấy lưu lượng tích chứa tích luỹ tầng chứa nước liên tục giảm đạt 2,88×109 m3 vào tháng 12/2020 Hình 9: Sự thay đổi lưu lượng tích chứa tích lũy tầng chứa nước 174 Hội thảo Khoa học Quốc gia 2023 Sự suy giảm mực nước đất Hình 10 minh họa chênh lệch mực nước đất ngày 01/01/2000 ngày 01/12/2020 Mực nước sụt giảm lớn hình thành khu vực có địa hình cao phía Tây Bắc Đông Bắc lưu vực 40 m nước sụt giảm mơ phía Đơng tầng chứa nước Ở khu vực nơi hoạt động trồng trọt bơm thâm canh diễn ra, m suy giảm cột nước đất mơ Hình 10: Sự suy giảm mực nước đất tầng chứa nước Kết luận Mơ hình dịng chảy nước đất biến đổi theo thời gian (2000-2020) xây dựng hiệu chỉnh với mực nước đất đo lưu vực Rio Cobre Ở lưu vực Thượng Rio Cobre, tầng chứa nước thoát lượng lớn nước đất sông, suối chiếm 33,7 % tổng lượng nước tự nhiên Mực nước đất thập kỷ qua biến động theo lượng bổ cập tự nhiên nước đất, nhiên khơng có xu hướng sụt giảm liên tục quan sát Tính tốn cân nước cho thấy lượng cung cấp thấm từ mưa vào tầng chứa nước nguồn bổ cập chính, sau lượng bổ cập từ nước sơng Tầng chứa nước lượng nhỏ sông suối biển với tỷ lệ nhỏ tổng lượng nước thoát (~26 %) Trữ lượng nước đất giảm chủ yếu năm khô hạn mưa dẫn tới lượng bổ cập giảm Tuy nghiên cứu cho thấy lưu vực không trạng thái suy thoái nguồn nước đất, bên liên quan quan quản lý nước Jamaica cần có hành động thích hợp để bảo tồn bền vững cho tầng chứa nước lưu vực sông Rio Cobre Bằng cách tích hợp mơ hình nước đất mô tạm thời với chiến lược quản lý tài nguyên nước bền vững khác, hệ thống nước đất Rio Cobre bảo đảm khỏi suy giảm mức Kiến nghị Mơ hình có số nhược điểm cần giải để cải thiện kết mô số thiếu liệu đầu vào hiệu chuẩn Trong mơ dịng chảy nước đất phụ thuộc thời gian, lưu lượng bơm từ giếng khoan nên thay đổi từ lưu lượng không thay đổi theo thời gian thành lưu lượng thay đổi theo thời gian để phản ánh điều kiện tự nhiên tốt Dữ liệu khai thác nên thu thập chi tiết bao gồm giếng cấp phép không cấp phép Mực nước đất quan sát không đầy đủ liệu, phương pháp nội suy ngoại suy cần thực để lấp đầy thiếu sót cải thiện mức độ xác mơ hình Sự chênh lệch cao độ mực nước đất số liệu quan trắc mô tầng chứa nước thứ ba đặc biệt khu vực có tầng địa chất phù sa lớn, cần tiếp tục tiến hành Hội thảo Khoa học Quốc gia 2023 175 nghiên cứu bổ sung liên quan đến đặc điểm thạch học địa chất thủy văn Các trạm khí tượng thủy văn cung cấp lượng mưa dòng chảy cần cải thiện hệ thống quản lý liệu, phần mềm hỗ trợ kỹ thuật khác cần cài đặt để theo dõi, lưu trữ xử lý liệu hiệu TÀI LIỆU THAM KHẢO [1] Taylor R., Tindimugaya C (2012) The impacts of climate change and rapid development on weathered crystalline rock aquifer systems in the humid tropics of sub-Saharan Africa: Evidence from south- western Uganda Climate Change Effects on Groundwater resources: A Global Synthesis of Findings and Recommendations, 17-32 [2] Houghton J.T., Ding Y., Griggs D.J., Noguer M., Van der Linden P.J., Dai X., Maskell K., Jonson C (2001) Contribution of Working group I to the third assessment report of the Intergovernmental panel on Climate Change Cambridge university press Climate Change [3] Peterson R.N., Burnett W.C., Glenn C.R., Johnson A.G (2009) Quantification of point - source groundwater discharges to the ocean from the shoreline of the Big Island, Hawaii Limnology and Oceanography, Vol 54, 890-904 [4] Water Resources Authority (2005) Water resources development master plan [5] Yihdego Y., Reta G., Becht R (2016) Human impact assessment through a transient numerical modelling on the UNESCO World Heritage Site, Lake Naivasha, Kenya Environmental Earth Sciences, Vvol 79 [6] Wang S., Shao J., Song X., Zhang Y., Huo Z., Zhou X (2008) Application of MODFLOW and geographic information system to groundwater flow simulation in North China Plain, China Environmental geology, Vol 55 [7] Aquaveo (2021) MODFLOW modelling with GMS [8] Henry Uton (2021) Groundwater resources Assessment of the Rio Cobre basin, Jamaica MSc thesis [9] Botbol (1982) Lower Rio Cobre limestone aquifer Kingston, Jamaica Hydrogeology [10] Harbaugh A.W., Banta E.R., Hill M.C., McDonald M.G (2000) Geological survey modular ground- water model-user giude to modularization concepts and the groundwater flow process [11] Kuang X., Jiao J.J., Zheng C., Cherry J.A., Li H (2020) A review of specific storage in aquifers Journal of Hydrology, Vol 581 BBT nhận bài: 27/7/2023; Chấp nhận đăng: 15/9/2023 176 Hội thảo Khoa học Quốc gia 2023