Luận văn Sự biến thiên của front nhiệt mặt biển trong biển đông

1 MỤC LỤC MỞ ĐẦU .................................................................................................. 2 CHƢƠNG 1: GIỚI THIỆU CHUNG..................................................... 3 1.1 Khái niệm................................................................................. 3 1.2 Các loại front........................................................................... 5 CHƢƠNG 2: TỔNG QUAN VỀ NGHIÊN CỨU FRONT .................. 8 CHƢƠNG 3: HỆ THỐNG VỆ TINH VÀ DỮ LIỆU ẢNH VỆ TINH 16 3.1 Hệ thống các vệ tinh NOAA-AVHRR ................................. 3.2 Cơ sở dữ liệu ảnh NOAA-AVHRR ...................................... 16 17 CHƢƠNG 4: PHƢƠNG PHÁP XÁC ĐỊNH FRONT TỪ ẢNH VỆ TINH ........................................................................................................ 19 4.1 Thuật toán xác định ranh giới các ảnh nhiệt mặt biển (SST) của Cayula & Cornillon (1992) .................................................. 19 4.2 Xác định ranh giới sử dụng nhiều ảnh SST .............................. 4.3 Tính toán Gradient ..................................................................... 4.4 Quy trình tính toán front nhiệt mặt biển với phƣơng pháp của Cayula & Cornillon (1992) .................................................. 24 25 27 CHƢƠNG 5: XÁC ĐỊNH FRONT NHIỆT MẶT BIỂN KHU VỰC BIỂN ĐÔNG ............................................................................................ 42 5.1 Tiền xử lý ảnh vệ tinh ………………………………………… 5.2 Xác định biến thiên gradient nhiệt mặt biển bằng phƣơng pháp tính toán gradient …………............................................. 42 46 5.3 Nghiên cứu xác định vị trí, tần suất xuất hiện của front ....... 51 KẾT LUẬN .............................................................................................. 72 TÀI LIỆU THAM KHẢO ...................................................................... 73 2 MỞ ĐẦU Front đại dương là một trong những yếu tố rất quan trọng trong hải dương học. Front nhiệt mặt biển là một dạng cơ bản nhất của front đại dương. Trong tất cả các đại dương cũng như các biển, luôn luôn tồn tại sự chênh lệch nhiệt độ theo các khu vực riêng biệt. Sự chênh lệch này có thể là do sự xâm nhập của các khối nước, quá trình vận chuyển nước của các hoàn lưu hay do sự khác biệt của các yếu tố vật lý, dinh dưỡng. Chính vì thế, nghiên cứu front nhiệt có thể giúp xác định các điều kiện vật lý khác biệt của các khu vực nước rộng lớn, nhờ đó có thể xác định ranh giới của các khối nước hay ranh giới của các dòng chảy lớn... Nghiên cứu front rất quan trọng trong việc nghiên cứu cũng như kết hợp hoàn thiện các trường hải dương khác. Khu vực Biển Đông là khu vực có đặc trưng gió mùa: gió mùa Đông Bắc vào mùa đông và gió mùa Tây Nam vào mùa hè. Trường nhiệt trong Biển Đông cũng có sự thay đổi theo mùa: vào mùa đông, xuất hiện một lưỡi nước lạnh từ phía Bắc xâm nhập sâu xuống phía Nam; mùa hè thường xuất hiện một khu vực nước trồi ở khu vực Nam Trung Bộ. Hai hiện tượng đặc trưng mùa này làm thay đổi rõ rệt trường nhiệt trong Biển Đông. Chính điều này, front trong Biển Đông cũng có sự thay đổi về cường độ và vị trí theo mùa. Khu vực Nam Trung Bộ là khu vực có sự thay đổi nhiệt độ mặt biển rõ rệt cả hai xu thế của trường nhiệt trong Biển Đông. Do cơ sở dữ liệu nhiệt mặt biển từ nguồn dữ liệu vệ tinh rất phong phú nên luận văn này tập trung nghiên cứu xác định front trong toàn khu vực biển Đông (99E – 121E; 1N – 25 N ) với sự tập trung vào nghiên cứu xác định front khu vực Nam Trung Bộ. Sự mở rộng khu vực nghiên cứu này cho cái nhìn tổng quát và rộng mở hơn về sự biến thiên của front nhiệt mặt biển trong Biển Đông. 3 CHƢƠNG 1: GIỚI THIỆU CHUNG 1.1 Khái niệm Front đại dương là những dải tương đối hẹp ngăn cách những vùng không gian rộng của các khối nước khác nhau hay các khu vực có cấu trúc thẳng đứng khác nhau. Các front hầu như luôn đi kèm với sự tăng cường gradient ngang của các yếu tố như: nhiệt độ, độ muối, mật độ, chất dinh dưỡng và các đặc trưng khác (Fedorov, 1986; Belkin, 2003[2]). Các front và các dòng chảy liên quan rất quan trọng trong vận chuyển nhiệt, muối, tương tác biển – khí quyển và các hoạt động sinh thái khác [1]. Front xuất hiện với các quy mô khác nhau: từ vài trăm mét cho tới hàng nghìn kilomet. Một số front xuất hiện trong thời gian ngắn nhưng hầu hết đều tồn tại tựa ổn định hoặc tồn tại theo mùa. Các front thường xuất hiện hoặc biến mất ở cùng vị trí trong cùng một mùa trong năm. Hầu như tất cả các front lớn thường xuất hiện quanh năm. Sự biến thiên nhiệt, muối qua các front mạnh nhất có thể cao từ 10 0 C - 15 0 C và từ 2‰ - 3‰, tuy nhiên nhìn chung thường từ khoảng 50C và 1‰. Độ rộng của front biến thiên rất khác nhau, có những nơi front rất hẹp (dưới 100 m) nhưng cũng có các front lớn rộng từ 50 - 20 km. Front có thể mở rộng tới độ sâu hàng trăm mét; một số front lớn phát triển tớ độ sâu 2000 m. Các front chủ yếu xuất hiện trong các quá trình: tương tác đại dương – khí quyển hoặc đại dương với băng và đáy biển [1]: - Các front có liên quan đến dòng chảy do vậy hình dạng front thể hiện hình dạng dòng chảy - Theo dòng chảy dọc front có thể tính toán tổng thông lượng nước/nhiệt/muối - Các front phân chia các khối nước khác nhau và tạo nên các chu trình tương ứng với tổng thông lượng nước, nhiệt và muối dọc theo front và theo kinh độ. 4 - Thông thường, các front thường trùng khớp với các biên sinh – địa lý lớn với các dải tăng cường sinh khối. - Các thông lượng nhiệt trên mặt, ứng suất gió và các điều kiện khí tượng khác có thể khác biệt lớn giữa các khu vực lạnh và nóng của front. Các front tương tác mạnh mẽ với các lớp biên khí quyển và phân chia các khu vực có sự khác nhau các tác động của các lực khí quyển, vì thế chúng rất quan trọng đối với việc quan trắc và dự báo thời tiết biển. - Các front ở vĩ độ cao liên quan trực tiếp với các điều kiện băng biển, vì thế vị trí của các front có thể xác định quy mô cực đại của băng biển. - Các front có ảnh hưởng lớn đến môi trường sóng âm, do vậy khi nghiên cứu quá trình truyền âm trong đại dương cũng cần phải biết được vị trí và đặc trưng của các front. - Cơ chế vận chuyển trầm tích được xác định chủ yếu thông qua hoàn lưu, do vậy việc giải đoán các thông tin cổ xưa về đại dương và khí quyển yêu cầu trước tiên là cần phải biết về vị trí của front hiện đại. - Do các front có liên quan với các dòng chảy hội tụ, sự ô nhiễm đại dương và cửa sông có thể tập trung với nồng độ lớn gấp hàng nghìn lần ở các dải front, do vậy đe dọa đời sống của các loài cá, động vật sinh sống trong các dải front. Các dữ liệu nhiệt độ mặt biển (SST) thu được từ vệ tinh đã được sử dụng để nghiên cứu các front ở Thái Bình Dương từ những năm 1970 (Legeckis, 1978). Các nghiên cứu trước đây được tập trung chủ yếu vào các front liên quan với các dòng chảy biên phía Tây như Kuroshio (Qiu và cs, 1990) và dòng phía Đông Australia (Nilsson và Cresswell, 1980). Các front gắn với các dòng chảy biên phía Đông và các khu vực nước trồi ven bờ đã được nghiên cứu ở các khu vực giống như Dòng chảy California (Strub và cs, 1991; Strub và James, 1995) và dòng chảy Peru – Chile (Brink và cs, 1983; Fonseca, 1989) [2]. 5 Các nghiên cứu trên rất quan trọng trong việc nghiên cứu các tính chất vật lý và địa lý của các đặc trưng khu vực sử dụng các phương pháp khác nhau; Hầu hết các nghiên cứu trên đều thực hiện trên các bộ số liệu có thời gian tương đối hạn chế. 1.2 Các loại front Các front được hình thành bởi các quá trình khác nhau. Theo các quá trình đó, chúng ta có thể chia thành các loại front sau [1]: - Front cửa sông: Những front này hình thành trên ranh giới giữa dòng nước ngọt từ sông đổ ra và nước biển. Do vậy, chúng chủ yếu là front độ muối mặc dù trong hầu hết các trường hợp có tồn tại chênh lệch nhiệt độ ngang front độ muối. Khi các dòng sông mang trầm tích lớn ra cửa sông, hầu hết các front lúc này là front độ đục; các front này có thể quan trắc dễ dàng từ không gian do các gradient có màu riêng biệt các dải front này. - Front thềm lục địa: Front thềm lục địa có thể quan trắc được tại những khu vực giữa thềm lục địa, ở phía trong của đứt gãy thềm lục địa. Những front này có nguồn gốc liên quan đến các dòng chảy ngược chiều trên thềm lục địa với các dòng chảy xa bờ. Điển hình của những dòng chảy này là dòng chảy phía Bắc (Bắc bán cầu) và dòng chảy ven bờ chảy xuống phía Nam (như trong vịnh Nam Địa Trung Hải). Ví dụ về các dòng chảy điển hình: Các front Biển Bering, Vịnh Nam Địa Trung Hải hay các front ở biển Đông Trung Hoa. - Front do sự xáo trộn thủy triều: Trong các khu vực nước nông, thủy triều bị tiêu tán một lượng lớn năng lượng do ma sát đáy. Trong các khu vực này, cột nước có thể được xáo trộn hoàn toàn do thủy triều. Độ sâu cột nước cực đại thủy triều có thể xáo trộn hoàn toàn là xấp xỉ 100m và trung bình là gần 50m. Chính nguyên nhân này, các front có thể hình thành giữa các khối nước xáo trộn toàn bộ ở khu vực nước nông và các khối nước phân tầng ở các độ sâu lớn hơn. Những front này có thể quan trắc ở nhiều thềm lục địa như: Biển Bắc, Biển Celtic, Bering, Đông Trung Hoa, Biển Hoàng Hải và Biển Okhotsk. Front do sự xáo trộn triều cũng có thể hình thành xung quanh các bờ biển, các hòn đảo (Đảo Pribilof trong Biển 6 Bering) và các bán đảo (Bán đảo Shandong, Biển Hoàng hải hay bán đảo Valde’s, thềm lục đại Patagonian). - Các front tại sườn lục địa: Đây là dạng front phổ biến nhất, những front này thường nằm trùng vị trí với sườn lục địa và ngăn các các khối nước ven bờ với khối nước xa bờ. Những front này là những front ngăn cách hai khối nước riêng biệt, gần bờ và xa bờ. - Các front nước trồi ven bờ: Các khối nước trồi (do gió, địa hình hay do thủy triều) mang nước từ dưới tầng sâu lên bề mặt. Do các khối nước phía dưới bề mặt thường lạnh hơn nước trên mặt biển nên các front xuất hiện giữa các khối nước trồi ven bờ lạnh với khối nước xa bờ ấm hơn. Những front đặc trưng có thể quan trắc được ở ngoài khơi Washington-Oregon-California, Peru-Chile, Tây Bắc Phi, Angola-Namibia-Nam Phi, Bán đảo Yucatan.... - Các front nước trồi xích đạo: Lực coriolis gây nên sự phân kỳ ở xích đạo với các khối nước mặt di chuyển xa Xích đạo trong cả hai bán cầu. Nước ở dưới sâu lạnh hơn đi lên trên mặt biển tạo nên một dải gradient nhiệt độ hay độ muối tại nơi giao lưu với nước ấm trên mặt biển. Những front này chủ yếu chiếm ưu thế ở trung tâm và phía Đông Thái Bình Dương và Đại Tây Dương. Sự bất đối xứng của bờ tại các biên Ấn Độ Dương có thể là nguyên nhân chính giải thích tại sao front nước trồi Xích đạo không phổ biến trong Ấn Độ Dương. - Các front dòng chảy biên phía Tây: Các front liên quan đến các dòng Gulf Stream, Kuroshio, Agulhas, Brazil và Đông Australia mạnh nhất thế giới. Các front này có dải biến thiên ngang của các đặc trưng hải dương học, độ mở rộng theo chiều sâu, vận tốc dòng chảy dọc front và xuất vận chuyển cũng như các đặc tính khác. Những front này có thể xâm nhập sâu trong các biển mở hàng nghìn km. Vào mùa đông, hai front khác biệt ở hai phía của dòng chảy biên phía Tây bởi vì mỗi dòng chảy mang theo nước nhiệt đới ấm hơn nước xung quanh các dòng chảy gần bờ và xa bờ. - Các front hội tụ cận nhiệt đới: Những front này hình thành do hội tụ gió Ekman mang nước với các nhiệt độ khác nhau (bởi vì gradient nhiệt quy mô lớn 7 Bắc – Nam trong đại dương) và giữ tại một khu vực tạo thành front. Những front này có thể quan trắc ở Biển Sargasso (Bắc Địa Trung Hải), Bắc Thái Bình Dương và phần phía Nam của Địa Trung Hải, Ấn Độ Dương và Thái Bình Dương. - Các front ở gần các khu vực băng đá: Những front này có liên quan đến các rìa băng biển và các quá trình riêng trong dải này. Trong suốt quá trình hình thành băng biển, quá trình giải phóng muối gây nên hiện tượng đối lưu muối biển. Trong quá trình tan băng, nước ngọt trên mặt tạo lên front muối khu vực giữa các khối nước do băng tan độ muối thấp và nước xung quanh. Vào mùa xuân, khi các khối nước do băng tan hấp thụ bức xạ mặt trời, một gradient nhiệt độ phát triển ngang front muối ban đầu. Những ví dụ điển hình của loại front này là: Labrador, Greenland, Barrents, Biển Bering và Nam Băng Dương. - Các front đại dương phía Nam: Các front này chỉ xuất hiện ở đại dương gần cực như: Cận Nam cực, Cực và các front do dòng chảy gần Nam cực. Nguồn gốc của các front trên liên quan đến các quá trình động lực của dòng chảy cận Nam cực. Các front sườn lục địa Nam cực hình thành giữa các khối nước thềm lục địa Nam cực và các khối nước ngoài khơi. Chính điều này lý giải việc cũng có thể coi front này là front tại sườn lục địa. Tuy nhiên, cấu trúc và quá trình động lực của front này khác xa so với front tại nơi đứt gãy thềm lục địa. Gió thổi từ lục địa Nam cực giữ duy trì dòng chảy hướng Tây trên thềm lục địa Nam cực đóng vai trò quan trọng trong quá trình tồn tại của front sườn lục địa Nam cực. 8 CHƢƠNG 2: TỔNG QUAN VỀ NGHIÊN CỨU FRONT Front là các dải có gradient cao, do vậy hầu hết các phương pháp xác định front tiến hành trên máy tính đều dựa trên các tính toán gradient (Kazmin và Rienecker, 1996; Yuan và Talley, 1996; Nakamura và Kazmin, 2003 [6]). Cách tiếp cận trong nghiên cứu này dựa trên phân tích phân bố cấp độ xám (histogram). Do mỗi front phân chia hai khối nước tương đối đồng nhất, do vậy phân bố cấp độ xám của bất kỳ một yếu tố hải dương học nào như SST, khu vực gần front có hai chế độ (mode) tần số tương ứng với các khối nước ngăn cách bởi front trong khi khối nước sau tương ứng với tần số thấp nhất giữa hai cấp độ (mode). Phân tách và xác định front được thể hiện ở 3 cấp độ: cửa sổ, ảnh và kết quả của việc chồng lấp các ảnh. Kích thước cửa sổ tối ưu được xác định thực nghiệm là 32 x 32 pixel đối với ảnh AVHRR độ phân giải 9 km (Cayula và Cornillon, 1992). Thuật toán xác định front được sử dụng dựa trên các giá trị SST của tất cả các pixel trong mỗi cửa sổ để tính toán histogram của tần số SST đối với 1 cửa sổ cho trước. Đối với mỗi cửa sổ chứa front, histogram SST tương ứng có tần số nhỏ nhất sẽ xác định front. Ý tưởng cơ bản này đã được thực hiện bởi Cayula và cs (1991), Cayula và Cornillon (1992,1995,1996) và Ullman và Cornillon (1999,2000,2001). Giai đoạn 1985 – 1996, các front đã thu được từ việc phân tích các trường SST Pathfinder NOAA/NASA. Những trường này thu được từ các số liệu bao phủ toàn cầu độ phân giải cao (AVHRR) (2 trường có độ phân giải 9.28 km một ngày) có sẵn từ phòng thí nghiệm Jet Propulsion. Các front nhiệt độ bề mặt thu được từ các trường nhiệt độ mặt biển đã loại bỏ mây với thuật toán tách ranh giới đa ảnh (Cayula và Cornillon, 1996; Ullman và Cornillon, 1999, 2000, 2001). Thuật toán tách front và tách mây đã được ứng dụng với ảnh nhiệt độ mặt biển (SST) trung bình tháng của bộ số liệu 22 năm. Số liệu front được tập hợp theo từng tháng (ví dụ: nhóm các số liệu của 22 tháng 1 với nhau). 9 Biển Nam Trung Hoa (Biển Đông Việt Nam) [2] Biển Nam Trung Hoa đặc trưng bởi số lượng lớn các front và phần lớn các front này nằm trùng vị trí với sườn lục địa. Ít nhất có 1 front (phía Tây Luzon) mở rộng xuống tới các lớp nước sâu. Front phía Tây Luzon trong trường nhiệt nằm trùng với sự nở hoa của các thực vật nổi (phytoplankton) quan trắc được từ các số liệu CZCS xuất hiện chủ yếu liên quan đến hiện tượng nước trồi (Tang & cs, 1999) hoặc liên quan đến xoáy thuận ổn định ngoài khơi Tây Bắc Luzon (Qu, 2000). Front này và dải front ngoài khơi phía Đông Việt Nam trùng với 3 khu vực hoạt động động lực mạnh có thể quan trắc được từ dữ liệu độ cao mực biển TOPEX/Poseidon (Ho & cs, 2000). Thông qua các dữ liệu TOPEX/Poseidon Wang và Cs (2000) đã chứng minh được rằng Front ven bờ Việt Nam gần với các dải hẹp có sự biến thiên theo mùa đáng kể với cực trị biến thiên lớn nhất. Chủ yếu các front nhiệt bề mặt (SST) thời kỳ mùa đông dọc theo sườn lục địa ngoài khơi phía Nam Trung Hoa và ngang cửa vịnh Thái Lan và Tây Bắc đảo Luzon. Hình 1: Các front SST chu kỳ dài biển Đông tháng hai giai đoạn 1985 – 1996[2] 10 Khu vực này có chế độ biến thiên theo mùa rất mạnh của các front mùa chủ yếu liên quan đến hiện tượng gió mùa. Các thành phần biến thiên mùa khác nhau rõ rệt ngang theo cả biển Đông. Các front đặc trưng trong biển Đông là: Front ven bờ Nam Trung Hoa (South China Coastal Front – SCCF) ở phía Bắc biển Đông, Front Vịnh Bắc Bộ (Bakbo Bay Front – BBF) giữa Việt Nam và đảo Hải Nam, Front Vịnh Thái Lan (Gulf of Thailand Front – GTF) và front phía Tây đảo Luzon (West of Luzon – WLF). Năm 2005, Belkin I.M đã hoàn thiện sơ đồ các dải front lớn trong Biển Đông (Hình 2). Hình 2: Sơ đồ phân bố front trong biển Đông (theo Belkin I.M)[6] Front ven bờ Nam Trung Hoa (SCCF) chủ yếu xuất hiện trong mùa đông (từ tháng 10 – 11 cho tới tháng 4 – 5) và biến mất trong mùa hè (Từ tháng 6 – 7 tới tháng 9 – 11). Tương tự như SCCF, front trong vịnh Bắc Bộ (BBF) chủ yếu xuất hiện trong mùa đông (từ tháng 11 – 12 cho tới tháng 4 - 5) và biến mất hoàn toàn trong mùa hè (từ tháng 5 – 6 đến tháng 10 – 11). Năm 2001, Wang D. cũng đã 11 nghiên cứu front với cơ sở dữ liệu từ năm 1985 – 1996 cũng đã đưa ra được bức tranh phân bố front theo các mùa trong năm (hình 3) Hình 3: Bản đồ phân bố tần suất xuất hiện front theo mùa từ năm 1991 – 1998 A: Mùa đông, B: mùa xuân, C: mùa hè và D: mùa thu [10] Các front xuất hiện đồng thời khắp nơi trong mùa đông ở phía Bắc biển Đông trong khi gần như không xuất hiện ở khu vực phía Nam ngoại trừ ở cửa vịnh Thái Lan. Cửa vịnh Thái Lan là nơi các front sườn lục địa theo mùa phát triển mạnh từ tháng 1 – tháng 2. Front phía Tây Luzon (WLF) quan trắc được ở phía Tây đảo Luzon (mở rộng xuống các khối nước sâu) trong mùa đông (từ tháng 10 – 11 đến tháng 4 – 5). Front này mở rộng theo hướng Bắc theo tiến trình mùa cho tới Eo Bashi và biến mất trong mùa hè (từ tháng 5 – 6 tới tháng 9 – 10). Front ngoài khơi Tây Bắc Boneo xuất hiện trong mùa xuân, mùa hè và mùa thu và biến mất trong mùa đông. So với các biển khác, các front trong biển Đông được nghiên cứu khá ít, tuy nhiên mối quan hệ của các front với nước trồi do gió, xáo trộn triều và lưu lượng nước sông đổ vào đã được nghiên cứu trong một số công trình của Kwan (1978), Kester và Fox (1993), Huang và Cs (1994), Li (1996), Tang & Ni (1996), su (1998), Tang & Cs (1999). 12 Vịnh Thái Lan [2] Vịnh Thái Lan là một vịnh nông (độ sâu cực đại nhỏ hơn 80 m) với diện tích khoảng 350.000 km2. Vịnh nhận được một lượng lớn nước lợ chủ yếu từ sông Mekong với lưu lượng dòng chảy đạt đỉnh trên 30.000 m3/s từ tháng 9 đến tháng 12, do vậy phân bố mật độ chủ yếu bị ảnh hưởng do phân bố độ muối và về cơ bản Vịnh chịu ảnh hưởng của cửa sông (Wattayakorn và cs, 1998; Yanagi và Takao, 1998). Ở cửa vịnh, một front muối phát triển mạnh trong mùa đông (từ tháng 12 – tháng 1) với khoảng biến thiên độ muối trên mặt dS = 32.0 – 33.8 ‰ (Yanagi và Takao, 1998; Yanagi & cs, 2001). Một front nhiệt liên quan đến front muối (gọi là Front Vịnh Thái Lan – Gulf of Thailand Front, GTF) phát triển từ tháng 12 – tháng 2 xuất hiện rõ rệt thông qua dữ liệu AVHRR và dữ liệu hải văn từ tháng 3 – tháng 4 khi khoảng biến thiên nhiệt độ ngang front (dT) là lớn nhất, dT = 28.5 – 30.50C. Front này ngăn cách giữa khối nước khá ngọt trong vịnh đang ấm lên hơn so với các khối nước của biển Đông ngang với GTF, vì vậy khoảng biến thiên mật độ bề mặt (dD) ngang GTF là tương đối lớn, dD = 19.4 – 21.4 (Yanagi & Takao, 1998). Front nhiệt vịnh Thái Lan biến mất hoàn toàn vào tháng 5. Các dữ liệu hải văn đã cho thấy rằng không xuất hiện front nhiệt trong tháng 6 trong khi vẫn tồn tại front muối trong vịnh với khoảng biến thiên độ muối bề mặt giảm xuống tới dS = 32.8 – 33.2 (Yanagi & Takao, 1998). Các dòng triều trong vịnh Thái Lan đều song song với bờ với dòng triều trung bình tháng nhỏ khoảng dưới 0.07 m/s (Wattayakorn và cs, 1998). Điều này chứng tỏ thủy triều gần như không có vai trò đáng kể trong việc hình thành và tồn tại của front trong vịnh Thái Lan. Trường gió trong vịnh không đồng nhất và biến đổi (Yanagi & Takao, 1998) và do dó ứng suất gió là âm và góp phần quan trọng trong việc hình thành GTF. Lượng nước sông Mekong đổ vào vịnh chính là nguyên nhân chính hình thành front vịnh Thái Lan (GTF). Chính vì thế, có thể dự báo được sự thay đổi đáng kể theo mùa của các đặc trưng GTF như vị trí front (do GTF không bị ảnh hưởng bởi địa hình) liên quan đến các biến trình mùa của giáng thủy và lưu lượng nước sông Mekong. Hơn nữa, hoàn lưu trong biển Đông cũng mang nước 13 sông Mekong vào vịnh cũng có thể giải thích cho sự hình thành và biến đổi của GTF. Hình 4: Front trong vịnh Thái Lan (theo Belkin I.M) [6] Tại Việt Nam, nghiên cứu front và sự hiện diện của các khối nước đã được đề tài KT03-10 mô tả, chỉ ra vị trí cho từng tháng, từng mùa làm cơ sở dữ liệu quan trọng trong dự báo phân bố nguồn lợi cá. Với kết quả mô hình hóa, mô hình KHCN- 06-02 cũng đã khẳng định sự hiện diện của các khối nước và front. Mô hình kết hợp thủy nhiệt động lực và sinh thái cũng đã giúp lý giải rõ hơn về các khu vực tập trung cá liên quan đến hoạt động của nước trồi, các front và quá trình dịch chuyển chúng do hoàn lưu chung của biển. Mô hình cũng đã chỉ ra sự hiện diện của khối nước lạnh mùa đông ven bờ Trung Quốc với nhiệt độ tương tự như khối nước lạnh ven bờ Tây vịnh Bắc Bộ. 14 Trong mùa đông, giới hạn lan truyền của khối nước nước mặt Bắc Biển Đông có thể kéo dài đến tận 5oN và uốn về hướng Đông Nam đến 112 – 1130E và dải front nhiệt kết hợp với front muối phía Đông-Nam Phú Quý có thể có vị trí trung bình theo hướng Bắc Nam. [16] Hình 5:Sơ đồ phân bố các khối nước và front trên mặt biển trong mùa đông (Đề tài KT 03-10)[15] Trong mùa hè, số lượng các khối nước vẫn giữ nguyên như trong báo cáo của đề tài KT03-10, vị trí các vùng front cũng không có sự thay đổi nhiều nhưng nguồn gốc của các front có thể được làm sáng tỏ hơn. Front tổng hợp nhiệt và muối 15 nằm giữa khối nước ấm và tương đối nhạt do kết quả tương tác với nước sông đổ ra từ vùng bờ Nam Bộ với nước lạnh và mặn khu vực nước trồi. Front tại phía Bắc Đèo Ngang cơ bản là front nhiệt do kết quả xâm nhập của nước tầng sâu đi lên tại vùng tiếp giáp với nước bị đốt nóng mạnh mùa hè do gió khô nóng tại các vùng biển nông ven bờ Hà Tĩnh, Nghệ An. Dải front phía cửa sông Hồng lại là front muối do xâm nhập nước sông ra biển. [16] Hình 6:Sơ đồ phân bố các khối nước và front trên mặt biển trong mùa hè (Đề tài KT 03-10) [15] 16 CHƢƠNG 3 HỆ THỐNG VỆ TINH VÀ DỮ LIỆU ẢNH VỆ TINH 3.1 Hệ thống các vệ tinh NOAA-AVHRR Dự án thu ảnh độ phân giải cao 4 km sử dụng thuật toán Pathfinder phiên bản 5 (Pathfinder 5) là dự án tái phân tích mới bộ số liệu độ phân giải cao (AVHRR) được phát triển bởi trường đại học Khoa học Khí quyển và Đại đương Rosenstiel thuộc đại học Miami (Rosenstiel School of Marine and Atmospheric Science – RSMAS) và Trung tâm dữ liệu hải dương học quốc gia Mỹ NOAA (NODC). Với nhiều năm thực hiện và phát triển các bộ dữ liệu Pathfinder trước đó, Trung tâm thu dữ liệu chủ động, phân phối dữ liệu hải dương học vật lý của NASA (NASA's Physical Oceanography Distributed Active Archive Center - PO.DAAC) cũng đã cùng tham gia thực hiện dự án này. Với sự phát triển đó, dự án này đã được hy vọng sẽ mở rộng hệ thống ghi dữ liệu nhiệt mặt biển (SST) toàn cầu. Chương trình Pathfinder hải dương học với ảnh độ phân giải cao (AVHRR) của NOAA/NASA đã phát triển đối với nhiệt mặt biển toàn cầu độ phân giải 9.28 km từ đầu thập kỷ 90 với kết quả là bộ dữ liệu phiên bản 4.x (V4). Mặt dù đã có những thành công nhất định, nhưng bộ dữ liệu này cũng có những hạn chế nhất định, ví dụ như độ lệch đối với một khu vực bị tác động bởi các bất đồng nhất khí quyển (atmospheric aerosols) lớn. Hơn nữa, khả năng phân tách đất khá thô là một vấn đề lớn đối với các ứng dụng khu vực ven biển và làm hạn chế các thông tin về sự không tồn tại của các băng biển trong việc sử dụng các dữ liệu vĩ độ cao. Với 3 năm phát triển, RSMAS/NODC đã sử dụng thuật toán tiền xử lý mới đã giải quyết được những hạn chế của phiên bản 4.x và đưa ra bộ dữ liệu nhiệt mặt biển chính xác hơn với độ phân giải cao hơn. Quá trình tiền xử lý dữ liệu vệ tinh sử dụng phiên bản nâng cao của thuật toán Pathfinder và các bước xử lý để đưa ra 2 ảnh nhiệt mặt biển toàn cầu trong 1 ngày và các thông số liên quan từ những năm 1985. Với độ phân giải xấp xỉ 4 km, bộ số liệu này được coi là có độ phân giải cao nhất có thể đối với dữ liệu độ phân giải cao 17 (AVHRR) toàn cầu. Không những thế, bộ dữ liệu này còn bao gồm các dữ liệu trung bình giai đoạn 5 ngày, 7 ngày, 8 ngày, trung bình tháng và trung bình năm. So với bộ dữ liệu nhiệt mặt biển Pathfinder độ phân giải 9 km, bộ dữ liệu này đã có những tiến bộ hơn đó là sự chính xác hơn, khả năng phân tách đất tốt hơn, độ phân giải cao hơn và bao gồm cả thông tin về băng biển. Bảng 1: Các vệ tinh NOAA được sử dụng trong bộ dữ liệu SST Pathfinder 5.0 STT Tên vệ tinh Thời gian của dữ liệu Pathfinder 5.0 1 NOAA-9 04/01/1985 – 07/11/1988 2 NOAA-11 08/11/1988 – 13/09/1994 3 NOAA-9 14/09/1994 – 21/01/1995 4 NOAA-14 22/01/1995 – 11/10/2000 5 NOAA-16 12/10/2000 – 31/12/2002 6 NOAA-17 01/01/2003 – 04/06/2005 7 NOAA-18 05/06/2005 – 31/12/2006 3.2 Cơ sở dữ liệu ảnh NOAA-AVHRR Cơ sở dữ liệu ảnh NOAA-AVHRR là cơ sở dữ liệu đã bắt đầu từ nhiều năm trước và đã có những nghiên cứu độ chính xác của cơ sở dữ liệu này. Với hàm hồi quy: y = 0.075*x – 3.0, Lee và cs (2005) đã nghiên cứu độ chính xác của cơ sở dữ liệu này đối với khu vực biển xung quanh Đài loan. Nghiên cứu này đã cho thấy độ chính xác của CSDL này là 0.60C [12]. Barton (1995 và Kearns & cs (2000) đã chứng minh rằng độ chính xác của CSDL này là từ 0 – 0.240C[10]. Trong khi thực hiện luận văn này, tác giả cũng đã liên hệ với TS. Jason Roberts – Đại học Duke, TS. Jason Roberts (2002) cũng đã nghiên cứu so sánh với các số liệu của các trạm phao khu vực Thái Bình Dương cũng cho độ lệch chuẩn là 0.790C. Những số liệu này có thể cho thấy được rằng, CSDL ảnh NOAA-AVHRR là tin cậy và có thể sử dụng để nghiên cứu cường độ và vị trí front trong biển Đông. 18 Luận văn sử dụng các số liệu ảnh NOAA trung bình tháng từ năm 1985 – 2006 để nghiên cứu và xác định vị trí của front cũng như tính toán tần suất xuất hiện của front. Đồng thời, luận văn sử dụng số liệu trung bình 5 ngày và trung bình tháng năm 2006 theo cả số liệu chụp ban ngày (10h địa phương – vệ tinh đi xuống) cũng như ban đêm (22h địa phương – vệ tinh đi lên) để nghiên cứu chi tiết sự biến động của cường độ front theo ngày đêm, 5 ngày và theo tháng. 19 CHƢƠNG 4 PHƢƠNG PHÁP XÁC ĐỊNH FRONT TỪ ẢNH VỆ TINH 4.1 Thuật toán xác định ranh giới các ảnh nhiệt mặt biển (SST) của Cayula & Cornillon (1992) [7] - Cần phân chia ảnh toàn cảnh thành các cửa sổ với kích thước cho trước - Mỗi cửa sổ sẽ được xử lý độc lập - Các cửa sổ cần thiết kế có phần chồng lên nhau để có thể liên kết các ranh giới xác định được trên các cửa sổ khác nhau - Kích thước cửa sổ có thể lựa chọn theo: 16x16, 32x32,… tùy theo kích thước các ranh giới cần xác định cũng như độ phân giải của ảnh, nhiễu… - Ranh giới xác định được là biên của hai khu vực ω1 và ω2 có nhiệt độ gần như đồng nhất a) Phân tích Histogram: Hình 7: Histogram của cửa sổ có ranh giới thể hiện yếu (không rõ ràng), thông qua sự phân bố các khoảng cực đại phân tách yếu Hình 8: Histogram của cửa sổ có ranh giới thể hiện mạnh (rõ ràng), thông qua sự phân tách rõ các khoảng cực trị 20 Hình 9: Các histogram tương ứng được thể hiện trong window 2 và window 1 trong hình trên Trong phép phân tích histogram, cần xác định giá trị ngưỡng để phân tách hai miền giá trị nhiệt có sự phân dị đáng kể. Trong trường hợp có hai miền nhiệt phân dị, cần sử dụng phương pháp thống kê để xác định (phân loại) chúng. Giả sử: + ω1 và ω2 là hai miền nhiệt xác định + x là điểm ảnh thuộc miền xác định X (X là toàn bộ các điểm trong cửa sổ đã chọn) + t(x) là giá trị nhiệt tại x Khi đó, mật độ xác suất cho các giá trị nhiệt trên 2 miền ω1 và ω2 được xác định theo công thức: 𝑝 𝑡 𝑥 = 𝑝[𝑡(𝑥) 𝑥 ∈ 𝜔𝑖]𝑃(𝑥 ∈ 𝜔𝑖) 2 𝑖=1 Nếu giả thiết, hai miền nhiệt thay đổi, có thể xác định được giá trị ngưỡng lựa chọn T để đánh giá sai số phân loại theo tiêu chí sai số cực tiểu Bayes, khi đó: Hàm mật độ xác suất đánh giá sai số P(error) được mô tả theo công thức: 21 𝑃 𝑒𝑟𝑟𝑜𝑟 = 𝑃 𝑡 𝑥 ≥ 𝜏 𝑥 ∈ 𝜔1]𝑃( 𝑥 ∈ 𝜔𝑖1 + 𝑃 𝑡 𝑥 < 𝜏 𝑥 ∈ 𝜔2]𝑃( 𝑥 ∈ 𝜔2 Để đơn giản, giả thiết P(error) là không đáng kể và có thể loại bỏ. Khi đó, nếu xác định được τ tối ưu, các tham số kỳ vọng và phương sai của hai miền nhiệt sẽ xác định được theo công thức: 𝜇1 𝜏 = 𝑡ℎ(𝑡)𝑡<𝜏 ℎ(𝑡)𝑡<𝜏 𝜇2 𝜏 = 𝑡ℎ(𝑡)𝑡≥𝜏 ℎ(𝑡)𝑡≥𝜏 𝜎1 2 𝜏 = (𝑡−𝜇1𝑡<𝜏 ) 2ℎ(𝑡) ℎ(𝑡)𝑡<𝜏 𝜎2 2 𝜏 = (𝑡−𝜇2𝑡≥𝜏 ) 2ℎ(𝑡) ℎ(𝑡)𝑡≥𝜏 Trong đó: t: Nhiệt độ h(t): Giá trị histogram tại t (h(t)=[0,n]) n: Tổng số giá trị nhiệt trong miền xác định Khi đó, hàm mật độ xác suất tính theo ngưỡng tối ưu τ sẽ được chọn sao cho J(τ) cực tiểu: 𝐽 𝜏 = {𝑝[𝑡 𝜇1 𝜏 , 𝑛 𝑡=0 𝜇2 𝜏 ,𝜎1 𝜏 ,𝜎2 𝜏 ]− ℎ 𝑡 } 2 Sau phân loại theo hướng τ, nếu đặt Je(τ) là phương sai trong mỗi miền và Jb(τ) là phương sai theo 2 miền đã phân loại. Jtot tổng là: Jtot = Je(τ) + Jb(τ) Trong đó: 𝐽𝑒 𝜏 = 𝑁1 𝑁1 + 𝑁2 𝑆1 𝜏 + 𝑁2 𝑁1 + 𝑁2 𝑆2 𝜏 𝑆1 𝜏 = 𝑡 − 𝜇1 𝜏 2ℎ(𝑡)𝑡<𝑟 𝑁1 𝑆2 𝜏 = 𝑡 − 𝜇2 𝜏 2ℎ(𝑡)𝑡≥𝑟 𝑁2 22 Và 𝑁1 = ℎ(𝑡)𝑡<𝑟 𝑁2 = ℎ(𝑡)𝑡≥𝑟 𝐽𝑏 𝜏 = 𝑁1𝑁2 (𝑁1 + 𝑁2)2 [𝜇1 𝜏 − 𝜇1 𝜏 ] 2 Nếu chọn τ tối ưu, thì khi đó: Jb(τ) >>Je(τ) b) Thuật toán liên kết (Cohension Algorithm) Sử dụng trong trường hợp tính đến sự phân bố không gian của dữ liệu Giả sử: + ω1’ và ω2’ là hai miền xác định + x: là điểm trên ảnh xác định + t(x) là giá trị nhiệt độ tương ứng + τopt: giá trị ngưỡng tối ưu phân chia hai miền ω1’ và ω2’ Hệ số liên kết (cohension cofficient) đối với ω1’ và ω2’ cho toàn bộ miền xác định được xác định như sau: 𝐶1 = 𝑅1 𝑇1 𝐶2 = 𝑅2 𝑇2 𝐶 = 𝑅1+𝑅2 𝑇1+𝑇2 Trong đó: T1: Tổng số điểm so sánh giữa nhiệt độ tại x các giá trị thuộc ω1’ và giá trị thuộc các miền nhiệt khác, được biểu diễn như sau: 𝑇1 = { 𝑥,𝑦 , 𝑦 ∈ 𝑁 𝑥 ∩ 𝜔1 ′ ,∀𝑥 ∈ 𝜔1 ′ } R1: Tổng số điểm so sánh giữa nhiệt độ tại x và các điểm xung quanh thuộc ω1’ 𝑅1 = { 𝑥,𝑦 , 𝑦 ∈ 𝑁 𝑥 ∩ 𝜔1 ′ ,∀𝑥 ∈ 𝜔1 ′ } Tương tự: với T2 và R2 Nếu xi,j là điểm trung tâm thì: 23 N(xi,j) = {xi,j+1,xi,j-1,xi+1,j,xi+1,j+1} Hệ số liên kết có giá trị lớn khi điểm xi,j không gần với ranh giới và các điểm xung quanh đều thuộc về một miền nhiệt Trong trường hợp lý tưởng, các miền nhiệt gần như đồng nhất, hệ số liên kết: C=0.9, C1=0.92, C2=0.92. c) Vị trí các ranh giới: Sau khi đã xác định được ranh giới theo từng cửa sổ. Bước tiếp theo là xác định và khẳng định các front. Ảnh front được thể hiện dưới dạng các pixel không có front sẽ nhận giá trị 0. Nếu giá trị nhiệt tại miền ω1’ và ω2’ được thể hiện theo hàm Ω(x), khi đó: ∀𝑥 ∈ 𝝌,Ω 𝑥 = 0 𝑛ế𝑢 𝑥 ∈ 𝜔1 ′ 1 𝑛ế𝑢 𝑥 ∈ 𝜔2 ′ Các điểm ranh giới được thể hiện dưới dạng: ∀𝑥 ∈ 𝝌 nếu tồn tại 𝑦 ∈ 𝑁(𝑥) ′ sao cho Ω(𝑥) ≠ Ω(𝑦) Thì khi đó => x là vị trí ranh giới 24 4.2 Xác định ranh giới sử dụng nhiều ảnh SST Thuật toán xác định ranh giới sử dụng nhiều ảnh thực chất là tổ hợp nhiều thuật toán gồm: + Xác định mây theo đơn ảnh + Xác định mây theo đa ảnh + Xác định ranh giới sử dụng ảnh đơn lẻ + Xác định các ranh giới Trình tự thực hiện như sau: - Sử dụng phương pháp xác định ranh giới theo đơn ảnh (Cayula & Cornillon, 1992) và áp dụng cho nhiều ảnh đơn lẻ - Liên kết các ranh giới theo từng ảnh - Tối thiểu hóa độ rộng của ranh giới (Xử lý để ranh giới chỉ tồn tại trên 1 điểm) - Kết hợp ranh giới đã xác định với ảnh ban đầu để xác định các ranh giới yếu - Thành lập mask cho các khu vực mây che phủ * Xác định các front: Đường C có độ dài phân bố trên N điểm được xác định bởi tập hợp các pixel như sau: 𝐶 = {(𝑥𝑝 , 𝑦𝑝 𝑝 ∈ 1,𝑁 } Mỗi pixel thuộc C đều bắt đầu 1 đoạn có độ dài λ. Độ dài này được lựa chọn theo đặc trưng của vùng nghiên cứu, độ phân giải ảnh và khoảng các liên kết. Đoạn Sq thuộc C được xác định trên miền 𝑞 ∈ [1,𝑁] bởi: 𝑆𝑞 = {(𝑥𝑝 ,𝑦𝑝 𝑝 ∈ 𝑞, 𝑞 + 𝜆 − 1 ∩ [1,𝑁]} Phương pháp liên kết các đoạn Sq thuộc contour của các ảnh khác nhau được thông qua trường gradient (tốt hơn sử dụng trường nhiệt) do trường nhiệt các ảnh thường khác nhau. Khi đó, hàm gradient m(grad1,grad2) được xác định như sau: 25 𝑚 𝑔𝑟𝑎𝑑1,𝑔𝑟𝑎𝑑2 = 0 𝑛ế𝑢 𝑔𝑟𝑎𝑑1.𝑔𝑟𝑎𝑑2 ≤ 0 𝑔𝑟𝑎𝑑1.𝑔𝑟𝑎𝑑2 𝑔𝑟𝑎𝑑1 2 𝑛ế𝑢 𝑔𝑟𝑎𝑑1 > 𝑔𝑟𝑎𝑑2 𝑔𝑟𝑎𝑑1.𝑔𝑟𝑎𝑑2 𝑔𝑟𝑎𝑑2 2 𝑛ế𝑢 𝑔𝑟𝑎𝑑1 ≤ 𝑔𝑟𝑎𝑑2 Kết quả sau liên kết Sq thông qua t = (tx,ty) 𝑀 𝑡 = 𝑚[𝑔𝑟𝑎𝑑𝑁 𝑥𝑝 ,𝑦𝑝 ,𝑔𝑟𝑎𝑑𝐶( (𝑥𝑝 ,𝑦𝑝 )∈𝑆𝑞 𝑥𝑝 + 𝑡𝑥 ,𝑦𝑝 + 𝑡𝑦)] Hàm grad(i,j) được xác định như sau: 𝑔𝑟𝑎𝑑 𝑖, 𝑗 = 𝑔𝑟𝑎𝑑𝑥(𝑖, 𝑗) 𝑔𝑟𝑎𝑑 𝑦 (𝑖, 𝑗) = 𝜏𝑖+1,𝑗−𝜏𝑖−1,𝑗 𝜏𝑖 ,𝑗+1 − 𝜏𝑖 ,𝑗−1 4.3 Tính toán Gradient Tính toán Gradient sử dụng các ma trận tính gradient có công thức tổng quát như sau: GX={GXdimension, direction (+) 0 GXdimension, direction (-)} GY={GYdimension, direction (+) 0 GYdimension, direction (-)} với dimension là kích thước ma trận tính gradient, ví dụ 3x3, 5x5... direction là hướng tính gradient, có thể là hai hay nhiều hướng Luận văn sử dụng vec tơ tính gradient được tính toán theo toán tử Prewitt bao gồm 2 ma trận kích thước 3x3, tính theo hai hướng NS và EW như sau: 𝐺𝑋 = −1 0 1 −1 0 1 −1 0 1 𝐺𝑌 = 1 1 1 0 0 0 −1 −1 −1 Các ảnh được tính toán với hai ma trận này, kết quả cho ta các giá trị theo các hướng X,Y. 26 Giả sử ma trận của ảnh ban đầu là A ta có: Gx=GX*A Gy=GY*A Tại mỗi điểm ảnh, cường độ gradient được tính theo công thức: 𝐺𝑀 = 𝐺𝑥2 + 𝐺𝑦2 [3] Toán tử Prewitt là một dạng ma trận đơn giản, được sử dụng rộng rãi trong tính toán gradient Phương pháp gradient cho cái nhìn tổng quan về cường độ của front trong khu vực biển rộng lớn. Đây là phương pháp phổ biến hiện nay để xác định các ngưỡng phổ biến nhất của front từ đó tùy theo mục đích nghiên cứu front có thể định ra các bản đồ front theo ngưỡng. 27 4.4 Quy trình tính toán front nhiệt mặt biển với phƣơng pháp của Cayula & Cornillon (1992) Hình 10: Sơ đồ khối của quy trình tính toán front theo phương pháp của Cayula & Cornillon(1992) 28 Các quy trình tính toán trên được thực hiện với sự hỗ trợ của phần mềm MGET (Marine Geology and Environment Tools). Đây là một bộ gồm hơn 150 công cụ hữu ích với nhiều chức năng khác nhau. Những công cụ tiêu biểu của MGET đó là chuyển đổi các dữ liệu hải dương học sang format của phần mềm ArcGIS, xác định các front bằng ảnh nhiệt bề mặt hay các chương trình thống kê... Đây là một gói các chương trình nguồn mở được viết bằng ngôn ngữ Python được trường Đại học Duke, Mỹ phát triển cho các nhà phân tích GIS và nghiên cứu biển. Bộ chương trình này có thể download miễn phí dễ dàng từ website của nhóm phát triển ( [11]. Mặc dù hiện nay, MGET được viết độc lập với nền hoạt động và chỉ hỗ trợ trên nền Microsoft Windows XP và Vista nhưng nhiều công cụ vẫn yêu cầu nền ArcGIS Desktop. MGET hoạt động trên nền ArcGIS như là một công cụ của ArcToolbox [11]. Ta thực hiện theo các bước như sau: - Bước 1: Chuyển đổi các file dạng HDF sang dạng ArcGIS Raster - Bước 2: Tính toán các giá trị nhiệt mặt biển từ các giá trị hồng ngoại - Bước 3: Loại bỏ mây - Bước 4: Tính toán phân bố của front nhiệt mặt biển - Bước 5: Chuyển đổi các file dạng này sang dạng ASCII - Bước 6: Tính toán tần suất xuất hiện và chuyển từ ASCII sang dạng ArcGIS Raster để hiển thị số liệu Bước 1: Chuyển đổi các file dạng HDF sang dạng ArcGIS Raster Sau khi download MGET về và cài đặt (Yêu cầu trước đó phải cài đặt ArcGIS Desktop – Chú ý chỉ có từ version 9.3 mới hoạt động được trên nền Window Vista) chúng ta cần download các ảnh hồng ngoại nhiệt lưu dưới dạng *.HDF từ trang web của NOAA. 29 Luận văn này trình bày kỹ thuật tính toán với một loạt các ảnh để đơn giản hơn các bước thực hiện khi làm với nhiều ảnh. Trước tiên chúng ta cần sử dụng Microsoft Excel tạo một bảng gồm 2 cột như sau: datain dataout E:\ThesisYear8506\198505.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\sd198505 E:\ThesisYear8506\198605.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\sd198605 E:\ThesisYear8506\198705.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\sd198705 E:\ThesisYear8506\198805.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\sd198805 E:\ThesisYear8506\198905.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\sd198905 E:\ThesisYear8506\199005.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\sd199005 Trong đó, cột thứ nhất là cột các file dữ liệu ảnh vệ tinh thu được, cột này cũng có thể là các file quality chứa các thông tin về mây. Cột thứ 2 là các Raster thu được. Tạo các bảng chứa các Raster mang thông tin về mây: datain dataout E:\ThesisYear8506\198505.m04m3pfv50-qual.hdf E:\ThesisYear8506\May\qld198505 E:\ThesisYear8506\198605.m04m3pfv50-qual.hdf E:\ThesisYear8506\May\qld198605 E:\ThesisYear8506\198705.m04m3pfv50-qual.hdf E:\ThesisYear8506\May\qld198705 E:\ThesisYear8506\198805.m04m3pfv50-qual.hdf E:\ThesisYear8506\May\qld198805 E:\ThesisYear8506\198905.m04m3pfv50-qual.hdf E:\ThesisYear8506\May\qld198905 E:\ThesisYear8506\199005.m04m3pfv50-qual.hdf E:\ThesisYear8506\May\qld199005 E:\ThesisYear8506\199105.m04m3pfv50-qual.hdf E:\ThesisYear8506\May\qld199105 E:\ThesisYear8506\199205.m04m3pfv50-qual.hdf E:\ThesisYear8506\May\qld199205 E:\ThesisYear8506\199305.m04m3pfv50-qual.hdf E:\ThesisYear8506\May\qld199305 và các Raster mang giá trị thực của nhiệt: datain dataout E:\ThesisYear8506\198505.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\Td198505 E:\ThesisYear8506\198605.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\Td198605 E:\ThesisYear8506\198705.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\Td198705 E:\ThesisYear8506\198805.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\Td198805 E:\ThesisYear8506\198905.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\Td198905 E:\ThesisYear8506\199005.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\Td199005 E:\ThesisYear8506\199105.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\Td199105 E:\ThesisYear8506\199205.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\Td199205 E:\ThesisYear8506\199305.s04m3pfv50-sst-16b.hdf E:\ThesisYear8506\May\Td199305 30 Chúng ta sử dụng công cụ chuyển đổi từ HDF -> Raster Hình 11: Cửa sổ giao diện của MGET với công cụ chuyển từ HDF sang ArcGIS Raster 31 Công cụ chuyển đổi các file dạng HDF sang Raster có các thông số: Hình 12: Cửa sổ các thông số khi chuyển từ HDF -> ArcGIS Raster Với cửa sổ này, chúng ta cần nhập các thông số để xác định dạng vào ra của công cụ chuyển đổi này: - Tên file chứa các bảng: Cần nhập tên bảng vừa tạo được ở trên - Nhập tên các trường: Chọn đầu vào, ra (input/output) tương ứng với tên hai cột vừa tạo trong bảng - Tên của cơ sở dữ liệu HDF: Với nhiệt độ là sst - Nhập vị trí góc dưới cùng bên trái của CSDL - Nhập độ phân giải của dữ liệu - Nhập hệ quy chiếu: Các Raster trong luận văn này có hệ quy chiếu là WGS- 84 32 - Nhập tọa độ khu vực nghiên cứu: Kinh độ: 99 E – 121 E Vĩ độ: 1 N – 25 N Thực hiện tương tự với một chuỗi các file quality để loại bỏ mây. Đây là các file chứa thông tin về mây với các giá trị qual = 1,..,7. Bước 2: Tính toán các giá trị thực của nhiệt mặt biển Riêng với tính toán các giá trị thực của nhiệt mặt biển chúng ta cần nhập hàm hồi quy như sau: Hình 13: Cửa sổ các thông số khi tính toán giá trị thực nhiệt mặt biển 33 Với các bước thực hiện như trên chúng ta đã chuyển đổi thành công các dữ liệu từ dạng file HDF sang dạng ArcGIS Raster và tính toán được các giá trị thực nhiệt mặt biển. Bước 3: Loại bỏ mây (SetNull) Như đã đề cập ở trên, trong quá trình loại bỏ mây chúng ta cần sử dụng 2 dạng số liệu, đó là số liệu hồng ngoại nhiệt (hoặc nhiệt) và các thông tin về mây của mỗi file. Để coi một pixel cho trước là giá trị chứa mây hay đất (chuyển thành dạng không có dữ liệu cần quan tâm NoData) chúng ta cần cho các giá trị pixel của file quality nhỏ hơn một giá trị nào đó, các giá trị này từ 0 (chất lượng thấp) đến 7 (chất lượng cao). Để giảm thiểu các giá trị front bị tính sai do ảnh hưởng của mây, luận văn đã chọn giá trị tốt nhất là qual < 7. Tiến sĩ Jason Roberts – người đứng đầu nhóm nghiên cứu và tạo lập nên công cụ MGET, trường Đại học Duke Mỹ đã nghiên cứu so sánh các giá trị AVHRR SST với các giá trị SST đo được bởi ngư dân ở phía Đông Thái Bình Dương cận nhiệt đới, thu được các giá trị như sau: Bảng 2: Sai khác giữa các giá trị AVHRR SST và các giá trị in-situ SST đối với các ngưỡng quality khác nhau (tổng giá trị điểm ảnh = 858) Các giá trị quality n Độ lệch SST trung bình (°C) Độ lệch tiêu chuẩn (°C) 0 657 11.21 10.83 1 21 7.09 6.42 2 34 7.36 4.17 3 12 0.51 1.77 4 5 0.38 1.31 5 26 0.34 0.74 6 18 -0.31 0.54 7 85 0.02 0.79 Bảng trên là sự so sánh giữa các giá trị in-situ SST và ảnh SST hàng ngày. Với chất lượng tốt nhất là qual < 7, chúng ta chỉ còn 85 điểm ảnh (pixel) không có mây chiếm 10% với độ lệch SST trung bình là 0.02oC và độ lệch tiêu chuẩn là 0.79 oC. Đối với các số liệu trung bình tháng sử dụng trong luận văn nên với qual < 34 7 đã chứa rất ít mây có thể sử dụng để tính toán rất tốt. Với số liệu hồng ngoại nhiệt trung bình 5 ngày, số điểm ảnh không mây thu được sẽ nhỏ hơn nhưng đáp ứng tốt việc nghiên cứu front nhiệt mặt biển. Để loại bỏ mây, sử dụng hàm SetNull trong công cụ Spatial Analyst, cũng tương tự như chuyển các file từ HDF sang ArcGIS Raster cần tạo các bảng: E:\ThesisYear8506\12\qld198512 E:\ThesisYear8506\12\Td198512 E:\ThesisYear8506\12\fnd198512 "VALUE" <7 E:\ThesisYear8506\12\qld198612 E:\ThesisYear8506\12\Td198612 E:\ThesisYear8506\12\fnd198612 "VALUE" <7 E:\ThesisYear8506\12\qld198712 E:\ThesisYear8506\12\Td198712 E:\ThesisYear8506\12\fnd198712 "VALUE" <7 E:\ThesisYear8506\12\qld198812 E:\ThesisYear8506\12\Td198812 E:\ThesisYear8506\12\fnd198812 "VALUE" <7 E:\ThesisYear8506\12\qld198912 E:\ThesisYear8506\12\Td198912 E:\ThesisYear8506\12\fnd198912 "VALUE" <7 E:\ThesisYear8506\12\qld199012 E:\ThesisYear8506\12\Td199012 E:\ThesisYear8506\12\fnd199012 "VALUE" <7 Trong đó cột SST Raster là cột chứa các file quality, cột thứ 2 là cột chứa các Raster hồng ngoại nhiệt (nhiệt), cột thứ 3 là cột chứa các Raster kết quả và cột cuối cùng là cột điều kiện để loại bỏ mây. Tiến hành như sau: Hình 14: Cửa sổ giao diện khi loại bỏ mây 35 Hình 15: Giao diện cửa sổ các thông số đầu vào khi loại bỏ mây Sau đó chỉ cần chọn OK là chúng ta đã setNull được hàng loạt các file dữ liệu hồng ngoại vệ tinh hay dữ liệu các giá trị thực của nhiệt độ. Áp dụng các bước thực hiện trên ta sẽ có thể dễ dàng thu được một CSDL nhiệt độ mặt biển phong phú. Bước 4: Tính toán phân bố front nhiệt mặt biển Sau bước 3 đã thu được một CSDL các Raster nhiệt mặt biển và hồng ngoại nhiệt không chứa mây. Tính toán front trên các Raster hồng ngoại nhiệt đó ta sẽ thu được phân bố front. Để tính toán với một loạt các Raster chúng ta cũng cần tạo các bảng chứa các cột Raster đầu vào và kết quả. Tuy nhiên, một điểm khác biệt đó là chúng ta cần tạo một bảng dạng DBF. Để tạo nên bảng này cần sử dụng chương trình Find Rasters trong bộ công cụ MGET, các thông số cần để chạy chương trình như sau: 36 - Thư mục để tìm Raster: Chọn thư mục chứa các Raster hồng ngoại nhiệt - Thư mục chứa kết quả - Tên file DBF kết quả Hình 16: Giao diện cửa sổ chứa các thông số tìm kiếm Raster Một file DBF có tên là AUGtable.DBF sẽ được tạo ra với một cột chứa tất cả các Raster mà công cụ tìm kiếm Raster (Find Rasters) đã tìm được trong thư mục chứa các Raster hồng ngoại. Muốn tạo cột chứa các Raster kết quả cần thêm 1 trường mới có tên là Outfront với công cụ AddField trong bộ công cụ của ArcGIS. Tạo trường này với dạng dữ liệu là TEXT, độ dài 250. 37 Hình 17: Giao diện cửa sổ thêm trường vào một bảng cho trước Để gán các tên kết quả cho trường Outfront vừa tạo được sử dụng công cụ tính toán các trường (Calculate Field) với việc sử dụng một lệnh VBScript: " E:\ThesisYear8506\Aug\f" + Right([Image], 14) Toàn bộ lệnh này có ý nghĩa tạo một loạt các file Raster trong trường cho trước là Outfront với địa chỉ thư mục là E:\ThesisYear8506\Aug\ sau đó lấy tên file hồng ngoại nhiệt 14 ký tự tính từ bên phải của chuỗi và thêm vào ký tự f ở đầu tên. 38 Hình 18: Giao diện cửa sổ tính toán với các trường 39 Hình 19: Giao diện công cụ xác định front với một bảng các Raster Sau khi đã có một bảng các Raster đầu vào/ra của chương trình tính front, chạy trực tiếp công cụ tính toán front với một loạt các Raster (Cayula-Cornillon Fronts in ArcGIS Rasters in Table) với các thông số như sau: - Kích thước cửa sổ lọc: 3 - Kích thước cửa sổ histogram: 16 - Số lượng pixel được cửa sổ lọc lặp đi lặp lại: 1 - Tỷ lệ tối thiểu của pixel với dữ liệu: 0.65 . Đây chính là tỷ lệ của không phải là NoData (chứa đất/mây) với dữ liệu. Ví dụ nếu tỷ lệ này là 0.75 và kích thước cửa sổ histogram là 20x20 thì khi đó phải có ít nhất 300 pixel trong 400 pixel phải không chứa giá trị Nodata 40 - Độ lệch giá trị hồng ngoại nhiệt tối thiểu để xác định front: 6.666667. Đây chính là cơ sở để xác định front theo một ngưỡng nào đó. Do cơ sở dữ liệu của chúng ta là dạng dữ liệu 16 bit do vậy mỗi một giá trị sai khác của giá trị hồng ngoại nhiệt sẽ tương ứng với 0.0750C. Ta xác định ngưỡng của chúng ta là 0.5 0 C nên giá trị này sẽ tính bằng: 0.5/0.075 = 6.666667 - Các thông số còn lại có thể chọn theo các giá trị của lý thuyết của Cayula & Cornillon (1992) Hình 20: Giao diện cửa sổ các thông số tính toán front Sau bước tính toán này chúng ta đã thu được một cơ sở dữ liệu bao gồm các Raster chứa các thông tin về phân bố front. Sau đó, chuyển các Raster này sang dạng ASCII bằng công cụ chuyển từ Raster sang ASCII sẵn có trong ArcGIS rồi lập trình một chương trình FORTRAN tính toán tần suất xuất hiện của front trong 22 năm từ 1985 – 2006 rồi sử dụng công cụ chuyển từ ASCII sang ArcGIS Raster chúng ta thu được các Raster chứa thông tin về tần suất xuất hiện của front nhiệt mặt biển. Tần suất xuất hiện front chính là 41 tỷ số giữa số lần xuất hiện front tại một pixel cho trước với số lần pixel cho trước đó không chứa mây, theo công thức: 𝐹 = 𝑁 𝐶 Trong đó: N: Số lần xuất hiện front tại một pixel cho trước C: Số lần pixel cho trước đó không chứa mây F: Tần suất xuất hiện (%) 42 CHƢƠNG 5: XÁC ĐỊNH FRONT NHIỆT MẶT BIỂN KHU VỰC BIỂN ĐÔNG 5.1 Tiền xử lý ảnh vệ tinh Các ảnh sau khi được download về được ứng dụng trực tiếp tính toán front theo phương pháp của Cayula & Cornillon. Tuy nhiên, với phương pháp gradient, để đơn giản hơn trong việc tính toán và xác định ngưỡng biến đổi của nhiệt mặt biển cần tiến hành thực hiện tiền xử lý ảnh. Trước tiên là lựa chọn ảnh: Hình 21: Ảnh hồng ngoại nhiệt tháng 01 năm 1985 43 Hình 22: Ảnh hồng ngoại nhiệt tháng 04 năm 2006 44 Hình 23: Ảnh hồng ngoại nhiệt tháng 08 năm 2006 Hình 24: Histogram ảnh hồng ngoại nhiệt 01/1985 Hình 25: Histogram ảnh hồng ngoại nhiệt 01/1985 Hình 26: Histogram ảnh hồng ngoại nhiệt 01/1985 45 Bảng 3: Bảng thống kê cấp độ xám 01/1985 04/2006 08/2006 Cực tiểu 0 0 0 Cực đại 255 255 255 Đỉnh cực trị lớn 202, 215, 225, 245 & 254 243 246 Độ lệch chuẩn 112.11 110.92 119.41 Giá trị trung bình 132.10 132.60 135.86 Trung vị 200 200 235 Với hai ảnh hồng ngoại nhiệt như tháng 01 năm 1985 và tháng 4 năm 2006, các ảnh đều bao gồm đầy đủ các cấp độ xám từ 0 – 255, các cực đại thể hiện được ranh giới các địa vật khác nhau.Với tháng 01 năm 1985, cấp độ xám từ khoảng 0 – 30 chiếm 40,27% là khoảng giá trị cấp độ xám của các khu vực bao gồm đất và mây. Khoảng giá trị từ 162 – 255 chiếm 54,99% chứa các thông tin nhiệt mặt biển. Một điểm rất đáng chú ý đó là trong dải cấp độ xám này lại chứa tới 5 đỉnh (peak) cấp độ xám khác nhau tại các giá trị: 202, 215, 225, 245 & 254. Điều này đã chứng tỏ rằng ảnh hồng ngoại nhiệt này chứa nhiều ngưỡng biến đổi nhiệt khác nhau. Điều này cho thấy ảnh hồng ngoại nhiệt tháng 01 năm 1985 sẽ chứa rất nhiều các cấp độ khác nhau của giá trị gradient nhiệt mặt biển. Với ảnh hồng ngoại nhiệt tháng 4 và tháng 8 năm 2006, các giá trị cấp độ xám biến đổi ít hơn, phân bố cấp độ xám cũng chia thành 2 khoảng như tháng 01 năm 1985 nhưng trong khoảng giá trị cấp độ xám của nhiệt độ có ít đỉnh cực trị hơn hoặc các đỉnh phân bố gần nhau và ít tập trung hơn. Điều này cho thấy với ảnh hồng ngoại nhiệt tháng 4 và tháng 8 có ít sự biến đổi nhiệt hơn, ít có ngưỡng biến đổi rõ rệt hơn đồng nghĩa với việc các dải front sẽ ít xuất hiện và không liên tục như tháng 01 năm 1985. 46 5.2 Xác định biến thiên gradient nhiệt mặt biển bằng phƣơng pháp tính toán gradient Tính toán gradient với phương pháp tính toán gradient cho các ảnh hồng ngoại năm 2006 đối với cả các số liệu trung bình tháng và trung bình 5 ngày. Với các số liệu trung bình 5 ngày và trung bình tháng năm 2006, qua tính toán gradient đã xác định được ngưỡng biến thiên các giá trị gradient nhiệt mặt biển chủ yếu nằm trong khoảng 0.125 – 10C/4km. Hình 27: Phân bố gradient tháng 2 năm 2006 47 Hình 28: Phân bố gradient nhiệt độ tháng 4 năm 2006 48 Hình 29: Phân bố gradient nhiệt độ tháng 8 năm 2006 49 Hình 30: Phân bố gradient tháng 11 năm 2006 Thông qua việc tính toán giá trị gradient các tháng năm 2006 cho thấy: Các front xuất hiện không đồng nhất trong biển Đông, gradient của các front có giá trị mạnh, yếu nằm xen kẽ nhau. Tại các khu vực có chênh lệch nhiệt độ lớn, giá trị gradient xuất hiện lớn hơn rõ rệt. Nghiên cứu chi tiết hơn về sự biến động của front theo trung bình 5 ngày chúng ta thu được: Trong thời gian 5 ngày, có thể thấy được rằng các dải front có sự thay đổi rõ rệt, chúng có thể di chuyển, biến mất hay hình thành front mới. Điều này rất quan trọng trong việc nghiên cứu các trường hải dương học khác. Trong hình 30,31 các dải front ở phía Nam đảo Hải Nam đã thể hiện rất rõ hiện tượng này. 50 Hình 31: Phân bố cường độ gradient trung bình 5 ngày từ 11-15/11/2006 (Trung bình ban đêm) Hình 32: Phân bố cường độ gradient trung bình 5 ngày từ 11-15/11/2006 (Trung bình ban ngày) 51 5.3 Nghiên cứu xác định vị trí, tần suất xuất hiện của front Sử dụng ảnh trung bình tháng từ năm 1985 – 2006, tính toán theo phương pháp xác định front của Cayula và Cornillon (1992), sau đó tính toán tần suất xuất hiện của front. Phương pháp của Cayula và Cornillon là phương pháp tương đối đơn giản và hiệu quả. Phương pháp của Cayula và Cornillon định ra một tuyến các giá trị bằng 1 nếu pixel đó chứa thông tin về front và bằng 0 nếu không chứa các thông tin về front. Tuy nhiên, phương pháp này cũng bị hạn chế do không định rõ chính xác cường độ front của một front nào đó mà xác định được vị trí của front với một cường độ đã biết trước. 52 Hình 33: Trường nhiệt mặt biển (SST) tháng 01 năm 1985 High : 30.75 Low : 13.275 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N High : 30.75 Low : 13.275 53 Hình 34: Trường nhiệt mặt biển (SST) tháng 8 năm 2006 Với phương trình hồi quy cho trước của ảnh NOAA-AVHRR: y = 0.075 * x – 3.0 Trong đó: - x: Giá trị cấp độ xám - y: Giá trị nhiệt độ mặt biển Có thể áp dụng trực tiếp hàm này với tất cả các ảnh sẽ thu được một cơ sở dữ liệu (CSDL) nhiệt tương đối phong phú. Hình 33 & 34 chính là hai ví dụ phân bố 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N High : 32.325 Low : 26.55 54 nhiệt độ khu vực nghiên cứu. Với các kết quả đó, ta thấy được xu thế biến đổi của nhiệt mặt biển. Tháng 01 năm 1985 trường nhiệt thể hiện xu thế của trường nhiệt vào mùa đông: Hình thành một dải lưỡi lạnh từ phía Bắc xâm nhập sâu xuống phía Nam. Trường nhiệt tháng 08 năm 2006 thể hiện rõ xu thế của trường nhiệt mùa hè: Xuất hiện một khu vực nước trồi hoạt động ở khu vực Nam Trung Bộ có xu thế lan dần ra phía giữa Biển Đông. 55 Hình 35:Phân bố front nhiệt mặt biển và phân bố nhiệt độ tại các mặt cắt tháng 01 năm 1985 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 25 26 27 28 29 1 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 1 0 1 1 1 1 1 2 1 1 3 1 0C SST 24 25 26 27 1 9 1 7 2 5 3 3 4 1 4 9 5 7 6 5 7 3 8 1 8 9 9 7 0C SST 56 Hình 36:Phân bố front nhiệt mặt biển và phân bố nhiệt độ tại các mặt cắt tháng 01 năm 1985 Phân bố front nhiệt mặt biển (~0.5oC/4km) cũng thể hiện rõ xu thế biến đổi của trường nhiệt mặt biển. Tháng 01 năm 1985 các front nhiệt mặt biển đều có xu hướng ảnh hưởng của lưỡi lạnh nhiệt độ. Các dải front đều chủ yếu có hướng Tây Bắc – Đông Nam. Tại cửa vịnh Thái Lan hình thành một dải front lớn chắn ngang cửa vịnh. Front này hình thành là do nước lạnh từ trên phía Bắc đi xuống gặp nước nóng hơn trong vịnh Thái Lan tạo thành. Dải front này cũng đã được công bố trong các bài báo của Belkin I.M [2,6]. Trong hình 35, để xác định rõ hơn xu thế biến đổi 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 57 của trường nhiệt qua các front, tác giả đã cắt 2 mặt cắt để nghiên cứu chi tiết hơn về xu hướng biến đổi này. Tại cửa vịnh Thái Lan, nhiệt độ cao hơn ở trong vịnh và giảm dần ra cửa. Tại khu vực Nam Trung Bộ thấy rõ được xu thế biến đổi của lưỡi nước lạnh, nhiệt độ thấp ở khu vực lưỡi nước lạnh xâm nhập và tăng dần ra giữa biển. Phân bố front nhiệt mặt biển tháng 8 năm 2006 tương đối phong phú. Xuất hiện nhiều dải front ở khu vực Nam Trung Bộ hơn do ảnh hưởng của hiện tượng nước trồi. Tại đỉnh vịnh Thái Lan cũng hình thành một front ngăn cách trong vịnh. Điều này là do vào thời gian này trong năm đỉnh vịnh Thái Lan có một lượng lớn nước ngọt từ sông Chao Phraya đổ vào. Tuy nhiên theo như tính toán giá trị gradient (hình 29), gradient của front nhiệt mặt biển tháng 8 yếu hơn nhiều so với gradient các tháng trong mùa đông. Vào mùa hè, sự xuất hiện và mất đi của front nhiệt mặt biển khu vực Nam Trung Bộ thường diễn ra trong thời gian khoảng từ 1 tuần cho tới hơn 10 ngày. Trên hình 37, đường màu đỏ là front từ ngày 10-14/7/2006 đã phát triển ra xa hơn so với front từ này 05-09/7/2006 (đường màu đen) (hình 37 a). Tuy nhiên, những ngày tiếp theo dải front này đã không còn phát triển nữa mà đã biến mất, thay vào đó là một dải front chạy sát bờ hơn (hình 37 b). 58 (a) (b) Hình 37: Phân bố front nhiệt mặt biển từ 05-09/7/2006 (đường màu đen) và 10-14/7/2006 (đường màu đỏ) (hình a) từ 10-14/7/2006 (đường màu đỏ) và 15-19/7/2006 (đường màu xanh đậm) (hình b) Nghiên cứu sự biến động của front theo các tháng trong năm với các bản đồ phân bố tần suất xuất hiện front từ 1985 – 2006. Tần suất xuất hiện front chính là tỷ số giữa số lần xuất hiện front tại một pixel cho trước với số lần pixel cho trước đó không chứa mây, theo công thức: 𝐹 = 𝑁 𝐶 Trong đó: N: Số lần xuất hiện front tại một pixel cho trước C: Số lần pixel cho trước đó không chứa mây F: Tần suất xuất hiện (%) 114°0'0"E 114°0'0"E 112°0'0"E 112°0'0"E 110°0'0"E 110°0'0"E 108°0'0"E 108°0'0"E 106°0'0"E 106°0'0"E 14°0'0"N 14°0'0"N 12°0'0"N 12°0'0"N 10°0'0"N 10°0'0"N 8°0'0"N 8°0'0"N 6°0'0"N 6°0'0"N 114°0'0"E 114°0'0"E 112°0'0"E 112°0'0"E 110°0'0"E 110°0'0"E 108°0'0"E 108°0'0"E 106°0'0"E 106°0'0"E 14°0'0"N 14°0'0"N 12°0'0"N 12°0'0"N 10°0'0"N 10°0'0"N 8°0'0"N 8°0'0"N 6°0'0"N 6°0'0"N 59 Hình 38: Phân bố tần suất xuất hiện của front nhiệt mặt biển trong tháng 1 (từ 1985 – 2006) 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 0 0 - 0.05 0.05 - 0.1 0.1 - 0.15 0.15 - 0.2 0.2 - 0.25 0.25 - 0.3 0.3 - 0.35 0.35 - 0.4 0.4 - 0.45 0.45 - 0.5 > 0.5 µ 60 Hình 39: Phân bố tần suất xuất hiện của front nhiệt mặt biển trong tháng 2 (từ 1985 – 2006) 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 0 0 - 0.05 0.05 - 0.1 0.1 - 0.15 0.15 - 0.2 0.2 - 0.25 0.25 - 0.3 0.3 - 0.35 0.35 - 0.4 0.4 - 0.45 0.45 - 0.5 > 0.5 µ 61 Hình 40: Phân bố tần suất xuất hiện của front nhiệt mặt biển trong tháng 3 (từ 1985 – 2006) 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 0 0 - 0.05 0.05 - 0.1 0.1 - 0.15 0.15 - 0.2 0.2 - 0.25 0.25 - 0.3 0.3 - 0.35 0.35 - 0.4 0.4 - 0.45 0.45 - 0.5 > 0.5 µ 62 Hình 41: Phân bố tần suất xuất hiện của front nhiệt mặt biển trong tháng 4 (từ 1985 – 2006) 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 0 0 - 0.05 0.05 - 0.1 0.1 - 0.15 0.15 - 0.2 0.2 - 0.25 0.25 - 0.3 0.3 - 0.35 0.35 - 0.4 0.4 - 0.45 0.45 - 0.5 > 0.5 µ 63 Hình 42: Phân bố tần suất xuất hiện của front nhiệt mặt biển trong tháng 5 (từ 1985 – 2006) 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 0 0 - 0.05 0.05 - 0.1 0.1 - 0.15 0.15 - 0.2 0.2 - 0.25 0.25 - 0.3 0.3 - 0.35 0.35 - 0.4 0.4 - 0.45 0.45 - 0.5 > 0.5 µ 64 Hình 43: Phân bố tần suất xuất hiện của front nhiệt mặt biển trong tháng 6 (từ 1985 – 2006) 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 0 0 - 0.05 0.05 - 0.1 0.1 - 0.15 0.15 - 0.2 0.2 - 0.25 0.25 - 0.3 0.3 - 0.35 0.35 - 0.4 0.4 - 0.45 0.45 - 0.5 > 0.5 µ 65 Hình 44: Phân bố tần suất xuất hiện của front nhiệt mặt biển trong tháng 7 (từ 1985 – 2006) 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 0 0 - 0.05 0.05 - 0.1 0.1 - 0.15 0.15 - 0.2 0.2 - 0.25 0.25 - 0.3 0.3 - 0.35 0.35 - 0.4 0.4 - 0.45 0.45 - 0.5 > 0.5 µ 66 Hình 45: Phân bố tần suất xuất hiện của front nhiệt mặt biển trong tháng 8 (từ 1985 – 2006) 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 0 0 - 0.05 0.05 - 0.1 0.1 - 0.15 0.15 - 0.2 0.2 - 0.25 0.25 - 0.3 0.3 - 0.35 0.35 - 0.4 0.4 - 0.45 0.45 - 0.5 > 0.5 µ 67 Hình 46: Phân bố tần suất xuất hiện của front nhiệt mặt biển trong tháng 9 (từ 1985 – 2006) 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 0 0 - 0.05 0.05 - 0.1 0.1 - 0.15 0.15 - 0.2 0.2 - 0.25 0.25 - 0.3 0.3 - 0.35 0.35 - 0.4 0.4 - 0.45 0.45 - 0.5 > 0.5 µ 68 Hình 47: Phân bố tần suất xuất hiện của front nhiệt mặt biển trong tháng 10 (từ 1985 – 2006) 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 0 0 - 0.05 0.05 - 0.1 0.1 - 0.15 0.15 - 0.2 0.2 - 0.25 0.25 - 0.3 0.3 - 0.35 0.35 - 0.4 0.4 - 0.45 0.45 - 0.5 > 0.5 µ 69 Hình 48: Phân bố tần suất xuất hiện của front nhiệt mặt biển trong tháng 11 (từ 1985 – 2006) 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 0 0 - 0.05 0.05 - 0.1 0.1 - 0.15 0.15 - 0.2 0.2 - 0.25 0.25 - 0.3 0.3 - 0.35 0.35 - 0.4 0.4 - 0.45 0.45 - 0.5 > 0.5 µ 70 Hình 49: Phân bố tần suất xuất hiện của front nhiệt mặt biển trong tháng 12 (từ 1985 – 2006) Vào mùa Đông, các front chủ yếu liên quan đến hiện tượng lưỡi nước lạnh xâm nhập sâu xuống phía Nam, hình thành các khu vực front chủ yếu: Tây Philippin, eo Đài Loan, Phía Đông đảo Hải Nam, Vịnh Bắc Bộ, Nam Trung Bộ và - Cửa vịnh Thái Lan. Vào mùa hè: Front xuất hiện chính ở: Eo Đài loan, quanh đảo Hải Nam và Nam Trung Bộ. Tại Vịnh Bắc Bộ, vào mùa đông front thường xuất hiện thành một khu vực riêng biệt ở giữa vịnh do tồn tại một khu vực nước ấm ở giữa vịnh. Đặc điểm này 120°0'0"E 120°0'0"E 115°0'0"E 115°0'0"E 110°0'0"E 110°0'0"E 105°0'0"E 105°0'0"E 100°0'0"E 100°0'0"E 25°0'0"N 20°0'0"N 20°0'0"N 15°0'0"N 15°0'0"N 10°0'0"N 10°0'0"N 5°0'0"N 5°0'0"N 0 0 - 0.05 0.05 - 0.1 0.1 - 0.15 0.15 - 0.2 0.2 - 0.25 0.25 - 0.3 0.3 - 0.35 0.35 - 0.4 0.4 - 0.45 0.45 - 0.5 > 0.5 µ 71 khác biệt với những nghiên cứu của Belkin [2,6] và khá tương đồng với kết quả nghiên cứu của Wang vào mùa Đông [10], so với những kết quả của đề tài KT03-10 kết quả này chỉ ra rõ hơn khu vực có tần suất xuất hiện front là một khu vực giữa vịnh, kéo dài hơn so với sơ đồ phân bố front của đề tài KT03-10. Tuy nhiên, luận văn này đưa ra phân bố front theo một khu vực rộng lớn hơn so với Wang và chi tiết hơn so với Belkin. Đây là những đặc điểm cần tiếp tục được nghiên cứu sâu hơn nữa. Với số liệu trung bình tháng, các front do hiện tượng nước trồi không thể hiện rõ, tuy nhiên với các số liệu trung bình 5 ngày đã cho thấy rõ biến đổi của các front nhiệt quanh khu vực này. 72 KẾT LUẬN VÀ KIẾN NGHỊ 1. Luận văn đã nghiên cứu với một CSDL lớn về các trường nhiệt mặt biển, đã xác định được các phương pháp nghiên cứu xác định front cũng như thu thập, đưa ra những khái niệm cơ bản của front. Luận văn đã thu được một cơ sở dữ liệu lớn về trường nhiệt mặt biển cũng như phân bố của front trong biển Đông. 2. Front trong biển Đông chủ yếu có giá trị gradient nằm trong khoảng 0.125 – 1 0 C/4km, mỗi front lại có giá trị gradient khác nhau, các giá trị này nằm xen kẽ nhau. Cường độ gradient có giá trị lớn hơn vào mùa đông so với mùa hè. Đây là những kết luận rất quan trọng trong nghiên cứu front trong biển Đông. 3. Luận văn cũng đã tính toán được phân bố tần suất xuất hiện của front trong tất cả các tháng với số liệu trung bình tháng từ năm 1985 – 2006 và đưa ra được 12 bản đồ phân bố tần suất xuất hiện front theo các tháng trong năm. Vào mùa Đông, các front chủ yếu liên quan đến hiện tượng lưỡi nước lạnh xâm nhập sâu xuống phía Nam, hình thành các khu vực front chủ yếu: Tây Philippin, eo Đài Loan, Phía Đông đảo Hải Nam, Vịnh Bắc Bộ, Nam Trung Bộ và - Cửa vịnh Thái Lan. Vào mùa hè: Front xuất hiện chính ở: Eo Đài loan, quanh đảo Hải Nam và Nam Trung Bộ. Front do nước trồi xuất hiện ở khu vực Nam Trung Bộ có thời gian tồn tại từ 5 ngày cho tới hơn 10 ngày. Đây là một điểm cần chú ý khi nghiên cứu chi tiết hơn về loại front này ở khu vực Nam Trung Bộ. 4. Trong tương lai, để phát triển thêm những nghiên cứu về front nhiệt mặt biển cần nghiên cứu sâu hơn với độ phân giải thời gian ngắn hơn (5 ngày) cũng như cần có sự kết hợp với các số liệu thực đo để xác định rõ hơn phân bố theo chiều ngang cũng như theo chiều thẳng đứng của front. Đồng thời cũng cần hoàn thiện hơn nữa quy trình nghiên cứu và tính toán front nhiệt mặt biển. 73 TÀI LIỆU THAM KHẢO Tài liệu tham khảo tiếng Anh 1. Belkin I.M, 2003. (Graduate School of Oceanography, University of Rhode Island). Front Encyclopedia, Interdisciplinary Encyclopedia of Marine Sciences. 2. Belkin I.M, Cornillon P. (2003). SST fronts of the Pacific coastal and marginal seas, Physical Oceanography, Vol. 1, No. 2, 3. Belkin I.M , John E. O'Reilly (2009). An algorithm for oceanic front detection in chlorophyll and SST satellite imagery, Journal of Marine Systems, doi:10.1016/j.jmarsys.2008.11.018 4. Belkin I.M., Cornillon P. and Ullman D. (2003). Ocean fronts around Alaska from satellite SST data, Proceedings of the Amer. Met. Soc. 7 th Conf. on the Polar Meteorology and Oceanography, Hyannis, MA, Paper 12.7, 15pp 5. Belkin I.M, Cornillon P.C. (2004). Surface thermal fronts of the Okhotsk Sea, PACIFIC OCEANOGRAPHY, Vol. 2, No. 1-2. 6. Belkin I.M, Cornillon Peter C., and Kenneth Sherman (2008). Fronts in Large Marine Ecosystems of the World Ocean, Re-submitted to the Progress in Oceanography’ Special Issue on Large Marine Ecosystems. 7. Cayula, J.-F. and Cornillon P. (1992). Edge detection algorithm for SST images. Journal of Atmospheric and Oceanic Technology 9: 67-80. 8. Cayula, J.-F. and Cornillon P. (1995). Multi-Image edge Detection for SST images, Journal of atmospheric and oceanic technology, Vol 12. 9. Cayula, J.-F. and Cornillon P. (1996). Cloud detection from a sequence of SST images, Remote Sensing of Environment, 55(1), pp. 80–88. 10. Dongxiao Wang, Yun Liu, Yiquan Qi, Ping Shi (2001). Seasonal Variability of Thermal Fronts in the Northern South China Sea from Satellite Data, Geophysical research letters, VOL. 28, NO. 20, PAGES 3963-3966. 11. MGET website: 74 12. Kuo-Wei Lan, Hiroshi Kawamura, Ming-An Lee, Yi Chang, Jui-Wen Chan, Cheng-Hsin Liao (2009). Summertime sea surface temperature fronts associated with upwelling around the Taiwan Bank, Continental Shelf Research, doi:10.1016/j.csr.2009.01.015 13. Teruhisa Shimada , Futoki Sakaida, Hiroshi Kawamura, Toshiaki Okumura (2005). Application of an edge detection method to satellite images for distinguishing sea surface temperature fronts near the Japanese coast, Remote Sensing of Environment 98, 21–34 14. Y. Chang; M. -A. Lee; T. Shimada; F. Sakaida; H. Kawamura; J. -W. Chan; H. -J. Lu (2008). Wintertime high-resolution features of sea surface temperature and chlorophylla fields associated with oceanic fronts in the southern East China Sea, International Journal of Remote Sensing. Tài liệu tham khảo tiếng Việt 15. Cơ sở khoa học cho việc dự báo phân bố và biến động nguồn lợi cá biển, 1995, Báo cáo Tổng kết đề tài KT 03-10, do GS Lê Đức Tố chủ trì 16. Đinh Văn Ưu & cs (2000), Báo cáo chuyên đề 7 đề tài KHCN 06-02: Các đặc trưng chế độ hải dương học nghề cá Biển Đông.