Bài viết Kết hợp dữ liệu vệ tinh trọng lực và thủy văn bề mặt để theo dõi diễn biến trữ lượng nước ngầm tại khu vực Bắc Trung Bộ phân tích số liệu cho thấy giải pháp thu được thông qua việc loại bỏ độ ẩm đất (SM) trong mô hình GLDAS là có đủ cơ sở.
KẾT HỢP DỮ LIỆU VỆ TINH TRỌNG LỰC VÀ THỦY VĂN BỀ MẶT ĐỂ THEO DÕI DIỄN BIẾN TRỮ LƯỢNG NƯỚC NGẦM TẠI KHU VỰC BẮC TRUNG BỘ Lê Tiến Duy, Trần Thành Lê, Hoàng Văn Tuấn, Lê Thị Liên, Lê Thị Hoa Huệ Trường Đại học Tài nguyên Mơi trường Hà Nội Tóm tắt Nhiệm vụ vệ tinh Phục hồi trọng lực thí nghiệm khí hậu (GRACE) áp dụng để xác định biến động trữ lượng nước (Terrestrial Water Storage - TWS), chủ yếu bao gồm nước ngầm (Ground Water - GW), nước mặt (Surface Water - SW), tức độ ẩm đất (Soil Moisture - SM) nước tuyết tương đương (Snow Water Equivalent - SWE) Để theo dõi biến động trữ lượng nước ngầm, nhóm tác giả xác định biến động trữ lượng nước (TWS) dựa vào liệu vệ tinh trọng lực GRACE Sau đó, tiến hành loại bỏ yếu tố nước mặt mơ theo mơ hình bề mặt đất Hệ thống đồng hóa liệu đất tồn cầu GLDAS (Global Land Data Assimilation System) Biến động trữ lượng nước (TWS) xác định từ sản phẩm GRACE (CSR, GFZ JPL) mơ hình GLDAS Phân tích số liệu cho thấy giải pháp thu thơng qua việc loại bỏ độ ẩm đất (SM) mô hình GLDAS có đủ sở Dựa quy trình xử lí liệu vệ tinh GRACE WeiFeng, sau đó, loại bỏ biến động liệu nước mặt từ NAOH mơ hình GLDAS, chúng tơi tiến hành ước tính biến động nước ngầm khu vực nghiên cứu, đồng thời, thu thập số liệu nước ngầm quan trắc chỗ tầng chứa nước không áp 03 tỉnh: Thanh Hóa, Nghệ An Hà Tĩnh, đại diện cho khu vực nghiên cứu để theo dõi xu hướng tăng giảm mực nước ngầm, nhằm đưa kết luận dự báo cho toàn vùng Dựa xu hướng biến động trữ lượng nước hàng năm nước mặt CSR - Nước mặt, GFZ - Nước mặt, JPL - Nước mặt cho thấy, mơ hình CSR -Nước mặt phù hợp khu vực nghiên cứu Trên sở đó, kết luận biến động trữ lượng nước ngầm khu vực Bắc Trung Bộ có xu hướng giảm nhẹ: - 0,09 ± 0,25 cm, giai đoạn từ năm 2011 - 2017 Từ khóa: Vệ tinh Phục hồi trọng lực Thí nghiệm khí hậu; Hệ thống đồng hóa liệu đất tồn cầu; Tổng trữ lượng nước; Nước ngầm; Nước mặt; Độ ẩm đất Abstract Combining gravity satellite data and surface hydrology to monitor changes in groundwater reserves in the North central region The Gravity Recovery and Climate Experiment (GRACE) satellite mission is applied to determine the variation of water storage (TWS), mainly consisting of groundwater (GW), surface water (SW), ie soil moisture (SM) and snow water equivalent (SWE) To monitor groundwater reserve fluctuations, we determine changes in terrestrial water storage (TWS) based on GRACE gravity satellite data, then proceed to remove surface water simulated by the Global Land Data Assimilation System model (GLDAS) Changes in terrestrial water storage were determined from GRACE products (CSR, GFZ, JPL) and GLDAS model, data analysis showed that the solution was obtained through soil moisture removal (SM) in the GLDAS model is enough basis Based on WeiFeng’s GRACE satellite data processing process, then removing surface water data fluctuations from NAOH in the GLDAS model, we estimate groundwater fluctuations in the study area, and at the same time collect Collecting groundwater data for in situ monitoring in the unconfined aquifer in provinces of Thanh Hoa, Nghe An, Ha Tinh representing the study area to monitor the trend of increasing and decreasing groundwater level, for the purpose of providing forecast conclusions for the whole region Based on the annual trend of water reserve variation of Surface Water and 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 461 CSR - Surface Water, GFZ - Surface Water, JPL - Surface Water, the CSR - Surface Water model is the most suitable in the study area On that basis, it can be concluded that the fluctuation of groundwater reserves in the North Central region trends to decrease slightly: - 0.09 ± 0.25cm in the period from 2011 to 2017 Keywords: Gravity Recovery Satellite and Climate Experiment; Global Land Data Assimilation System; Total water reserves; Groundwater; Surface water; Soil moisture Đặt vấn đề Biến đổi khí hậu nguyên nhân khiến cho mối quan hệ yếu tố cân vốn có Điều mối đe dọa an ninh lương thực, chất lượng nguồn nước nguồn lượng Biến đổi khí hậu tượng mà người phủ nhận tác động hiển nhiên [1] Tài nguyên nước (TNN) là nguồn nước mà người sử dụng sử dụng vào mục đích khác nhau, nước dùng hoạt động nông nghiệp, công nghiệp, dân dụng, giải trí và mơi trường, hầu hết hoạt động cần nước Vấn đề lớn với việc quản lý tài nguyên nước theo dõi, kiểm soát nguồn nước định kì Sự sụt giảm nguồn nước ngầm vơ hình hậu thường thể vào nhiều năm sau Dễ thấy rằng, mơ hình hạn hán gắn chặt với cạn kiện nước ngầm Nước bề mặt phục hồi sau có nhiều mưa lớn nhiều để phục hồi lượng nước bị hút khỏi tầng ngậm nước Sự suy giảm tài nguyên nước đất khu vực Bắc Trung Bộ vấn đề ban ngành quan tâm, khu vực có hai tầng chứa nước khai thác sử dụng, tầng chứa nước Holocen (qh) Pleistocen (qp) Xuất phát từ tầm quan trọng việc dự báo TNN, lựa chọn phương pháp sử dụng liệu vệ tinh Phục hồi trọng lực Thí nghiệm khí hậu (GRACE) Hệ thống đồng hóa liệu đất tồn cầu GLDAS để theo dõi biến động trữ lượng nước ngầm khu vực Bắc Trung Bộ Sứ mệnh GRACE bao gồm 02 vệ tinh thiết kế để lập đồ chuẩn xác biến đổi trọng trường Trái đất Hai vệ tinh giống hệt bay quỹ đạo cực, độ cao khoảng 500 km, phía Trái đất, đặt cách khoảng cách 220 km quỹ đạo Khoảng cách 02 vệ tinh đo liên tục GPS radar Dữ liệu GRACE sử dụng nghiên cứu bao gồm 60 dị thường TWS hàng tháng loại RL05, phát hành từ tháng 03 năm 2011 đến tháng 06 năm 2017, Trung tâm Nghiên cứu vũ trụ (Center for Space Research - CSR), Trung tâm Nghiên cứu Khoa học Địa chất Đức (German Research Center for Geosciences - GFZ) từ Phịng thí nghiệm Động phản lực NASA (Jet Propulsion Laboratory - JPL) cung cấp Theo số nghiên cứu gần chứng minh rằng, mơ hình GLDAS phù hợp mơ hình thủy văn khác ước lượng GRACE-GLDAS tương ứng tốt với ước lượng từ việc theo dõi giếng quan trắc [2, 3, 4] Vì vậy, nghiên cứu này, sử dụng mô hình GLDAS kết hợp với liệu vệ tinh 03 tổ chức CSR, GFZ JPL để theo dõi diễn biến trữ lượng nước ngầm cho tỉnh khu vực Bắc Trung Bộ Cơ sở lí thuyết phương pháp nghiên cứu 2.1 Giới thiệu phạm vi nghiên cứu Bắc Trung Bộ là phần phía Bắc của Trung Bộ, kéo dài từ tỉnh Thanh Hóa đến phía Bắc của đèo Hải Vân, vùng kinh tế trọng điểm miền Trung Với vị trí địa lý sau: phía Bắc giáp 462 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 Trung du miền núi Bắc Bộ, Đồng Sơng Hồng; phía Nam giáp Duyên hải Nam Trung Bộ; phía Tây giáp dãy Trường Sơn Lào; phía Đơng Biển Đơng (Vịnh Bắc Bộ) Khu vực bao gồm Trung du miền núi, địa hình phân dị phức tạp, thời tiết khắc nghiệt, nhiều biến động Phạm vi nghiên cứu bao gồm 06 tỉnh: Thanh Hóa, Nghệ An, Hà Tĩnh, Quảng Bình, Quảng Trị, Thừa Thiên Huế, giới hạn từ 16 đến 20.75 vĩ độ Bắc từ 103.75 đến 108.25 kinh độ Đơng (Hình 1) Hình 1: Sơ đồ khu vực Bắc Trung Bộ 2.2 Cơ sở lí thuyết phương pháp nghiên cứu Vệ tinh Phục hồi Trọng lực Thí nghiệm Khí hậu (GRACE) áp dụng để theo dõi diễn biến thay đổi tổng lượng nước (TWS), chủ yếu bao gồm nước ngầm, độ ẩm đất tuyết [5] ∆TWS = ∆GW + ∆SM + ∆SWE (1) Khu vực nghiên cứu khơng có tuyết nên phương trình (1) viết lại: ∆TWS = ∆GW + ∆SM (2) Trong đó: ∆GW thay đổi nước ngầm; ∆TWS thay đổi tổng lượng nước; ∆SM thay đổi độ ẩm đất Trong đó, biến thể từ SM mơ từ mơ hình thủy văn bề mặt đất Do đó, ước tính thay đổi mực nước ngầm cách loại bỏ biến thể nước mặt khỏi tổng lượng nước có nguồn gốc từ GRACE 2.2.1 Thu thập xử lý liệu GRACE Dữ liệu GRACE sử dụng để ước tính TWS nghiên cứu bao gồm 60 trường trọng lực Release (RL05) hàng tháng, từ tháng năm 2011 đến tháng năm 2017, cung cấp Trung tâm Nghiên cứu Không gian (CSR), Trung tâm Nghiên cứu Khoa học Địa chất Đức (GFZ) từ Phòng thí nghiệm Sức đẩy Phản lực NASA (JPL) Kết TWS có 16 tháng bị thiếu liệu quản lý pin (tháng 06/2011; 05/2012; 10/2012; 03/2013; 08/2013; 09/2013; 02/2014; 07/2014; 12/2014; 06/2015; 10/2015; 11/2015; 04/2016; 09/2016; 10/2016; 02/2017) xác định phép nội suy tuyến tính 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 463 Dựa nghiên cứu Cheng Tapley, đề xuất thay hệ số C20 ước tính từ giải pháp định vị Laser vệ tinh (SLR), đó, giá trị thu từ quan sát GRACE có độ không đảm bảo so với giá trị SLR [6], thay đổi theo mùa sóng hài hình cầu bậc đại diện cho biến thể Geocenter Earth, không cung cấp riêng GRACE Sử dụng lọc Gaussian với bán kính làm mịn 300 km lọc tương quan “P2M5” (đối với hệ số hài hình cầu bậc trở lên, đa thức bậc lắp vào cặp chẵn sau đa thức phù hợp loại bỏ khỏi hệ số, tương tự áp dụng với cặp lẻ) để giảm thiểu ảnh hưởng lỗi “sọc” [7, 8, 9] Sau xử lý có liên quan, mơ hình tính tốn lưu trữ nước cạn (EWH) tính theo [10, 11] (3) Trong đó: : Khối lượng riêng trung bình Trái đất (5517 kg/m3); : Khối lượng riêng nước (giả định 1000 kg / m3); : Diện tích góc vùng; kl: Hệ số love number bậc l, đại diện cho phản ứng Trái đất tải trọng bề mặt thu [12]; : Là hệ số hài hình cầu mơ tả hình dạng khu vực; W tương ứng với toán tử làm trơn Gaussian; : Là phần dư hệ số điều hịa hình cầu trường trọng lực, đó, giá trị trung bình dài hạn bị loại bỏ Tồn q trình xử lý liệu GRACE áp dụng theo quy trình FengWei chạy phần mềm Matlab 2.2.2 Dữ liệu thủy văn bề mặt Được phát triển chung Trung tâm bay vũ trụ Goddard NASA (GSFC) Trung tâm Dự báo Môi trường Đại dương Quốc gia (NOAA) (NCEP), GLDAS trình bày biến thể địa hình gần thời gian thực tối ưu hóa cung cấp thơng tin lượng độ ẩm đất bề mặt toàn cầu Sự phát triển mơ hình thủy văn bề mặt đất giúp thu phân bố theo thời gian không gian biến thể độ ẩm quy mô lớn Trong nghiên cứu này, sử dụng độ ẩm đất trung bình tương đương từ mơ hình GLDAS toàn cầu, sau lọc liệu độ ẩm đất phần mềm Panoply, thu liệu nước mặt tháng 06 tỉnh Bắc Trung Bộ Kết thảo luận 3.1 Các biến thể lưu trữ nước cạn Vì biến thể tích trữ nước có tín hiệu theo mùa nên phân tích hồi quy đa tuyến tính áp dụng để kiểm tra biến thiên theo thời gian đại lượng thủy văn, chẳng hạn TWS ước tính Do đó, chuỗi thời gian định, mơ hình sử dụng cơng việc tính đến [16]: 464 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 (4) Trong đó: t thời gian; a số; , liên quan đến biên độ, pha tần số; k đại diện cho bậc sóng hài hòa (k = k = tương ứng với thành phần hàng năm nửa năm tương ứng); thay đổi lại liệu, chủ yếu nhiễu với số tín hiệu dư Quan sát Hình Bảng 1, mơ hình từ GRACE GLDAS có biên độ rõ rệt, theo chu kì hàng năm, trữ lượng nước thường tăng cao vào tháng 8, 9, 10 giảm thấp vào tháng 02, 3, Duy có mơ hình GLDAS thể dị thường khác biệt vào thời điểm cuối năm 2016, thời điểm lũ lụt diễn cực mạnh tỉnh Bắc Trung Bộ Hình 2: Biến động trữ nước có nguồn gốc từ GRACE (từ tổ chức CSR, GFZ JPL) từ mơ hình GLDAS khu vực Bắc Trung Bộ Bảng Biên độ xu hướng hàng năm trữ lượng nước ước tính từ liệu vệ tinh trọng lực GRACE Nước mặt từ mô hình GLDAS Biến số CSR GFZ JPL Nước mặt Biên độ hàng năm (cm) 7,42 ± 0,62 6,1 ± 0,71 5,91 ± 0,59 4,01 ± 0,52 Xu hướng (cm) -0,15 ± 0,24 -0,18 ± 0,28 -0,05 ± 0,23 -0,06 ± 0,2 RMSD (cm) 1,93 2,21 1,84 1,6 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 465 Hình 3: Biến động mực nước ngầm xu hướng hàng năm ước tín sau loại bỏ biến động Nước mặt từ liệu vệ tinh GRACE - CSR Hình 4: Biến động mực nước ngầm xu hướng hàng năm ước tính sau loại bỏ biến động Nước mặt từ liệu vệ tinh GRACE - GFZ Hình 5: Biến động mực nước ngầm xu hướng hàng năm ước tính sau loại bỏ biến động nước mặt từ liệu vệ tinh GRACE - JPL 466 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 Bảng Biên độ xu hướng biến động trữ lượng nước vệ tinh trọng lực GRACE nước mặt từ mơ hình thủy văn GLDAS - NAOH Biến số CSR - Nước mặt GFZ - Nước mặt JPL - Nước mặt Biên độ hàng năm (cm) 3,41 ± 0,64 2,09 ± 0,72 1,9 ± 0,64 Xu hướng (cm) -0,09 ± 0,25 -0,12 ± 0,28 0,02 ± 0,25 RMSD (cm) 1,97 2,25 2,0 Sau loại bỏ liệu nước mặt từ mơ hình GRACE, nhóm tác giả nhận thấy nước ngầm khu vực có xu hướng giảm, có số liệu tính theo JPL có dấu hiệu tăng nhẹ (Bảng 2) Quan sát Bảng Hình 3, 4, cho thấy, mơ hình CSR - Nước mặt có độ lệch chuẩn biên độ lỗi trung bình bình phương (RMSD) hay độ lệch chuẩn phần dư bé, điều chứng tỏ mơ hình CSR - SM phù hợp 03 mơ hình Tuy nhiên, để có số liệu kiểm chứng kết nghiên cứu, nhóm tác giả tiến hành thu thập số liệu quan trắc nước ngầm 03 tỉnh Thanh Hóa, Nghệ An Hà Tĩnh với thời gian quan trắc từ tháng năm 2011 đến tháng năm 2017 3.2 Phân tích so sánh biến thể nước ngầm Số liệu quan trắc nước ngầm 03 tỉnh Thanh Hóa, Nghệ An Hà Tĩnh, thời gian quan trắc từ tháng năm 2011 đến tháng năm 2017 (đây khoảng thời gian địa điểm thu thập liệu từ giếng quan trắc chỗ, địa phương khác hầu hết đến năm 2018 xây dựng giếng quan trắc) 21 20.5 20 Lat 19.5 19 18.5 18 17.5 103.5 104 104.5 105 Long 105.5 106 106.5 107 Tồn vị trí giếng xác định nằm 08 lưới, kích thước 0,25 × 0,25 độ, với tổng số 18 giếng quan trắc, số liệu thống kê theo độ cao mực nước trung bình ô lưới tầng chứa nước không áp (qh) (Hình 6) Hình 6: Sơ đồ khu vực 03 tỉnh Thanh Hóa, Nghệ An Hà Tĩnh (bao quanh đường cong màu xanh) phân bố giếng quan trắc nước ngầm (chấm đỏ) Các hình chữ nhật màu đen đại diện cho lưới có độ rộng kinh độ vĩ độ (0,250 × 0,250) Quan sát Hình nhóm tác giả nhận thấy, biến động nước ngầm hàng năm tầng chứa nước không áp khu vực 03 tỉnh Thanh Hóa, Nghệ An Hà Tĩnh có xu hướng giảm (- 0,53 ± 0,18 cm) Nếu nhìn vào kết quả, ước tính biến động vùng Bắc Trung Bộ là: - 0,09 ± 0,25 cm (Bảng Hình 3) lí giải tương đồng xu hướng giảm nhẹ hàng năm sở địa hình 03 tỉnh Thanh Hóa, Nghệ An Hà Tĩnh - đại diện cho vùng Bắc Trung Bộ 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 467 Hình 7: Biến động mực nước ngầm xu hướng hàng năm khu vực 03 tỉnh Thanh Hóa, Nghệ An Hà Tĩnh theo số liệu quan trắc chỗ giếng tầng chứa nước không áp (qh) Kết luận đề xuất Căn theo ước tính biến động trữ lượng nước bề mặt đất từ liệu vệ tinh trọng lực GRACE mơ hình thủy văn bề mặt GLDAS, hồn tồn có đủ sở để ước tính biến đổi trữ lượng nước ngầm phạm vi đủ lớn, tài ngun vơ q giá nằm lòng đất mà từ trước tới nhiều cơng sức, chi phí để theo dõi quản lí Căn số liệu Bảng 2, chúng tơi thấy mơ hình CSR phù hợp với khu vực nghiên cứu, sở kết luận nước ngầm khu vực Bắc Trung Bộ có xu hướng cạn kiệt nhẹ khoảng thời gian từ năm 2011 đến năm 2017 ước tính - 0,09 ± 0,25 cm Theo dõi biến động trữ lượng nước ngầm dựa cảm nhận từ xa thay đổi trữ lượng nước mặt đất, ước tính từ GRACE mơ hình GLDAS thực cần thiết Hiện nay, giới ứng dụng mạnh mẽ liệu vệ tinh trọng lực GRACE theo dõi quản lí tài nguyên nước, đặc biệt nước ngầm phạm vi lớn, Việt Nam có khó khăn định để ứng dụng loại liệu nghiên cứu khoa học như: địa hình vùng miền tương đối hẹp; giếng quan trắc nước ngầm chỗ gần tăng cường xây dựng bổ sung, nên, chưa có nhiều nghiên cứu ứng dụng lĩnh vực Khoa học Trái đất, đặc biệt lĩnh vực Tài nguyên Môi trường Với xu hướng phát triển lĩnh vực Khoa học Vũ trụ NASA, chắn sứ mệnh vệ tinh lực GRACE tương lai tăng cường độ không gian Vậy nên, mong muốn trọng thời gian tới có nhiều nghiên cứu ứng dụng vệ tinh trọng lực GRACE, GRACE - FO, nhằm phát triển lĩnh vực Khoa học Trái đất, phục vụ cho chuyên ngành đặc thù, đặc biệt ngành Tài nguyên Môi trường Hy vọng kết nghiên cứu góp phần nhỏ vào nghiên cứu sau thuộc lĩnh vực Tài nguyên Môi trường TÀI LIỆU THAM KHẢO [1] Bộ Tài nguyên Môi trường (2016) Kịch biến đổi khí hậu, nước biển dâng cho Việt Nam Nhà xuất Tài nguyên Môi trường Bản đồ Việt Nam [2] Jin, S.; Feng, G (2013) Large-scale global groundwater variations from satellite gravimetry and 468 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 hydrological models, 2002 - 2012 Glob Planet Chang, 106, 20 - 30 [3] Jin, S G.; Hassan, A A.; Feng, G P (2012) Assessment of terrestrial water contributions to polar motion from GRACE and hydrological models J Geodyn, 62, 40 - 48 [4] Rodell, M.; Velicogna, I.; Famiglietti, J S (2009) Satellite-based estimates of groundwater depletion in India Nature, 460, 999 - 1002 [5] Rodell, M.; Chen, J.; Kato, H.; Famiglietti, J S.; Nigro, J.; Wilson, C R (2007) Estimating groundwater storage changes in the Mississippi River basin (USA) using GRACE Hydrogeol J, 15, 159 - 166 [6] Cheng, M K.; Tapley, B D (2004) Variations in the Earth’s oblateness during the past 28 years J Geophys Res.: Solid Earth Doi:10.1029/2004JB003028 [7] Swenson, S.; Wahr, J (2006) Post-processing removal of correlated errors in GRACE data Geophys Res Lett, doi:10.1029/2005GL025285 [8] Landerer, F W.; Swenson, S C (2012) Accuracy of scaled GRACE terrestrial water storage estimates Water Resour Res, 48, W04531 [9] Duan, X J.; Guo, J Y.; Shum, C K.; Van der Wal, W (2009) On the postprocessing removal of correlated errors in GRACE temporal gravity field solutions J Geod, 83, 1095 - 1106 [10] Wahr, J.; Molenaar, M.; Bryan, F (1998) Time variability of the earth’s gravity field: Hydrological and oceanic effects and their possible detection using GRACE J Geophys Res.: Solid Earth, 103, 30205 - 30229 [11] Chao, B F (2005) On inversion for mass distribution from global (time-variable) gravity field J Geodyn, 39, 223 - 230 [12] Han, D.; Wahr, J (1995) The viscoelastic relaxation of a realistically stratified earth, and a further analysis of postglacial rebound Geophys J Int, 120, 287 - 311 [13] Fetter, C W (1994) Applied Hydrogeology 3rd ed.; Prentice Hall: Upper Saddle River, NJ, USA [14] Sun, A Y.; Green, R.; Rodell, M.; Swenson, S (2010) Inferring aquifer storage parameters using satellite and in situ measurements: Estimation under uncertainty Geophys Res Lett, doi:10.1029/2010GL043231 [15] Viện Khoa học vật liệu - Viện Khoa học Công nghệ Việt Nam (2012) Đặc điểm phân bố hình thành trữ lượng nước đất vùng cát ven biển Quảng Bình, Phan Văn Trường [16] Ruya Xiao, Xiufeng He, Yonglei Zhang, Vagner G Ferreira, Liang Chang (2015) Monitoring Groundwater Variations from Satellite Gravimetry and Hydrological Models: A Comparison with in-situ Measurements in the Mid-Atlantic Region of the United States Remote Sens Ngày chấp nhận đăng: 10/11/2021 Người phản biện: TS Trần Thành Lê 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 469 PHỤ LỤC: BIẾN ĐỘNG MỰC NƯỚC TRUNG TẠI KHU VỰC 03 TỈNH THANH HÓA, NGHỆ AN, HÀ TĨNH TỪ THÁNG 3/2011 ĐẾN THÁNG 6/2017 TT 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 470 Thời gian 2011.205479 2011.289041 2011.372603 2011.458904 2011.545205 2011.624658 2011.708219 2011.791781 2011.832877 2011.991781 2012.043836 2012.126027 2012.208219 2012.258904 2012.358904 2012.458904 2012.542466 2012.627397 2012.704110 2012.794521 2012.884932 2012.961644 2013.043836 2013.121918 2013.212329 2013.302740 2013.372603 2013.456164 2013.539726 2013.623744 2013.707763 2013.791781 2013.875342 2013.958904 2014.023288 2014.114384 2014.205479 h(m) -0.010183152 -0.01317033 -0.018845713 -0.028447796 -0.007703554 0.013259604 0.040959232 0.047203037 0.025358037 0.01693762 0.009669287 0.005313165 -0.001315094 -0.014954471 -0.029868908 -0.014912701 -0.021164268 -0.020851968 0.049131483 0.054084617 0.038753988 0.035006576 0.029991684 0.018505836 0.009419816 0.001101773 -0.005172095 -0.017151396 -0.001640319 0.033600247 0.057877604 0.062797874 0.048773826 0.033933845 0.017115437 0.013383479 0.008007358 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 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 2014.289041 2014.372603 2014.447945 2014.536301 2014.624658 2014.708219 2014.791781 2014.875342 2014.967808 2015.060274 2015.124658 2015.205479 2015.289041 2015.319178 2015.428082 2015.536986 2015.624658 2015.704110 2015.795434 2015.886758 2015.978082 2016.043836 2016.121918 2016.208219 2016.296575 2016.384932 2016.458904 2016.539726 2016.641096 2016.730137 2016.819178 2016.908219 2016.983562 2017.056164 2017.151370 2017.246575 2017.312329 2017.361644 2017.441096 0.0036434 -0.006668612 -0.006843608 -0.002923882 -0.006217015 0.011912713 0.021827315 0.020896089 0.002107624 -0.014553296 -0.01254988 -0.018503407 -0.016274148 -0.027989769 -0.043495692 -0.035935587 -0.03011881 -0.01179474 -0.000957425 0.005607912 0.001540409 -0.008534021 -0.007477267 -0.017214536 -0.023132369 -0.02417283 -0.037106682 -0.042343395 -0.036293867 0.0031664 0.027815574 0.017664978 -0.008694853 -0.013031247 -0.021343996 -0.023646256 -0.029193185 -0.027814635 -0.026161063 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 471 ... Tĩnh theo số liệu quan trắc chỗ giếng tầng chứa nước không áp (qh) Kết luận đề xuất Căn theo ước tính biến động trữ lượng nước bề mặt đất từ liệu vệ tinh trọng lực GRACE mơ hình thủy văn bề mặt. .. liệu vệ tinh 03 tổ chức CSR, GFZ JPL để theo dõi diễn biến trữ lượng nước ngầm cho tỉnh khu vực Bắc Trung Bộ Cơ sở lí thuyết phương pháp nghiên cứu 2.1 Giới thiệu phạm vi nghiên cứu Bắc Trung Bộ? ?là... phù hợp với khu vực nghiên cứu, sở kết luận nước ngầm khu vực Bắc Trung Bộ có xu hướng cạn kiệt nhẹ khoảng thời gian từ năm 2011 đến năm 2017 ước tính - 0,09 ± 0,25 cm Theo dõi biến động trữ lượng