Xây dựng bản đồ phát thải khí nhà kính trên đất phù sa trồng lúa vùng Bắc Trung Bộ

Nghiên cứu này trình bày kết quả xây dựng bản đồ phát thải khí nhà kính (KNK) thông qua kỹ thuật GIS từ kết quả mô phỏng phát thải khí nhà kính CH4 , N2 O (GHG) từ đất canh tác lúa trên đất phù sa của khu vực Bắc Trung Bộ theo mô hình tính toán (DNDC).

Tạp chí Khoa học Cơng nghệ Nơng nghiệp Việt Nam - Số 06(127)/2021 VN are de cient in some typical macro and medium-nutrients (K, S) and micronutrients (B, Cu) To balance for these nutrients de ciencies, it is necessary to apply 90 kg K + 30 kg S + 0.25 kg B + 2.5 kg Cu per To improve the physio-chemical properties of sandy soils, bentonite is one of the soil amendments helping to increase water holding capacity, increase CEC, leading to increase crop productivity Adoption of enhanced irrigation technologies (sprinkler + mini-pan for peanut; dropping + mini-pan for mango) saved 30 - 70% amount of water as compared to traditional irrigation methods of farmers; increased crop yield (peanut, mango) from 12 to 30%; increased economic e ciency from 20 to 70% Assessment of groundwater quality in some coastal areas (Ninh uan and Phu Yen provinces) showed that excessive uses of inorganic fertilizers, residues from livestock farms and ooding irrigation were the main reasons causing the pollution of groundwater Keywords: ACIAR, sandy soils, SCC VN, mini-pan, mango, peanut Ngày nhận bài: 05/12/2020 Ngày phản biện: 15/12/2020 Người phản biện: TS Ngô Đức Minh Ngày duyệt đăng: 22/12/2020 XÂY DỰNG BẢN ĐỒ PHÁT THẢI KHÍ NHÀ KÍNH TRÊN ĐẤT PHÙ SA TRỒNG LÚA VÙNG BẮC TRUNG BỘ Nguyễn Văn iết1, Bùi ị Phương Loan1 TÓM TẮT Nghiên cứu trình bày kết xây dựng đồ phát thải khí nhà kính (KNK) thơng qua kỹ thuật GIS từ kết mơ phát thải khí nhà kính CH4, N2O (GHG) từ đất canh tác lúa đất phù sa khu vực Bắc Trung Bộ theo mơ hình tính tốn (DNDC) Mơ hình hiệu chỉnh tham chiếu với số liệu quan trắc phát thải thực địa tỉnh Quảng Trị, Nghệ An ừa iên - Huế Nghiên cứu dự báo phát thải KNK từ đất lúa cho khu vực đến năm 2030 cách sử dụng liệu đất, khí hậu, sử dụng đất sở liệu quản lý canh tác lúa, kịch khí hậu nước biển dâng cho Viêt Nam Bộ Tài nguyên Môi trường năm 2016 Kết nghiên cứu cho thấy, lượng phát thải KNK vụ Mùa có xu hướng cao vụ Xuân tỉnh Quảng Bình, Quảng Trị ừa iên - Huế Hà Tĩnh, Nghệ An anh Hoá phát thải vụ Xuân lại cao vụ Mùa Dự báo đến năm 2030, ừa iên - Huế có phát thải KNK cao (tương đương 11,512 CO2/ha) Nghệ An phát thải KNK thấp Ngoài ra, phát thải KNK từ tỉnh khu vực Bắc Trung Bộ nước ta có xu hướng tăng vào năm 2030 Từ khóa: DNDC, đất phù sa trồng lúa, GHG I ĐẶT VẤN ĐỀ Tình trạng phát thải KNK từ hoạt động sản xuất nước nhiều năm qua gây hiệu ứng nhà kính phạm vi tồn cầu làm cho khí hậu trái đất thay đổi, nước biển dâng, thiên tai, dịch bệnh ngày nhiều (Nguyen Hong Son et al., 2014) Nông nghiệp Việt Nam ngày phải đối mặt với nhiều rủi ro từ biến đổi khí hậu (BĐKH), đặc biệt điều kiện Việt Nam nước chịu ảnh hưởng nặng nề BĐKH eo báo cáo kiểm kê KNK Bộ Tài ngun Mơi trường, nơng nghiệp Việt Nam đóng góp 43,1% vào tổng lượng phát thải KNK Việt Nam (Ministry of Natural Resources and Environment, 2010) Hiện nay, có nhiều mơ hình sử dụng để mơ tính tốn phát thải KNK như: Mơ hình phát thải Metan (MEM), Phát thải Metan từ hệ sinh thái lúa (MERES), Công cụ cân carbon EX-Ante (EX-ACT), mơ hình nơng nghiệp sử dụng đất (ALU)… Mơ hình DNDC mơ hình áp dụng rộng rãi giới để tính tốn phát thải từ ruộng lúa Nguyên lý mô dựa cacbon (C), nitơ (N) q trình sinh địa hóa hệ sinh thái nơng nghiệp Mơ hình sử dụng để dự đoán sinh trưởng trồng, nhiệt độ đất, độ ẩm, cacbon đất, phát thải khí nhà kính bao gồm N2O, NO, NH3, CH4 CO2 Viện Môi trường Nông nghiệp 89 Tạp chí Khoa học Cơng nghệ Nơng nghiệp Việt Nam - Số 06(127)/2021 Tính tốn phát thải KNK từ hoạt động sản xuất nơng nghiệp có vai trị quan trọng việc xác định cấu phát thải đề xuất giải pháp giảm phát thải KNK Tuy nhiên, nghiên cứu phát thải KNK canh tác lúa cho vùng Bắc Trung Bộ hạn chế Nghiên cứu có mục tiêu xây dựng đồ phát thải KNK dự báo phát thải canh tác lúa đất phù sa vùng Bắc Trung Bộ tới năm 2030 II VẬT LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU 2.1 Vật liệu nghiên cứu - Đất: Đất phù sa trồng Lúa vùng Bắc Trung Bộ - Lúa: Giống lúa Khang Dân 18 - Phân bón: NPK hỗn hợp 16:16:8 2.2 Phương pháp nghiên cứu 2.2.1 Mô hình DNCD Mơ hình DNDC mơ hình sinh địa hóa đất, cho phép dự báo lượng cacbon giữ lại đất, hàm lượng đạm bị phát thải số KNK CO2, CH4 từ hệ sinh thái nơng nghiệp Mơ hình xây dựng với thông số đầu vào gồm thông số tính chất lý hóa đất, điều kiện khí hậu nhiệt độ, độ ẩm, thơng số trồng lịch gieo trồng, thu hoạch, phương thức chăm bón (DNDC Guideline, 2012 University of New Hampshire) Cấu trúc mơ hình gồm: hợp phần phụ mơ hình khí hậu, đất, trồng mơ hình phân hủy dùng để đánh giá nhiệt độ, độ ẩm, oxy hóa khử đất q trình chuyển hóa đất, suất trồng, ước lượng hàm lượng cacbon đưa vào đất từ trồng Các thông số chịu tác động đặc trưng khí hậu, đất, trồng hoạt động người Hợp phần thứ hai gồm mơ hình nitrat hóa, khử nitrat mơ hình oxy hóa khử nhằm ước lượng phát thải khí CO2, CH4, NH3, NO, N2O, N2 từ hệ canh tác nơng nghiệp Mơ hình DNDC nhằm mơ lại mối quan hệ chu trình sinh địa hóa cacbon, nitơ yếu tố sinh thái (Hình 1) Hình Sơ đồ tính tốn lượng phát thải KNK dựa lựa chọn cacbon thấp mơ hình DNDC Các liệu đầu vào: Các liệu khí tượng thủy văn (nhiệt độ, lượng mưa, tốc độ gió, xạ mặt trời, độ ẩm); liệu canh tác (giống, thời gian gieo cấy, thu hoạch, bón phân, tưới nước, quản lý mùa vụ, cỏ hại); liệu đất đai (loại đất, pH, độ xốp, độ mặn, hàm lượng OC, NO3-, NH4+), kịch BĐKH Bộ Tài nguyên Môi trường 90 Các liệu đầu ra: Mô phát thải KNK (CH4, N2O) theo ngày năm suất tiềm trồng 2.2.2 Phương pháp hiệu chỉnh mô hình Để hiệu chỉnh mơ hình, mẫu khí phát thải từ ruộng lúa lấy vụ Xuân vụ Mùa từ 2015 - 2017 tỉnh Quảng Trị, Nghệ An Tạp chí Khoa học Cơng nghệ Nông nghiệp Việt Nam - Số 06(127)/2021 ừa iên - Huế phương pháp buồng kín Buồng thu khí có kích thước dài × rộng × cao 0,45 m × 0,4 m × 0,9 m Trong buồng thu khí có gắn nhiệt kế, hai quạt đảo khí để điều khiển áp suất, ống nhựa với đường kính 1,5 mm để trì áp suất khí bên bên ngồi hộp Mẫu khí lấy khoảng thời gian từ đến 10 sáng Mẫu khí rút xi lanh 60 mL sau 10 phút đóng nắp hộp (ở thời điểm 0, 10, 20, 30 phút) Sau lần lấy mẫu, ghi chép nhiệt độ hộp, thời gian lấy mẫu mực nước ruộng Các thông số đo đạc trường; mẫu khí thu mang phân tích để hiệu chỉnh mơ hình tham chiếu với kết mơ Mẫu khí lấy ruộng, tỉnh điểm với lần lặp lại, giai đoạn: bén rễ hồi xanh, bắt đầu đẻ nhánh, đẻ nhánh rộ, phân hóa địng, phát triển địng, trỗ bơng, chín sữa chín sáp Trong đó: Oi Pi: giá trị thực đo O P: giá trị trung bình chuỗi số liệu thực đo i = 1, 2, 3, , n: chiều dài mẫu Đánh giá kết mơ từ mơ hình thực thông qua hệ số xác định (R2 - Coe ent of determination) Giá trị R2 dao động từ đến 1, với giá trị gần cho thấy tương quan tốt liệu mô thực tế Với giá trị lớn 0,5 coi chấp nhận kết mơ 2.2.3 Phương pháp đồ, GIS Trong cách tiếp cận này, mơ hình DNDC, tích hợp GIS sử dụng cho mơ vùng Tồn sở liệu phân chia thành vùng mục tiêu tương ứng với đa giác ô lưới ứng với điều kiện khí hậu đất đai Độ phân giải hệ thống đa giác lưới phụ thuộc vào độ phân giải liệu thô liệu đầu vào Cơ sở liệu hỗ trợ mô khu vực DNDC bao gồm tệp liệu đồ chuyên đề đất, tệp liệu đồ chuyên đề khí hậu, Các tệp GIS chứa thơng tin vị trí, tệp khí hậu, thuộc tính đất, vùng, khu vực phương pháp quản lý canh tác Tiến trình mơ khu vực DNDC chạy cho hệ thống canh tác đơn vị lưới toàn miền mục tiêu Trong mô khu vực, DNDC đọc tất liệu đầu vào từ tệp GIS đặt trước cho lưới sau tổ chức lại thơng tin vào tệp để mơ (Hình 2) Hình Ngun lý tích hợp kết mơ hình DNDC vào GIS 91 Tạp chí Khoa học Công nghệ Nông nghiệp Việt Nam - Số 06(127)/2021 2.2.4 Phương pháp tính tốn phát thải KNK Tính tốn tổng lượng CH4, N2O phát thải vụ (Δd ngày từ cấy đến thu hoạch) theo công thức sau đây: CH4 (kg C/ha/vụ) = ∑CH4 (mg C/m2/hr) × 24 × Δd /100 N2O (kg N/ha/vụ) = ∑N2O (mg N/m2/hr) × 24 × Δd /100 Dựa vào cách tính IPCC (2007), tính tốn tiềm nóng lên tồn cầu thơng qua việc quy đổi tất loại khí CO2 tương đương (CO2e) Hệ số quy đổi CH4 CO2e = CH4 × 25; hệ số quy đổi N2O CO2e = N2O × 298 (IPCC, Climate Change 2007) Tổng lượng phát thải KNK từ ruộng lúa tính theo cơng thức sau: GWP = Phát thải CH4 × 25 + Phát thải N2O × 298 2.3 ời gian địa điểm nghiên cứu Nghiên cứu thực từ tháng 01 năm 2015 đến tháng 12 năm 2017 tỉnh Bắc Trung Bộ từ anh Hóa đến ừa iên Huế III KẾT QUẢ VÀ THẢO LUẬN 3.1 Kết kiểm định mơ hình Kết kiểm định mơ hình so sánh kết mô phát thải CH4 N2O với kết tính tốn trường để hiệu chỉnh số thơng số q trình sinh địa hóa đất sinh trưởng trồng Kết sau hiệu chỉnh mơ hình thể hình Hình Kết tốn phát thải CH4 N2O từ thực địa kết chạy mơ hình mơ Kết hình cho thấy, sau hiệu chỉnh kết tính tốn phát thái khí CH4 N2O mơ tính tốn thực địa có tương quan chặt Đối với CH4 hệ số tương quan R2 = 0,97 N2O có hệ số tương quan R2 = 0,85 3.2 Dự báo phát thải KNK vùng Bắc Trung Bộ tới năm 2030 Tổng diện tích đất gieo trồng lúa vùng Bắc Trung Bộ 702.700 ha, đất trồng lúa vụ Xuân 356.900 vụ Mùa 174.700 tổng diện tích đất trồng lúa diện tích đất phù sa 305.429 chiếm 85% so với diện tích đất trồng lúa vụ Xuân So sánh kết tính tốn lượng phát thải KNK tỉnh vùng Bắc Trung Bộ cho thấy, tỉnh giáp Bắc Bộ có lượng phát thải KNK vụ lúa Xuân cao vụ lúa Mùa ngược lại, tỉnh giáp Nam Trung Bộ lại có lượng phát thải KNK vụ lúa Mùa cao vụ lúa Xuân, tỉnh ừa iên - Huế có lượng phát thải KNK cao vùng Bắc Trung Bộ cho đơn vị diện tích vụ Xuân vụ Mùa, lượng phát thải vụ lúa hàng năm lên tới 11,51 CO2e/ha/năm (Hình 4) Hình Phát thải KNK trung bình giai đoạn 2005-2017 vùng Bắc Trung Bộ 92 Tạp chí Khoa học Công nghệ Nông nghiệp Việt Nam - Số 06(127)/2021 Kết tính tốn phát thải KNK dự báo đến năm 2030 dựa kịch phát thải thấp (RCP4.5) kịch phát thải cao (RCP8.5) đến năm 2030 Mơ mơ hình DNDC kết phát thải thể hình Hình Dự báo lượng phát thải KNK tỉnh vùng Bắc Trung Bộ năm 2030 So sánh kết dự báo lượng phát thải KNK đến năm 2030 tỉnh vùng Bắc Trung Bộ cho thấy, lượng phát thải chênh lệch nhiều vụ Xuân vụ Mùa tỉnh, vụ Xuân phát thải dao động 5,93 - 6,31 CO2e/ha/năm vụ lúa Mùa 4,96 - 6,31 CO2e/ha/năm Kết dự báo phát thải KNK thể tỉnh giáp Bắc Bộ có lượng phát thải KNK vụ lúa Xuân cao vụ lúa Mùa ngược lại, tỉnh giáp Nam Trung Bộ lại có lượng phát thải KNK theo xu hướng vụ lúa Mùa cao vụ lúa Xuân, tổng lượng phát thải vụ lúa năm 2030 dao động 10,89 - 12,53 CO2e/ha/năm 3.3 Xây dựng đồ dự báo phát thải vùng Bắc Trung Bộ cho canh tác lúa đất phù sa tới 2030 Trong giai đoạn 2015 - 2017: Phát thải KNK cho vụ lúa tỉnh vùng Bắc Trung Bộ dao động từ 9,125 đến 11,509 CO2e/ha/năm; cao tỉnh ừa iên - Huế thấp Nghệ An, phát thải tăng cao Nghệ An thấp ừa iên - Huế Phát thải năm dự báo cho năm 2030 tỉnh dao động từ 10,887 đến 12,526 tấn/ha/năm lượng phát thải năm 2030 tăng so với thời kỳ 2015 - 2017 phát thải tăng thấp tại ừa iên Huế (1,017 CO2e/ha/năm) cao Nghệ An (2,91 CO2e/ha/năm) (Hình 6) Hình So sánh lượng phát thải KNK trung bình/năm thời kỳ 2015 - 2017 2030 Lượng phát thải KNK tỉnh giáp Bắc Bộ vào vụ lúa Xuân cao vụ lúa Mùa ngược lại, tỉnh giáp Nam Trung Bộ lại có lượng phát thải KNK vụ lúa Mùa cao vụ lúa Xuân giai đoạn 2015 - 2017 năm dự báo 2030 (Hình 7) 93 Tạp chí Khoa học Cơng nghệ Nơng nghiệp Việt Nam - Số 06(127)/2021 Hình Bản đồ thể lượng phát thải KNK dự báo năm 2030 vùng Bắc Trung Bộ trung bình giai đoạn 2015 - 2017 IV KẾT LUẬN Bản đồ phát thải KNK biên tập, xây dựng đồ đơn tính bao gồm: đồ đất, đồ địa hình, đồ khí hậu, đồ trạng sử dụng đất tỷ lệ 1:50.000 Lớp thông tin phát thải KNK cập nhật trực tiếp đồ số chuẩn hóa thành sở liệu GIS lưu trữ định dạng SHP thuận lợi cho việc truy xuất tích hợp sở liệu với chương trình mơ hình hóa mơ Lượng phát thải KNK dự báo năm 2030 vụ lúa Xuân đất phù sa dao động 5,93 - 6,31 CO2e/ha/năm, tăng lên 0,63 - 1,25 CO2e/ha/năm so với trung bình năm giai đoạn 2005 - 2017 Lượng phát thải KNK dự báo đến năm 2030 vụ lúa Mùa đất phù sa dao động 4,96 - 6,31 CO2e/ha/năm, tăng lên 0,33 - 1,66 CO2e/ha/năm so với trung bình năm giai đoạn 2005 - 2017 Lượng phát thải KNK tỉnh giáp Bắc Bộ vào vụ lúa Xuân cao vụ lúa Mùa ngược lại, tỉnh giáp Nam Trung Bộ lại có lượng phát thải KNK vụ lúa Mùa 94 cao vụ lúa Xuân giai đoạn 2005 - 2017 năm dự báo 2030 Lượng KNK phát thải từ đất phù sa trồng lúa năm (2 vụ lúa) tỉnh vùng Bắc Trung Bộ trung bình năm giai đoạn 2015 - 2017 dao động 9,13 - 11,51 CO2e/ha, cao tỉnh ừa iên - Huế thấp Nghệ An Dự báo phát thải KNK cho năm 2030 dao động từ 10,89 - 12,53 CO2e/ha TÀI LIỆU THAM KHẢO 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 đồ Forster, P., V Ramaswamy, P Artaxo, T Berntsen, R Betts, D.W Fahey, J Haywood, J Lean, D.C Lowe, G Myhre, J Nganga, R Prinn, G Raga, M Schulz and R Van Dorland, 2007 Changes in Atmospheric Constituents and in Radiative Forcing In: Climate Change 2007: e Physical Science Basis Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change [Solomon, S., D Qin, M Manning, Tạp chí Khoa học Công nghệ Nông nghiệp Việt Nam - Số 06(127)/2021 Z Chen, M Marquis, K.B Averyt, M.Tignor and H.L Miller (eds.)] Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA: 129-234 IPCC, Climate Change, 2007 Synthesis Report https:// www.ipcc.ch/report/ar4/syr/ Liu J., H.J Zhang, L.S Pereira, 1998 Model validation and crop coe cients for irrigation schedualing in the north China planing Agricultural Water Manage, 36: 233-246 Ministry of Natural Resources and Environment, 2010 Vietnam’s second national communication to the united nations framework convention on climate change https://data.opendevelopmentmekong.net/ dataset/thong-bao-qu-c-gia-l-n-th-hai-cho-cong-u-ckhung-c-a-lien-h-p-qu-c-v-bi-n-d-i-khi-h-u Nguyen Hong Son, Bui i Phuong Loan, Tran Quoc Viet, Vu Duong Quynh, Syed Faiz-ul Islam, 2014 Báo cáo đo phát thải KNK từ ruộng lúa nước tỉnh Bình Định Quảng Bình, thuộc dự án SNV University of New Hampshire, 2012 User’s Guide for the DNDC Model version 9.5 Institute for the Study of Earth, Oceans and Space Mapping greenhouse gas emissions on paddy rice alluvial soils in the North Central Nguyen Van iet, Bui i Phuong Loan Abstract is study presents GHG emissions mapping using GIS techniques for the simulation results of CH4 and N2O (Greenhouse gas) emissions from rice cultivation alluvial soils of the North Central Region of Vietnam according to the De-nitri cation and Decomposition (DNDC) model e model was calibrated and validated with eld observation data in Quang Tri, Nghe An and ua ien Hue province e model was then used to estimate GHG emissions from rice elds in this area by 2030 using a compiled soil data, climate, land use and rice cultivation management database according to climate and sea-level rise in Viet Nam by the Ministry of Natural Resource and Environment in 2016 e results illustrated that the amount of GHGs emission in the summer season tended to be higher than that in the spring season in Quang Binh, Quang Trị and ua ien Hue Meanwhile, Ha Tinh, Nghe An and anh Hoa witnessed the higher GHGs emission in the spring season It is forecasted that by 2030 that ua ien - Hue will have the highest amount of GHGs emission from rice elds with 11.512 tons CO2 eq/ha However, rice elds in Nghe An are likely to emit the lowest amount of GHGs In addition, the GHGs emission from provinces in North Central Region tends to increase by the year 2030 Keywords: Alluvium soil, DNDC, greenhouse gas emissions Ngày nhận bài: 03/9/2020 Ngày phản biện: 15/9/2020 Người phản biện: PGS TS Phạm Quang Hà Ngày duyệt đăng: 02/10/2020 TIỀM NĂNG ỨC CHẾ THỰC VẬT CỦA CAO PHÂN ĐOẠN TỪ CÂY SÀI ĐẤT BA THÙY (Wedelia trilobata) Nguyễn Quốc Cường1, Đỗ Tấn Khang1, Nguyễn Văn Ây2, Trần anh Mến3, Phạm Văn Trọng Tính4, Trần Ngọc Quý4 TÓM TẮT Cây sài đất ba thùy (Wedelia trilobata (L.) Hitchc) nhiều nghiên cứu cho thấy có khả gây ức chế sinh học số loài thực vật hệ sinh thái Trong nghiên cứu này, tác động ức chế phân đoạn hexane, ethyl acetate, nước methanol thử nghiệm nảy mầm, sinh trưởng phát triển xà lách (Lactuca sativa L.) Kết cho thấy cao ethyl acetate có nồng độ từ 2,5 - 5,0 mg/mL làm giảm đáng kể tỷ lệ nảy mầm tiêu sinh trưởng phát triển hạt xà lách Kết định lượng diệp lục tố xử lý phân đoạn ethyl acetate làm giảm hàm lượng diệp lục tố a, b caratenoid nồng độ mg/mL 68,89%; 57,58% 52,63% so với đối chứng Hàm lượng polyphenol avonoid cao có phân đoạn Viện Nghiên cứu Phát triển Công nghệ Sinh học, Trường Đại học Cần Thơ Khoa Nông nghiệp, Trường Đại học Cần Thơ Khoa Khoa học Tự nhiên, Trường Đại học Cần Thơ Khoa Phát triển Nông thôn, Trường Đại học Cần Thơ 95 ... vùng Bắc Trung Bộ tới năm 2030 Tổng diện tích đất gieo trồng lúa vùng Bắc Trung Bộ 702.700 ha, đất trồng lúa vụ Xuân 356.900 vụ Mùa 174.700 tổng diện tích đất trồng lúa diện tích đất phù sa 305.429... 06(127)/2021 Hình Bản đồ thể lượng phát thải KNK dự báo năm 2030 vùng Bắc Trung Bộ trung bình giai đoạn 2015 - 2017 IV KẾT LUẬN Bản đồ phát thải KNK biên tập, xây dựng đồ đơn tính bao gồm: đồ đất, đồ địa... Nam Trung Bộ lại có lượng phát thải KNK vụ lúa Mùa 94 cao vụ lúa Xuân giai đoạn 2005 - 2017 năm dự báo 2030 Lượng KNK phát thải từ đất phù sa trồng lúa năm (2 vụ lúa) tỉnh vùng Bắc Trung Bộ trung