Sử dụng kỹ thuật 210pb trong tính toán tốc độ trầm tích tại khu vực ven biển cửa ba lạt (sông hồng)

Tuyển tập báo cáo Hội nghị Khoa học Công nghệ hạt nhân toàn quốc lần thứ 14 Proceedings of Vietnam conference on nuclear science and technology VINANST-14 SỬ DỤNG KỸ THUẬT 210Pb TRONG TÍNH TỐN TỐC ĐỘ TRẦM TÍCH TẠI KHU VỰC VEN BIỂN CỬA BA LẠT (SÔNG HỒNG) USE OF 210Pb TECHNIQUE FOR CALCULATING SEDIMENTATION RATE IN THE BA LAT COASTAL AREA (RED RIVER) BÙI ĐẮC DŨNG, DƯƠNG ĐỨC THẮNG, DƯƠNG VĂN THẮNG, VƯƠNG THU BẮC, NGUYỄN THỊ THU HÀ, ĐỒN THÚY HẬU, NGUYỄN THỊ OANH, LÊ ĐÌNH CƯỜNG Viện Khoa học Kỹ thuật hạt nhân, 179 Hoàng Quốc Việt – Cầu Giấy – Hà Nội Email: dacdung@gmail.com Tóm tắt: Thực Đề tài nghiên cứu cấp Bộ năm 2020-2021 “Nghiên cứu sử dụng kỹ thuật hạt nhân đồng vị để nhận diện nguồn gốc lịch sử phát sinh ô nhiễm hệ sinh thái vùng ven biển cửa Sông Hồng” lấy cột mẫu trầm tích khu vực ven biển cửa Ba Lạt (Sông Hồng) Các cột mẫu lấy khu vực Vườn Quốc gia Xuân Thủy, cắt lát thành lát mẫu đo xác định nhân phóng xạ U, Th, 226Ra, 137Cs 210Pb Từ kết phân tích chúng tơi sử dụng mơ hình tính tốn để xác định tuổi trầm tích tốc độ trầm tích cho cột mẫu Báo cáo trình bày kết áp dụng mơ hình tính tốn cột mẫu BL13 nhằm đưa phương pháp luận chung cho việc xác định tuổi cột trầm tích tốc độ trầm tích khu vực ven biển cửa Ba Lạt Kết cho thấy mơ hình tốc độ cung cấp khơng đổi (CRS) có lẽ mơ hình tốt để xác định tốc độ trầm tích Đối với cột mẫu BL13, việc sử dụng mơ hình hỗn hợp cách điều chỉnh lại kết mơ hình CRS kết hợp với mốc thời gian tháng năm 1971 (là thời gian xảy lũ lịch sử Sông Hồng) cho kết hợp lý Tốc độ trầm tích trung bình vị trí BL13 trước năm 1960 khoảng 0,5 cm/năm, từ năm 1960 đến tăng cao đáng kể, lên đến giá trị cao 1,5 cm/năm vào năm 1970, giảm xuống khoảng 1,0 đến 1,2 cm/năm từ năm 1980 đến Nguyên nhân thay đổi tốc độ trầm tích thay đổi dịng lũ lớn Sơng Hồng giai đoạn từ năm 1960 đến năm 1980 Từ khóa: cửa Ba Lạt, 210Pbex, CIC, CFCS, CRS, Mơ hình hỗn hợp Abstract: In implementing the 2020-2021 Ministerial-level research project "Research on the use of nuclear and isotope techniques to identify the origin and history of pollution in the coastal ecosystem of the Red River estuary", we have taken sediment cores in the coastal area of Ba Lat estuary (Red River) The cores were taken in Xuan Thuy National Park area, sliced into sample slices and numbered from the top to the bottom and then analyzing for radionuclides (U, Th, 226Ra, 137Cs and 210Pb Based on the analysis results, we have tried to use computational dating models to determine the age and sedimentation rate for each sediment core This report presents the results of applying computational models to the sediment core BL13 that could provide a general methodology for determining the age and the sedimentation rate of the sediment cores taken in the coastal area of Ba Lat estuary Results show that Constant rate of supply (CRS) model probably is the best model for determining the sedimentation rate For the sediment core BL13, the use of the Composite model by adjusting the results of the CRS model in combination with the time marker of August 1971 (the time of the historical flood event on the Red River) seems to give the most reasonable results The average sedimentation rate at the BL13 site before 1960 was about 0.5 cm/year, has increased significantly from 1960, up to the highest value of 1.5 cm/year in the 1970s and then decreased to about 1.0 to 1.2 cm/year from the 1980s to the present The change in the sedimentation rate could be caused by the main flow changes and large floods on the Red River in the years from 1960 to 1980 Keywords: Ba Lạt estuary, 210Pbex, CIC, CFCS, CRS, Composite Model MỞ ĐẦU Kỹ thuật hạt nhân kỹ thuật đồng vị ứng dụng rộng rãi để nghiên cứu ô nhiễm môi trường biển ô nhiễm kim loại nặng (KLN), chất dinh dưỡng, hữu trầm tích biển tiến hành rộng khắp giới [1-4] Để nghiên cứu lịch sử phát thải ô nhiễm môi trường biển, mẫu trầm tích phải lấy theo độ sâu dụng cụ, thiết bị lấy ngun cột trầm tích Sau cột trầm tích cắt lát mỏng, xử lý, phân tích tính tốn để xác định tuổi lát trầm tích Một phương pháp xác định tuổi trầm tích phương pháp xác định tuổi trầm tích kỹ thuật 210Pb Nguyên lý phương pháp xác định tuổi trầm tích kỹ thuật 210Pb tóm tắt sau: Đồng vị 210Pb số đồng vị dãy phân rã phóng xạ 238U Giữa 210Pb đồng vị mẹ 226Ra có cân khuếch tán đồng vị phóng xạ khí trung gian 222Rn từ mặt đất vào khơng khí bên Trong khí quyển, 222Rn phân rã thành 210Pb qua loạt đồng vị sống ngắn, sau 210Pb rơi trở lại mặt đất mặt nước Khi rơi lắng xuống mặt nước, 210Pb lắng xuống lòng hồ, đại dương với trầm tích Sự vượt trội 210Pb so với lượng 210Pb cân với 226Ra chỗ trầm tích ứng dụng để đánh giá tuổi trầm tích Lượng 210Pb vượt trội so với 226Ra thường ký hiệu 441 Tiểu ban D2: Ứng dụng kỹ thuật hạt nhân công nghiệp lĩnh vực khác Section D2: Application of nuclear techniques in industries and others Pb dư, 210Pbex (hoạt độ 210Pbex = hoạt độ 210Pb – hoạt độ 226Ra) Lớp trầm tích bề mặt tiếp xúc bùn/nước nhận bổ sung 210Pb từ khí bị lớp trầm tích phủ lên Khi khơng cịn bổ sung thêm 210Pb, lượng 210Pbex lớp trầm tích giảm theo quy luật hàm mũ với thời gian bán rã 22,26 năm Quy luật sử dụng để xác định tuổi lớp trầm tích biết giá trị hoạt độ ban đầu 210Pbex lớp trầm tích thời điểm hình thành 210 Đã có nhiều mơ hình tính tốn phát triển để xác định giá trị hoạt độ ban đầu 210Pbex qua xác định tuổi lát trầm tích tốc độ trầm tích [5] Kỹ thuật 210Pb dư cơng nhận áp dụng rộng rãi để đánh giá tốc độ lắng đọng trầm tích cho lịng hồ, đầm lầy, biển Kỹ thuật 210Pb dư trở thành phổ biến đóng vai trị quan trọng xác định tốc độ bồi lắng trầm tích đại khoảng 100 năm trở lại Trên sở tốc độ bồi lắng trầm tích theo thời gian khôi phục lịch sử ô nhiễm KLN chất ô nhiễm khác [6-9] Đề tài nghiên cứu cấp Bộ năm 2020-2021 “Nghiên cứu sử dụng kỹ thuật hạt nhân đồng vị để nhận diện nguồn gốc lịch sử phát sinh ô nhiễm hệ sinh thái vùng ven biển cửa Sông Hồng” lấy cột mẫu trầm tích khu vực ven biển cửa Ba Lạt (Sông Hồng) Các cột mẫu lấy khu vực Vườn Quốc gia Xuân Thủy, cắt lát thành lát mẫu dày cm đánh số từ xuống Các lát mẫu chuyển phịng thí nghiệm, đưa vào bát sứ để khô tự nhiên sau tiếp tục xử lý để phân tích Các lát mẫu phân tích xác định nhân phóng xạ U, Th, 226Ra, 137Cs phổ kế gamma [10] 210Pb (thông qua 210Po) phổ kế alpha [11] Từ kết phân tích chúng tơi sử dụng mơ hình tính tốn để xác định tuổi trầm tích tốc độ trầm tích cho cột mẫu Báo cáo trình bày kết áp dụng mơ hình tính tốn xác định tuổi trầm tích cột mẫu BL13 nhằm đưa phương pháp luận chung cho việc xác định tuổi tốc độ trầm tích cột trầm tích khu vực ven biển cửa Ba Lạt (Sông Hồng) NỘI DUNG 2.1 Đối tượng Phương pháp nghiên cứu Cột mẫu BL13 lấy khu vực ngập nước (Hình 1) Vườn Quốc gia Xuân Thủy, vùng bãi bồi bán ngập nước hình thành Sơng Hồng bồi tụ mở rộng diện tích phía biển cho đồng Bắc Bộ khu vực cửa Ba Lạt [12] Cột mẫu BL13 lấy dụng cụ lấy mẫu lõi trầm tích có đường kính 5,90 cm Các kết nghiên cứu trước [13-15] cho thấy tốc độ trầm tích trung bình Cửa Ba Lạt dao động khoảng 0,18 đến 3,04 cm/năm, trung bình 0,77 cm/năm Vì vậy, cột mẫu trầm lấy đến độ sâu 90-100 cm để có khoảng thời gian hình thành trầm tích vịng 100-120 năm trở lại Các cột mẫu cắt lát cm/lát để đảm bảo có 20 lát mẫu cột mẫu Các lát mẫu đánh số từ xuống theo độ sâu Các lát mẫu chuyển phòng thí nghiệm, đưa vào bát sứ để khơ tự nhiên sau tiếp tục xử lý phân tích Trong phịng thí nghiệm, lát mẫu phân tích xác định hoạt độ riêng nhân U, Th, 226Ra, 137Cs 210Pb Từ kết phân tích tính hoạt độ riêng 210Pb dư, 210Pbex (hoạt độ 210Pbex = hoạt độ Pb – hoạt độ 226Ra) lát mẫu Hoạt độ 210Pb dư sử dụng để tính tốn tuổi lát trầm tích mơ hình xác định tuổi trầm tích Các bước việc xác định tuổi trầm tích thực sau [9]: 210 a) Xây dựng đồ thị logarit tự nhiên hoạt độ 210Pb dư, ln(210Pbex), theo độ sâu khối lượng tích luỹ (mass accumulation depth - mi , kg m-2) b) Sử dụng Mơ hình thơng lượng tốc độ trầm tích khơng đồi (Constant Flux Constant Sedimentation - CFCS) để xác định tốc độ trầm tích trung bình Mơ hình sử dụng đồng thời 02 giả thiết mật độ 210Pbex rơi lắng bề mặt trầm tích (Bq m-2 năm-1) khơng đổi theo thời gian giả thiết tốc độ trầm tích khối lượng r (kg m-2 năm-1) số Khi lát trầm tích có hoạt độ 210Pbex ban đầu C0 (Bq kg-1) 442 Tuyển tập báo cáo Hội nghị Khoa học Công nghệ hạt nhân toàn quốc lần thứ 14 Proceedings of Vietnam conference on nuclear science and technology VINANST-14 BL13 Hình Vị trị lấy cột mẫu BL13 Do tương quan tuổi lát trầm tích t(i) độ sâu khối lượng tích luỹ mi xác định thơng qua qua phương trình: m t(i) = i (1) r Nên giả thiết mơ hình CFCS thỏa mãn, tăng độ sâu khối lượng tích luỹ, hoạt độ Pbex lát thứ i (Ci) giảm theo hàm số mũ theo phương trình: 210 λ Ci = C0 e−λt(i) = C0 e−λmi /r , hay ln(Ci) = ln(C0) − r mi (2) Trong đó: t(i) khoảng thời gian kể từ lát trầm tích thứ (i) hình thành (năm); λ số phân rã 210Pb 0.03118 năm-1; Co Hoạt độ 210Pbex lát trầm tích thứ i vào thời điểm lát hình thành (Bq kg-1); Ci Hoạt độ 210Pbex lát trầm tích thứ i (Bq kg-1) Ta giải phương trình (2) phương pháp hồi quy tuyến tính Đường hồi quy tuyến tính y = a + bx cắt trục tung a = ln(C0 ) hệ số góc đồ thị b = λ/ r Khi C0 = ea r = λ/b Độ không đảm bảo đo C0 r xác định từ độ không đảm bảo đo a b Tuổi lát trầm tích t(i) xác định thơng qua phương trình (1) Mơ hình sử dụng đồ thị thể mối quan hệ ln(Ci ) mi có dạng tuyến tính Khi đồ thị có dạng tuyến tính theo đoạn (có đoạn tuyến tính với hệ số góc khác nhau) ta xác định giá trị r trung bình đoạn đồ thị C0 từ giá trị r đoạn đồ thị tuyến tính trước Tuy vậy, cần phải thực phương pháp cách cẩn thận đồ thị có dạng tuyến tính theo đoạn có nghĩa mật độ 210Pbex rơi lắng bề mặt trầm tích C0 thay đổi c) Áp dụng mơ hình CFCS số lát cột mẫu để xác định hoạt độ 210Pbex ban đầu C0 Sử dụng Mô hình hoạt độ ban đầu khơng đổi (Constant Initial Concentration - CIC) để ước lượng tuổi lát trầm tích Tuổi lát trầm tích thứ (i) tích theo công thức sau: C t(i) = ln o (3) λ Ci 443 Tiểu ban D2: Ứng dụng kỹ thuật hạt nhân công nghiệp lĩnh vực khác Section D2: Application of nuclear techniques in industries and others Cần lưu ý 210Pbex có phân rã q trình hình thành lát trầm tích, hoạt độ 210Pbex trung bình lát thứ (C1 ) khơng phải lúc tỷ lệ với giá trị C0 Trên thực tế, xác định rõ giá trị C0 lớp trầm tích bề mặt Vấn đề giải cách áp dụng mơ hình CFCS số lát cột mẫu Đường hồi quy tuyến tính cắt trục tung điểm a = ln(C0 ), từ xác định C0 = ea giá trị sử dụng để xác định tuổi trầm tích phương trình (3) d) Áp dụng mơ hình CFCS số lát cột mẫu để xác định tổng mật độ tồn lưu A(j) 210Pbex (dưới lớp trầm tích thứ (j) cột trầm tích q ngắn đồ thị hoạt độ 210Pbex khơng đạt đến giá trị nền) Giả sử ta xác định giá trị r trung bình mơ hình CFCS, lớp đất sâu lõi trầm tích Từ kết đó, ta tính mật độ tồn lưu 210Pbex lớp trầm tích thứ (j) theo Phương trình sau: A(j)= r.C(j) λ Như vậy, tổng mật độ tồn lưu cột mẫu A(0) tính sau: A(0) = δA + A(j) (4) (5) Trong đó: δ(A) tổng mật độ tồn lưu từ bề mặt tới lát sâu lõi trầm tích Từ kết xác định A(0) xác định tuổi lát trầm tích theo Mơ hình tốc độ cung cấp khơng đổi (Constant rate of supply model - CRS) Giả thiết mơ hình giống giả thiết số mơ hình CFCS mật độ 210Pbex rơi lắng bề mặt trầm tích (Bq m-2 năm-1) không đổi theo thời gian, nhiên tốc độ trầm tích khối lượng r (kg m-2 năm-1) thay đổi Hoạt độ 210Pbex ban đầu tốc độ trầm tích khối lượng ri lớp trầm tích khác thay đổi ln tỷ lệ nghịch với Tại nhiều khu vực, gia tăng nguồn cung trầm tích (ví dụ, lũ lụt) khơng đồng nghĩa với gia tăng mật độ 210Pbex rơi lắng bề mặt trầm tích, hoạt độ 210Pbex thấp (do hiệu ứng pha loãng) Tuổi lát trầm tích thứ (i) tính theo cơng thức sau: A(0) λ A(i) t(i) = ln (6) Trong đó: A(0) tổng mật độ tồn lưu cột trầm tích (Bq m-2); A(i) tổng mật độ tồn lưu bên lát trầm tích thứ (i) (Bq m-2) e) Từ kết tính tốn tuổi lát trầm tích tính tốc độ trầm tích khối lượng lát trầm tích thứ (i) (ri , kg m-2 năm-1) sau: ri = λA(i) C(i) = Δm(i) Δt(i) (7) Và tốc độ trầm tích lát trầm tích thứ (i) (s, m năm-1) xác định sau: Si Δz(i) = Δt(i) (8) Trong đó: Δm(i) khối lượng lát mẫu (kg); Δz(i) bề dày lát mẫu (m); Δt(i) khoảng thời gian hình thành lát mẫu (năm), tính từ kết xác định tuổi lát trềm tích f) Sử dụng Mơ hình hỗn hợp (composite model) cách biến đổi mơ hình CRS (CM-CRS) Nếu biết xác hai mốc thời gian tham chiếu, kết tính tuổi kỹ thuật 210Pbex có sai lệch lớn so với mốc thời gian tham chiếu, ta sử dụng mơ hình CM-CRS để xác định thay đổi mật độ 210 Pbex rơi lắng bề mặt trầm tích theo thời gian, đưa [16] sau: Sự chênh lệch tuổi xác định từ mơ hình CRS mốc thời gian xác định từ đặc điểm địa tầng thể thay đổi tốc độ cung cấp 210Pbex địa điểm lấy cột mẫu trầm 444 Tuyển tập báo cáo Hội nghị Khoa học Công nghệ hạt nhân toàn quốc lần thứ 14 Proceedings of Vietnam conference on nuclear science and technology VINANST-14 tích Sự chênh lệch nhiều nguyên nhân như: kiện lũ lụt, sụt lún trầm tích, dịng nước xốy, thay đổi mục đích sử dụng đất Do đó, giai đoạn ngắn, tốc độ cung cấp 210Pbex cao nhiều so với tốc độ cung cấp giai đoạn lắng đọng trầm tích bình thường Trong trường hợp tốc độ cung cấp 210Pbex tăng vậy, mơ hình CRS CIC cho kết xác định tuổi khơng xác Khi có chênh lệch tuổi xác định từ mơ hình phương pháp nêu trên, tốc độ cung cấp Pbex tính từ mốc thời gian độc lập, xác định dấu hiệu đánh dấu địa tầng Gọi x x2 độ sâu hai lát lõi trầm tích biết tuổi t t2, mật độ rơi lắng 210Pbex trung bình giai đoạn tính cơng thức: 210 λΔA P = e−λt1 −e−λt2 (9) ΔA hiệu số mật độ tồn lưu 210Pbex x1 x2 So sánh với mật độ rơi lắng 210Pbex khơng khí suy đốn q trình dẫn đến thay đổi mật độ tồn lưu 210Pbex trầm tích Nếu giả thiết mật độ tồn lưu giai đoạn trầm tích khơng đổi, ta xác định tuổi trầm tích hiệu chỉnh tốc độ trầm tích áp dụng ngun tắc mơ hình CRS với mật độ tồn lưu 210 Pbex Từ phương trình mơ hình CRS giá trị P tính được, tính tuổi lớp trầm tích độ sâu x (sâu x1 nông x2) thông qua công thức sau: P −λt P e = e−λt1 + ΔA(x1 , x) (10) λ λ ΔA(x1 , x) hiệu số mật độ tồn lưu 210Pbex x1 x Vì để sử dụng cơng thức ta gán mốc thời gian t1 t2 ứng với độ sâu x1 x2, tuổi xác định lớp trầm tích nằm x1 x2 khơng có sai số lớn tốc độ cung cấp 210Pbex thay đổi Tốc độ trầm tích khối lượng lát trầm tích thứ (i) thời điểm t tính theo cơng thức sau: ri = Pe−λt (11) C(i) 2.2 Kết Kết phân tích xác định hoạt độ riêng nhân 226Ra, 137Cs 210Pb đưa Bảng Bảng Kết phân tích nhân phóng xạ cột mẫu BL13 STT 10 11 12 13 14 15 16 17 Tên lát mẫu BL13.01 BL13.02 BL13.03 BL13.04 BL13.05 BL13.06 BL13.07 BL13.08 BL13.09 BL13.10 BL13.11 BL13.12 BL13.13 BL13.14 BL13.15 BL13.16 BL13.17 Độ sâu lát mẫu (cm) 10 14 18 22 26 30 34 38 42 46 50 54 58 62 66 226 Ra (Bq/kg) 33.77 34.16 30.07 35.00 31.82 33.44 33.45 31.82 33.85 34.88 35.09 36.54 36.71 36.82 37.27 36.49 35.29 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 1.31 1.34 1.28 1.41 1.27 1.35 1.36 1.39 1.40 1.26 1.48 1.47 1.46 1.40 1.56 1.43 1.41 445 137Cs 1.44 1.30 1.44 1.39 1.51 1.78 1.41 1.59 1.45 1.36 1.09 1.38 1.57 1.48 1.66 1.41 1.34 (Bq/kg) ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 0.35 0.21 0.22 0.22 0.21 0.23 0.23 0.23 0.22 0.22 0.42 0.23 0.23 0.23 0.25 0.21 0.30 Pb tổng (Bq/kg) 210 67.23 67.18 59.13 62.96 54.26 51.76 53.81 50.12 49.66 45.78 48.80 46.34 46.82 44.43 47.79 48.74 49.15 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 2.11 2.11 2.29 2.19 3.17 1.95 2.96 2.01 3.11 1.94 3.31 3.25 2.80 1.89 2.06 1.99 3.10 Tiểu ban D2: Ứng dụng kỹ thuật hạt nhân công nghiệp lĩnh vực khác Section D2: Application of nuclear techniques in industries and others 18 19 20 21 22 23 24 BL13.18 BL13.19 BL13.20 BL13.21 BL13.22 BL13.23 BL13.24 70 74 78 82 86 90 94 36.40 34.22 36.17 34.84 36.18 36.30 35.94 ± ± ± ± ± ± ± 1.50 1.40 1.40 1.43 1.13 1.52 1.38 1.44 1.20 1.06 1.39 1.25 1.28 0.98 ± ± ± ± ± ± ± 0.23 0.32 0.36 0.22 0.16 0.24 0.19 49.42 42.85 48.61 40.17 43.16 38.91 39.23 ± ± ± ± ± ± ± 2.05 1.88 3.13 1.92 1.64 1.68 1.83 Hình Hình đồ thị biểu diễn thăng giáng hoạt độ 210Pbex 237Cs theo độ sâu Trên Hình nhận thấy hoạt độ 210Pbex giảm dần từ bề mặt đến độ sâu 54 cm (tương ứng với lát mẫu BL13.14), sau lại tăng giảm khơng có quy luật đạt gần giá trị (gần sai số) độ sâu 90-94 cm Trên đồ thị hoạt độ 237Cs theo độ sâu (Hình 3) nhận biết 02 đỉnh hoạt độ 237 Cs cao độ sâu 22 58 cm, tương ứng với lát mẫu BL13.06 BL13.15 Hình đồ thị biểu diễn thăng giáng logarit tự nhiên hoạt độ 210Pb dư, ln(210Pbex), theo độ sâu khối lượng tích luỹ kết áp dụng mơ hình CFCS cho cột mẫu BL13 Có thể thấy đồ thị thể mối quan hệ ln(Ci ) mi có hai đoạn tuyến tính với hệ số góc khác nhau, thể hai giai đoạn trầm tích với tốc độ trầm tích khác Đoạn tuyến tính (từ bề mặt đến độ sâu 58 cm, m i = 65,66 g cm-2) cho giá trị C0 37,88 Bq kg-1 tốc độ trầm tích khối lượng trung bình 1,34 g cm -2 năm-1 Đoạn tuyến tính bên cho giá trị tốc độ trầm tích khối lượng trung bình 0,78 g cm -2 năm-1, thấp gần nửa so với đoạn phía Sử dụng phương trình (1) cho giai đoạn trầm tích để xác định tuổi lát trầm tích, kết đưa Bảng Kết xác định giá trị C0 37,88 Bq kg-1 mơ hình CFCS cho đoạn tuyến tính bên (Hình 4) sử dụng mơ hình CIC để tính tuổi lát trầm tích theo phương trình (3) Kết đưa Bảng Từ kết áp dụng mơ hình CFCS cho đoạn tuyến tính bên (Hình 4) xác định giá trị tốc độ trầm tích khối lượng trung bình 0,78 g cm-2 năm-1 Kết tính tốn tổng mật độ tồn lưu A(24) 210Pbex (dưới lát trầm tích thứ 24) 372 Bq m-2 tổng mật độ tồn lưu cột trầm tích A(0) 17220 Bq m -2 Áp dụng phương trình (6) tính tuổi lát trầm tích, kết đưa Bảng 446 Tuyển tập báo cáo Hội nghị Khoa học Công nghệ hạt nhân toàn quốc lần thứ 14 Proceedings of Vietnam conference on nuclear science and technology VINANST-14 Hình Hoạt độ 237Cs theo độ sâu Hình Hoạt độ 210𝑃𝑏𝑒𝑥 theo độ sâu Hình Thăng giáng ln(210Pbex) theo độ sâu khối lượng tích luỹ kết áp dụng mơ hình CFCS cho cột mẫu BL13 Ngồi ra, đề tài thử nghiệm sử dụng mơ hình CM-CRS để điều chỉnh lại kết mơ hình CRS theo mốc thời gian biết Đối với cột mẫu BL13, đề tài thử nghiệm sử dụng mốc thời gian sau: - Hai đỉnh hoạt độ riêng 137Cs cao, tương ứng với năm 1986 1963 năm ghi nhận có thơng lượng rơi lắng 137Cs cao (Hình 3) - Đỉnh có tốc độ trầm tích cao từ Mơ hình CRS (Hình 5), tương ứng với mốc thời gian tháng năm 1971, thời gian xảy lũ lịch sử Sông Hồng [15], nguyên nhân làm tăng đột biến lượng trầm tích đổ cửa Ba Lạt Áp dụng phương trình (10) tính tuổi lát trầm tích cột mẫu BL13, kết đưa Bảng 447 ... ban đầu 210Pbex qua xác định tuổi lát trầm tích tốc độ trầm tích [5] Kỹ thuật 210Pb dư công nhận áp dụng rộng rãi để đánh giá tốc độ lắng đọng trầm tích cho lịng hồ, đầm lầy, biển Kỹ thuật 210Pb. .. “Nghiên cứu sử dụng kỹ thuật hạt nhân đồng vị để nhận diện nguồn gốc lịch sử phát sinh ô nhiễm hệ sinh thái vùng ven biển cửa Sông Hồng” lấy cột mẫu trầm tích khu vực ven biển cửa Ba Lạt (Sơng Hồng). .. 226Ra, 137Cs 210Pb Từ kết phân tích tính hoạt độ riêng 210Pb dư, 210Pbex (hoạt độ 210Pbex = hoạt độ Pb – hoạt độ 226Ra) lát mẫu Hoạt độ 210Pb dư sử dụng để tính tốn tuổi lát trầm tích mơ hình