Tổng quan luận án

Hiện tượng spread F xích đạo (Equatorial Spread F – ESF) là trạng thái bất ổn định plasma theo phương ngang xảy ra tại tầng điện ly lớp F vào ban đêm ở khu vực vĩ độ thấp. Hiện tượng này làm biến dạng cấu trúc phân tầng ngang, gây suy giảm phẩm chất hoặc làm gián đoạn các tín hiệu vô tuyến phản xạ từ tầng điện ly cũng như các tín hiệu vệ tinh truyền xuyên qua môi trường plasma. Do chưa có các giải pháp kỹ thuật triệt tiêu hoàn toàn sự can thiệp này, việc nâng cao năng lực dự báo thời điểm và không gian xuất hiện của ESF trở thành yêu cầu thiết yếu đối với các hệ thống thông tin liên lạc và định vị - dẫn đường vệ tinh.

Sau hơn tám thập kỷ nghiên cứu kể từ các quan trắc đầu tiên, nguyên nhân gây ra sự biến thiên hàng ngày của hiện tượng spread F vẫn chưa có một hệ thống lý thuyết hoàn chỉnh để giải thích thấu đáo. Các nghiên cứu hiện đại tập trung vào hai yếu tố điều khiển chính: sự nâng lên sau hoàng hôn của lớp F (Post-Sunset Rise of F layer – PSSR) và cấu trúc dạng sóng quy mô lớn (Large-Scale Wave Structure – LSWS). Trong khi PSSR đã được khảo sát sâu rộng nhờ chuỗi số liệu thăm dò mặt đất và vệ tinh phong phú, LSWS mới chỉ được chú ý trong khoảng mười năm gần đây với nguồn dữ liệu hạn chế từ radar ALTAIR và vệ tinh chuyên dụng C/NOFS. Bên cạnh đó, tại dải kinh độ bao gồm Việt Nam, Indonesia và Philippines, các dữ liệu quan trắc vệ tinh toàn cầu cho thấy sự xuất hiện của ESF mang những đặc tính dị biệt mà các mô hình áp dụng cho những khu vực khác chưa giải thích được. Trong bối cảnh số liệu quan trắc chuyên dụng tại Việt Nam còn thiếu thốn, việc kết hợp số liệu từ các trạm quan trắc quốc tế có mật độ đo đạc cao để hoàn thiện phương pháp luận, sau đó áp dụng phân tích chuỗi dữ liệu trong nước là hướng tiếp cận được lựa chọn.

Mục tiêu nghiên cứu của luận án được xác định gồm hai nội dung chính:

  1. Làm sáng tỏ phương pháp trích lọc thông tin từ thăm dò thẳng đứng (điện ly đồ) nhằm mô tả gián tiếp và quan trắc hai hình thái vào giai đoạn đầu tiên và cuối cùng của spread F trong hoàn cảnh không có số liệu mô tả trực tiếp, là cấu trúc dạng sóng quy mô lớn và bong bóng plasma xích đạo (Equatorial Plasma Bubble – EPB).
  2. Đóng góp thêm các bằng chứng nhằm tiếp tục làm sáng tỏ các quan điểm hiện tại về nguyên nhân gây biến thiên và đánh giá vai trò của các yếu tố điều khiển ESF tại khu vực nghiên cứu.

Đối tượng nghiên cứu của luận án bao gồm:

  1. Hiện tượng spread F và hai nguồn gây ra sự biến đổi hoạt động này theo thời gian trong phạm vi các đảo thuộc Trung tâm Thái Bình Dương.
  2. Hiện tượng spread F và cấu trúc bong bóng plasma xích đạo tại Việt Nam cùng khu vực lân cận.

Phạm vi nghiên cứu:

  • Không gian: Trọng tâm nghiên cứu thực nghiệm đặt tại khu vực đảo Kwajalein thuộc Trung tâm Thái Bình Dương (vùng xích đạo từ), có đối sánh với các trạm Kodaikanal, Trivandrum (Ấn Độ), Jicamarca (Peru), Agua Verde; mạng lưới quan trắc tại Việt Nam gồm trạm Phú Thụy, Hóc Môn, Bạc Liêu, liên kết đối sánh với Chumphon (Thái Lan) và Baguio (Philippines).
  • Thời gian: Khảo sát tập trung vào chu kỳ hoạt động Mặt Trời cao (năm 2000) và thấp (năm 2011, chi tiết từ tháng 01/2011 đến tháng 01/2012); các chuỗi số liệu lịch sử tại Phú Thụy (1962–1979 và 2007–2010), Hóc Môn (2002–2006), Baguio (1953–1955); dữ liệu đo đạc vệ tinh C/NOFS (giai đoạn 2009–2011) và vệ tinh DMSP (1999–2002).

Tổng quan tài liệu và vị trí của luận án

Hiện tượng ESF được ghi nhận lần đầu vào năm 1934 và chính thức công bố bởi Booker và Wells (1938) thông qua các quan sát điện ly đồ ghi nhận vết phản xạ khuếch tán theo độ cao và theo tần số tại đài quan trắc Huancayo (Peru). Về mặt cơ chế vật lý, Dungey (1956) đặt nền móng lý thuyết đầu tiên khi áp dụng cơ chế bất ổn định Rayleigh - Taylor hấp dẫn (Gravitational Rayleigh - Taylor – GRT) để giải thích sự hình thành các mảng bất đồng nhất plasma tại đáy lớp F. Do lý thuyết GRT tuyến tính cổ điển không thể giải thích sự thâm nhập của cấu trúc bất ổn định lên vùng đỉnh lớp F, các công trình tiếp theo của Scannapieco và Ossakow (1976), Woodman và La Hoz (1976) đã phát triển lý thuyết GRT phi tuyến, mô tả sự hình thành các cấu trúc dạng lông vũ (plume) và bong bóng plasma (EPB) di chuyển lên độ cao lớn. Đến các nghiên cứu của Kelley (1989, 2009) và Sultan (1996), lý thuyết GRT tổng quát được hoàn thiện bằng cách tích hợp thêm các yếu tố điện trường nền, chuyển động kéo theo của gió trung hòa và độ nghiêng của đáy lớp F.

Về yếu tố PSSR, Farley et al. (1986), Rishbeth (1971), Heelis et al. (1974) và Eccles (1998) đã nghiên cứu các cơ chế động lực học hình thành sự tăng cường của điện trường hướng Đông trước thời điểm đảo chiều (Pre-Reversal Enhancement – PRE). Fejer et al. (1979, 1989, 1991, 1999) đã hệ thống hóa các đặc trưng thực nghiệm và xây dựng mô hình kinh nghiệm tính toán vận tốc trượt thẳng đứng của plasma (mô hình Scherliess - Fejer 1999). Các tác giả Abdu et al. (1983), Anderson et al. (2004) đi sâu vào việc xác định các giá trị ngưỡng về độ cao và vận tốc nâng lên cực đại của lớp F cần thiết để kích hoạt bất ổn định. Tsunoda (1985) đề xuất lý thuyết về sự thẳng hàng giữa ranh giới ngày - đêm và đường sức từ trường Trái Đất (Solar Terminator - B Alignment – STBA), chứng minh rằng khi xảy ra hoàng hôn đồng thời ở hai vùng chân đường sức lớp E liên hợp, gradient độ dẫn Pedersen dọc đường sức đạt cực đại, tạo điều kiện thuận lợi nhất cho PSSR và ESF phát triển vào các kỳ phân điểm. Tuy nhiên, các quan trắc vệ tinh DMSP của Burke et al. (2004) và Gentile et al. (2006) chỉ ra sự tồn tại của cực đại ESF vào các kỳ chí điểm tại khu vực Trung tâm Thái Bình Dương, Châu Phi và Nam Mỹ, tạo nên mâu thuẫn mà mô hình STBA đơn thuần không giải thích được.

Về yếu tố nguồn kích thích ban đầu, Rottger (1973, 1976) phát hiện dấu hiệu của các mảng bất đồng nhất có tính chu kỳ thông qua phương pháp truyền sóng xuyên xích đạo (TEP). Hysell et al. (1990, 1994) và Tsunoda et al. (2010, 2011) mô tả trực tiếp các biến điệu dạng sóng ở đáy lớp F (LSWS) bằng radar ALTAIR và ảnh chụp quang học phát quang khí quyển (airglow). Cơ chế kích hoạt LSWS được Tsunoda (2010) giải thích qua lý thuyết thẳng hàng giữa đường sức từ trường và mặt đầu sóng trọng lực khí quyển (Gravity Wave - B Alignment – GWBA), theo đó sóng trọng lực (Atmospheric Gravity Wave – AGW) xuất phát từ các hoàn lưu đối lưu quy mô trung bình (MCC) trong dải hội tụ nội nhiệt đới (ITCZ) truyền lên tầng nhiệt quyển và kích thích plasma đáy lớp F. Thêm vào đó, Tsunoda et al. (2010, 2013) bước đầu liên hệ các vết phản xạ dị thường trên điện ly đồ như vết phản xạ nhiều lần (Multi-Reflected Echo – MRE) và vết phụ (Satellite Trace – ST) với sự xuất hiện của cấu trúc vòm nâng (upwelling).

Khoảng trống nghiên cứu mà luận án lựa chọn giải quyết là:

  • Đa số các công trình trước đây chỉ đánh giá đơn lẻ vai trò của PSSR hoặc LSWS; luận án tiến hành nghiên cứu đồng thời hai yếu tố này ở quy mô thời gian hàng tháng để đánh giá mức độ đóng góp của từng thành phần đối với sự phát triển của ESF.
  • Khắc phục sự phụ thuộc vào các thiết bị chuyên dụng đắt tiền (radar ALTAIR, vệ tinh C/NOFS) bằng cách chuẩn hóa phương pháp trích xuất thông tin về LSWS thông qua dấu hiệu vết phản xạ nhiều lần (MRE/MRD) trên điện ly đồ của máy thăm dò thẳng đứng thông thường.
  • Làm sáng tỏ bản chất của các cấu trúc spread F phức hợp ghi nhận tại trạm nằm xa xích đạo từ (Phú Thụy, Việt Nam), xác định mối liên hệ vật lý giữa các cấu trúc này với sự lan truyền của bong bóng plasma từ vùng xích đạo từ.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án dựa trên khung lý thuyết điện động lực học tầng điện ly vùng vĩ độ thấp. Chuyển động của môi trường plasma được mô tả bằng hệ phương trình Maxwell kết hợp phương trình chuyển động của hạt mang điện dưới tác dụng của lực Lorentz và trọng lực:

$$\vec{j} = \sigma_0 \vec{E}\parallel + \sigma_P \vec{E}\perp + \sigma_H \frac{\vec{B} \times \vec{E}}{B}$$

Trong đó, $\sigma_0$, $\sigma_P$, và $\sigma_H$ lần lượt là độ dẫn điện song song, độ dẫn Pedersen và độ dẫn Hall. Tại lớp E xích đạo vào ban ngày, sự ngăn chặn của dòng Hall theo phương thẳng đứng do ranh giới mật độ sinh ra điện trường phân cực, làm tăng độ dẫn tổng cộng lên thành độ dẫn Cowling $\sigma_C$, hình thành dòng điện xích đạo (Equatorial Electrojet – EEJ). Vào ban đêm, quá trình dynamo lớp F chiếm ưu thế khi lớp E suy giảm ion hóa ("hở mạch"), duy trì điện trường phân cực hướng xuống dưới và điều khiển dòng plasma trôi sang phía Đông.

Tốc độ tăng trưởng tuyến tính của bất ổn định Rayleigh - Taylor tổng quát ($\gamma_{GRT}$) được tác giả sử dụng dưới dạng biểu thức:

$$\gamma_{GRT} = \frac{1}{L} \left( \frac{g}{\nu_{in}} \cos\Delta + \frac{E_x - U_z B}{B} \cos\Delta + \frac{E_z + U_x B}{B} \sin\Delta \right)$$

Biểu thức này cho thấy tốc độ phát triển bất ổn định tỉ lệ nghịch với độ dài tỉ lệ gradient mật độ electron ($L = N / (\partial N/\partial z)$), tỉ lệ thuận với gia tốc trọng trường ($g$), vận tốc trượt thẳng đứng hướng lên của plasma trước khi đảo chiều ($V_z = E_x / B$), vận tốc gió trung hòa vĩ hướng ($U_x$) và góc nghiêng của đáy lớp F ($\Delta$).

Phương pháp nghiên cứu thực nghiệm trong luận án kết hợp kỹ thuật phân tích số liệu đa nguồn:

  • Phân tích điện ly đồ (Ionosonde): Xử lý tín hiệu quét tần số cao (1–30 MHz) từ máy thăm dò thẳng đứng. Các thông số được trích xuất gồm độ cao biểu kiến đáy lớp F ($h'F$), tần số tới hạn ($f_oF_2$), vận tốc nâng lên biểu kiến ($V_{PSSR} = \Delta h'F / \Delta t$). Nhận diện hiện tượng spread F trải rộng theo độ cao (Range Spread F – RSF) khi bề dày vết khuếch tán vượt quá 30 km và trải rộng tần số khi độ giãn vượt quá 0,3 MHz. Nhận diện vết phản xạ nhiều lần (MRE, MRD) và vết phụ (ST) dựa trên độ dốc và độ trễ phản xạ so với vết 1F.
  • Kỹ thuật tính nồng độ điện tử tổng cộng (TEC) từ sóng vô tuyến dẫn đường vệ tinh: Thu nhận tín hiệu kết hợp hai tần số kết hợp 150 MHz và 400 MHz phát từ thiết bị CERTO trên vệ tinh C/NOFS thông qua máy thu vô tuyến số GRBR. Sử dụng phương pháp một trạm của Yamamoto et al. và thuật toán bình phương tối thiểu phi tuyến để xác định độ lệch pha ban đầu và tính toán giá trị TEC thẳng đứng tuyệt đối. Thành phần nhiễu loạn PTEC được tách lọc bằng cách lấy hiệu số giữa giá trị TEC đo được với đường biến thiên nền tính bằng trung bình trượt trong 2,5 phút (tương ứng khoảng cách ~1000 km theo quỹ đạo bay).
  • Phân tích bức xạ sóng dài (OLR): Khai thác cơ sở dữ liệu bức xạ hồng ngoại nhiệt quét hai kênh từ vệ tinh khí tượng NOAA trên lưới tọa độ $2,5^\circ \times 2,5^\circ$. Luận án sử dụng ngưỡng thông lượng $230\text{ W/m}^2$ để phân lập các tâm đối lưu mạnh trong tầng đối lưu nhiệt đới, làm căn cứ định vị nguồn kích thích sóng trọng lực khí quyển.
  • Mô hình hóa kinh nghiệm: Sử dụng mô hình Scherliess - Fejer (SF99) để tính toán vận tốc trôi thẳng đứng trung bình của plasma lớp F và mô hình HWM93 để đối sánh thành phần vận tốc gió trung hòa vĩ hướng.

Bảng tổng hợp các nguồn tư liệu và dữ liệu quan trắc sử dụng trong luận án:

Nguồn dữ liệu / Thiết bị Vị trí địa lý / Nền tảng Tham số vật lý trích lọc Thời gian khảo sát
Máy thăm dò điện ly thẳng đứng Đài Kwajalein (Trung tâm Thái Bình Dương) Điện ly đồ: $h'F$, $f_oF_2$, RSF, MRE, MRD, ST, $V_{PSSR}$ 01/2011 – 01/2012
Máy thăm dò điện ly thẳng đứng Phú Thụy, Hóc Môn, Bạc Liêu (Việt Nam) Điện ly đồ: $h'F$, cấu trúc spread F phức hợp, vết phụ ST 1962–1979; 2002–2006; 2007–2010
Máy thăm dò điện ly thẳng đứng Trivandrum (Ấn Độ), Baguio (Philippines) Vận tốc trượt thẳng đứng, tần suất xuất hiện ESF ($POF_{SF}$) 1953–1955; 2004
Vệ tinh C/NOFS (đầu thu GRBR/CERTO) Quỹ đạo thấp (độ nghiêng $13^\circ$, độ cao 375–710 km) Pha tín hiệu 150/400 MHz, nồng độ điện tử tổng cộng (TEC, PTEC) 2009–2011
Vệ tinh DMSP Vệ tinh khí tượng quốc phòng Tần suất xuất hiện bong bóng plasma xích đạo (EPB) 1999–2002; 2011
Vệ tinh ROCSAT-1, AE-E Vệ tinh khoa học quỹ đạo thấp Tần suất xuất hiện bất ổn định mật độ plasma Các năm hoạt động Mặt Trời cao
Dữ liệu vệ tinh NOAA Toàn cầu (lưới $2,5^\circ \times 2,5^\circ$) Thông lượng bức xạ sóng dài (OLR), ngưỡng đối lưu $230\text{ W/m}^2$ 2011–2012
Mô hình SF99 và HWM93 Mô hình tính toán số trị Vận tốc trượt thẳng đứng lý thuyết ($V_{PSSR}$), gió trung hòa ($U$) Chu kỳ cực đại (2000) và cực tiểu (2011)

Nội dung chính theo từng chương

Chương 1: Tổng quan về điện ly vùng xích đạo và spread F

Chương 1 hệ thống hóa các đặc trưng cấu trúc vật lý và điện động lực học của tầng điện ly xích đạo. Tác giả phân tích sự hình thành các lớp ion hóa D, E, F theo độ cao dựa trên cân bằng quang hóa và động học va chạm giữa các hạt mang điện và hạt trung hòa. Tại khu vực xích đạo từ ($\pm 10^\circ$ vĩ độ từ), từ trường Trái Đất có phương nằm ngang tạo ra tính dị hướng cao trong độ dẫn điện. Ban ngày, dynamo lớp E tạo nên dòng điện xích đạo EEJ hướng Đông với độ dẫn Cowling tăng cường. Ban đêm, dynamo lớp F duy trì sự trôi dạt của plasma sang hướng Đông và hình thành điện trường phân cực hướng xuống dưới.

Nội dung trọng tâm của chương là phân tích cơ chế hình thành spread F xích đạo thông qua bất ổn định Rayleigh - Taylor hấp dẫn (GRT) phi tuyến và tổng quát. Luận án phân loại các hình thái bất ổn định quan trắc được thành: cấu trúc vòm nâng (upwelling) ở đáy lớp F có quy mô ngang 100–1000 km; cấu trúc dạng lông vũ (plume) phát triển ở sườn vòm nâng; và cấu trúc bong bóng plasma xích đạo (EPB) có kích thước 10–100 km thâm nhập xuyên qua cực đại lớp F lên độ cao hàng nghìn kilômét. Phân tích lý thuyết chỉ ra hai yếu tố điều khiển xuất hiện hàng ngày của hiện tượng là sự nâng lên sau hoàng hôn của lớp F (PSSR) và cấu trúc dạng sóng quy mô lớn (LSWS). Luận án cũng trình bày chi tiết các hạn chế của lý thuyết STBA đơn thuần khi không giải thích được cực đại xuất hiện ESF vào các kỳ chí điểm tại một số dải kinh độ, đồng thời giới thiệu mô hình tích hợp của Tsunoda về sự tương tác giữa PSSR và LSWS dưới tác động kích hoạt của sóng trọng lực khí quyển (AGW) sinh ra từ hoàn lưu đối lưu tầng đối lưu (ITCZ). Cuối chương, tác giả mô tả nguyên lý vận hành của các phương pháp thăm dò được sử dụng xuyên suốt luận án gồm ionosonde, hệ đo pha vệ tinh C/NOFS và dữ liệu OLR.

Chương 2: Cấu trúc dạng sóng quy mô lớn từ điện ly đồ và vai trò của nó đối với spread F

Chương 2 tập trung chứng minh phương pháp phát hiện cấu trúc dạng sóng quy mô lớn (LSWS) từ máy thăm dò điện ly thẳng đứng, giải quyết bài toán thiếu hụt thiết bị thăm dò không gian đắt tiền. Tác giả phân tích chi tiết diễn biến tầng điện ly trong đêm điển hình ngày 24/04/2011 tại trạm Kwajalein. Trên điện ly đồ, dấu hiệu ban đầu xuất hiện lúc 18h47 giờ địa phương (LT) – gần thời điểm hoàng hôn lớp E – dưới dạng các vết phản xạ nhiều lần (Multi-Reflected Echo – MRE) và vết kép phản xạ nhiều lần (Multi-Reflected Doublet – MRD), có độ dốc lớn hơn đáng kể so với vết 1F thông thường. Tiếp đó, vết phụ (ST) được ghi nhận lúc 20h32 LT trước khi hiện tượng spread F phát triển mạnh mẽ.

Để xác nhận bản chất vật lý của các vết phản xạ này, tác giả đối chiếu dữ liệu điện ly đồ với số liệu nồng độ điện tử tổng cộng (TEC) từ ba đường bay liên tiếp của vệ tinh C/NOFS qua khu vực đảo Kwajalein vào các thời điểm 17h47, 19h31 và 21h15 LT cùng đêm. Kết quả phân tích thành phần nhiễu loạn PTEC ghi nhận các dao động dạng sóng rõ rệt với bước sóng không gian khoảng 1000 km, biên độ biến thiên đạt tới vài đơn vị TECu. Quá trình mô hình hóa đường truyền tia vô tuyến khẳng định: khi đáy lớp F bị biến điệu dạng sóng tạo thành các cấu trúc vòm nâng và lòng máng lõm, các tia thăm dò bị phản xạ xiên và hội tụ pha. Nếu chu kỳ phát xung (IPP) của máy thăm dò ngắn hơn thời gian trễ của tín hiệu phản xạ bậc cao qua lại giữa tầng điện ly và mặt đất, tín hiệu phản xạ bậc $n$ sẽ bị biểu diễn sai lệch về độ cao, tạo thành vết MRE/MRD. Khi cấu trúc vòm nâng tiếp tục trôi dạt về phía Đông, tia phản xạ xiên từ đáy vòm nâng ghi nhận tách biệt với vết thẳng đứng 1F, tạo ra vết phụ ST. Kết quả này khẳng định MRE/MRD và ST trên điện ly đồ chính là dấu hiệu chỉ thị trực tiếp của LSWS ở đáy lớp F trước khi bất ổn định bùng phát.

Chương 3: Vai trò của sự nâng lên sau hoàng hôn của lớp F và cấu trúc dạng sóng quy mô lớn với spread F

Chương 3 trình bày các kết quả phân tích thống kê đồng thời vai trò của PSSR và LSWS đối với hoạt động của spread F ở quy mô hàng tháng tại khu vực đảo Kwajalein trong chu kỳ một năm (01/2011 – 01/2012), kết hợp so sánh với các chuỗi số liệu ở thời kỳ hoạt động Mặt Trời cao (năm 2000). Về tham số PSSR, tác giả trích xuất vận tốc nâng lên thẳng đứng cực đại ($V_{PSSR}$) từ điện ly đồ, đối chiếu với mô hình Scherliess - Fejer (SF99) và số liệu quan trắc tại trạm Trivandrum. Kết quả thực nghiệm chỉ ra rằng hoạt động của ESF có mối tương quan chặt chẽ với độ lớn của $V_{PSSR}$. Khi $V_{PSSR}$ vượt qua giá trị ngưỡng xác định, tốc độ phát triển bất ổn định $\gamma_{GRT}$ tăng nhanh, tạo điều kiện thuận lợi cho mảng plasma bùng phát.

Khảo sát tần suất xuất hiện của vết phản xạ nhiều lần (MRE) trong các tháng 4, 7 và 8 năm 2011 cho thấy MRE xuất hiện với tần suất cao trong khoảng thời gian từ hoàng hôn lớp E đến trước khi hình thành spread F. Khi kết hợp phân tích cả hai yếu tố, tác giả phát hiện: trong những đêm ghi nhận đồng thời giá trị $V_{PSSR}$ cao và sự hiện diện của MRE (chỉ thị LSWS), spread F luôn bùng phát mạnh; ngược lại, nếu chỉ có một trong hai yếu tố xuất hiện đơn lẻ, spread F phát triển yếu hoặc không hình thành. Đáng chú ý, kết quả nghiên cứu chỉ ra rằng vận tốc nâng lên cực đại $V_{PSSR}$ đóng vai trò kiểm soát quan trọng hơn so với điều kiện hình học của nút hoàng hôn (STBA). Về nguyên nhân điều khiển biến thiên của $V_{PSSR}$, luận án đưa ra bằng chứng cho thấy có sự phân hóa theo mùa: vào các kỳ phân điểm, biến thiên cực đại của $V_{PSSR}$ bị chi phối chủ yếu bởi cường độ của dòng điện xích đạo EEJ ban ngày; trong khi đó, vào các kỳ chí điểm (đặc biệt là hạ chí), biến thiên này chịu sự điều hòa mạnh mẽ của lực triều khí quyển số sóng 2 kết hợp với hoạt động đối lưu nhiệt đới (được định vị qua bản đồ bức xạ sóng dài OLR).

Chương 4: Cấu trúc bong bóng plasma xích đạo quan trắc tại Việt Nam

Chương 4 giải quyết bài toán thực tế về đặc điểm dị thường của spread F ghi nhận tại trạm Phú Thụy (Việt Nam), một trạm quan trắc nằm ở vĩ độ thấp nhưng cách xa xích đạo từ khoảng $10^\circ$. Tác giả thực hiện so sánh tổng hợp giữa tần suất xuất hiện spread F ($POF_{SF}$) tại Phú Thụy giai đoạn 2007–2010 với số liệu lịch sử tại trạm Baguio (Philippines, 1953–1955), trạm xích đạo Chumphon (Thái Lan), trạm Bạc Liêu, cùng các quan sát bong bóng plasma từ vệ tinh C/NOFS và DMSP trong cùng dải kinh độ.

Phân tích chuỗi điện ly đồ tại Phú Thụy cho thấy các vết spread F xuất hiện dưới dạng phân tầng phức hợp, có sự hiện diện của các vết phụ kéo dài và độ khuếch tán không đồng nhất, hoàn toàn khác biệt với dạng vết trải rộng độ cao (RSF) đồng nhất quan sát được tại các trạm nằm ngay trên xích đạo từ. Phân tích chi tiết hai chuỗi quan trắc ban đêm ngày 24/06/2007 và 25/07/2010 chỉ ra rằng các cấu trúc spread F phức hợp này không phải là bất ổn định phát sinh tại chỗ ở đáy lớp F phía trên trạm Phú Thụy. Thay vào đó, chúng là biểu hiện quang hình của các phản xạ sóng vô tuyến theo phương nghiêng hướng về các cấu trúc bong bóng plasma (EPB) hình thành tại vùng đỉnh xích đạo từ và phát triển dọc theo các đường sức từ trường Trái Đất xuống các vĩ độ thấp hơn. Sự dịch chuyển ngang của các vách bong bóng plasma qua búp sóng của máy thăm dò thẳng đứng đã tạo ra các vết phản xạ xiên với độ trễ biến thiên, giải thích thỏa đáng nguồn gốc các dạng spread F dị thường tại miền Bắc Việt Nam mà các nghiên cứu trước đây chưa làm rõ.

Kết quả và những đóng góp mới

Các đóng góp mới về mặt lý luận và khoa học của luận án bao gồm:

  • Khẳng định vết phản xạ nhiều lần (gồm dạng đơn MRE và dạng kép MRD) ghi nhận trên điện ly đồ là dấu hiệu nhận biết gián tiếp của cấu trúc dạng sóng quy mô lớn (LSWS) ở đáy lớp F. Bằng chứng này xác nhận cơ chế hội tụ phản xạ vô tuyến từ sườn cấu trúc vòm nâng dạng sóng trước thời điểm bùng phát bất ổn định.
  • Chứng minh rằng trong mối quan hệ giữa hai yếu tố điều khiển đối với biến thiên spread F ở quy mô hàng tháng, vận tốc nâng lên thẳng đứng cực đại của lớp F ($V_{PSSR}$) giữ vai trò chi phối quan trọng hơn gradient độ dẫn Pedersen dọc kinh tuyến (nút hoàng hôn trong mô hình STBA).
  • Đưa ra giả thuyết và luận cứ chỉ rõ sự tồn tại của hai cơ chế vật lý khác nhau điều khiển biến thiên cực đại của $V_{PSSR}$: sự biến đổi cường độ của dòng điện xích đạo EEJ kiểm soát chủ đạo trong các kỳ phân điểm, trong khi lực triều khí quyển số sóng 2 kết hợp năng lượng đối lưu nhiệt đới chi phối trong các kỳ chí điểm.
  • Xác định bản chất vật lý của các dạng spread F phức hợp quan trắc tại trạm Phú Thụy (vùng xa xích đạo từ) là dấu hiệu phản xạ từ các cấu trúc bong bóng plasma xích đạo (EPB) phát triển dọc theo đường sức từ từ xích đạo truyền tới, phân biệt rõ ràng với cơ chế bất ổn định tại chỗ ở vùng xích đạo từ.

Kiến nghị và giải pháp đề xuất từ luận án:

  • Kiến nghị thiết lập quy trình khai thác các dấu hiệu MRE, MRD và ST trên mạng lưới máy thăm dò điện ly thẳng đứng mặt đất như một chỉ số cảnh báo sớm (trước từ 30 phút đến 1 giờ) cho sự hình thành của bong bóng plasma và hiện tượng nhấp nháy sóng vô tuyến.
  • Đề xuất tăng cường kết hợp số liệu giữa máy thăm dò thẳng đứng hiện có tại Bạc Liêu với hệ thống máy thu tín hiệu GNSS/beacon đa trạm trên toàn quốc để xây dựng mô hình giám sát biến động tầng điện ly phục vụ công tác định vị và truyền thông không gian tại Việt Nam.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án:

  • Khảo sát thực nghiệm định lượng đồng thời hai yếu tố điều khiển (PSSR và LSWS) mới chỉ thực hiện tập trung trên chuỗi số liệu một năm (từ tháng 01/2011 đến tháng 01/2012) tại khu vực đảo Kwajalein do giới hạn về thời gian hoạt động của các thiết bị đo đạc chuyên dụng và dữ liệu vệ tinh C/NOFS.
  • Mạng lưới trạm thăm dò thẳng đứng hoạt động liên tục tại Việt Nam còn mỏng, số liệu sử dụng tại trạm Phú Thụy và Hóc Môn là các chuỗi đo đạc không liên tục trong quá khứ, gây khó khăn nhất định cho việc phân tích thống kê dài hạn theo nhiều chu kỳ hoạt động Mặt Trời khác nhau.

Hướng nghiên cứu tiếp theo:

  • Mở rộng phạm vi kiểm chứng mối tương quan giữa $V_{PSSR}$, LSWS và ESF trên các chuỗi số liệu đo đạc nhiều năm, bao quát trọn vẹn cả thời kỳ hoạt động Mặt Trời cực đại và cực tiểu của các chu kỳ tiếp theo.
  • Ứng dụng các mô hình mô phỏng số ba chiều hoàn chỉnh về động lực học tầng điện ly để định lượng hóa tác động truyền năng lượng của sóng trọng lực từ tầng đối lưu qua ranh giới trung hòa - ion hóa lên đáy lớp F.
  • Đầu tư khôi phục và mở rộng mạng lưới máy thăm dò điện ly số thế hệ mới kết hợp máy thu định vị đa vệ tinh tại các trạm trọng điểm ở Việt Nam nhằm nâng cao độ chính xác của bài toán dự báo thời tiết không gian.

Giá trị tham khảo

Luận án cung cấp hệ thống dữ liệu thực nghiệm và phương pháp luận hữu ích cho các đối tượng và lĩnh vực sau:

  • Các nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Vật lý địa cầu, Vật lý khí quyển, Khoa học vũ trụ và Vô tuyến điện tử: Luận án là tài liệu tham khảo chi tiết về điện động lực học tầng điện ly xích đạo, lý thuyết bất ổn định plasma (GRT), cơ chế dynamo và kỹ thuật phân tích điện ly đồ chuyên sâu.
  • Các kỹ sư và chuyên gia vận hành hệ thống thông tin liên lạc sóng ngắn (HF), hệ thống dẫn đường - định vị vệ tinh (GPS, GLONASS, Galileo) và truyền thông vệ tinh: Kết quả nghiên cứu cung cấp cơ sở để đánh giá mức độ suy giảm chất lượng tín hiệu do hiện tượng nhấp nháy ion hóa, hỗ trợ tối ưu hóa thuật toán hiệu chỉnh trễ tầng điện ly.
  • Các cơ quan quản lý và dự báo thời tiết không gian: Phương pháp trích xuất các dấu hiệu sớm (MRE, MRD, ST) trên điện ly đồ mở ra phương án chi phí thấp để thiết lập hệ thống cảnh báo sớm bất ổn định điện ly phục vụ an ninh quốc phòng và kinh tế biển.

Câu hỏi thường gặp

  1. Vết phản xạ nhiều lần (MRE) và vết kép phản xạ nhiều lần (MRD) trên điện ly đồ phản ánh hiện tượng vật lý gì? Vết MRE và MRD là các vết phản xạ có độ dốc lớn hơn bình thường, xuất hiện ở dải tần số thấp xung quanh thời điểm hoàng hôn lớp E. Về mặt vật lý, chúng là biểu hiện của các sóng phản xạ vô tuyến bị lệch góc và hội tụ pha khi gặp cấu trúc vòm nâng dạng lòng máng lõm của cấu trúc sóng quy mô lớn (LSWS) ở đáy lớp F. Khi chu kỳ phát xung (IPP) của máy thăm dò ngắn hơn thời gian trễ của tín hiệu phản xạ bậc cao giữa điện ly và mặt đất, tín hiệu phản xạ bị máy thu nhận diện sai về độ cao biểu kiến, hiển thị thành các vết nF dịch chuyển hoặc tách đôi thành vết kép MRD.

  2. Luận án xác định vai trò tương đối giữa PSSR và LSWS trong việc kích hoạt spread F như thế nào? Luận án chỉ ra rằng cả hai yếu tố đều tham gia vào quá trình hình thành ESF nhưng giữ vai trò khác nhau: LSWS cung cấp mầm mống nhiễu loạn ban đầu làm biến điệu đáy lớp F thành các vòm nâng, trong khi PSSR cung cấp điều kiện nền thông qua việc nâng lớp F lên độ cao có tần số va chạm ion - trung hòa thấp để vận tốc nâng lên cực đại ($V_{PSSR}$) khuếch đại tốc độ tăng trưởng bất ổn định Rayleigh - Taylor. Hoạt động của ESF đạt mức cao nhất khi có sự xuất hiện đồng thời của cả mầm mống LSWS và vận tốc $V_{PSSR}$ vượt ngưỡng; trong đó $V_{PSSR}$ giữ vai trò kiểm soát trực tiếp và quan trọng hơn so với điều kiện hình học nút hoàng hôn.

  3. Tại sao lý thuyết STBA đơn thuần không giải thích được sự xuất hiện của ESF vào các kỳ chí điểm tại khu vực Trung tâm Thái Bình Dương? Lý thuyết STBA giả định rằng ESF chỉ xuất hiện cực đại vào các kỳ phân điểm khi ranh giới ngày - đêm trùng khớp với đường sức từ, làm cho gradient độ dẫn Pedersen dốc nhất và cực đại hóa điện trường PRE. Tuy nhiên, quan trắc vệ tinh tại Trung tâm Thái Bình Dương ghi nhận cực đại ESF xuất hiện vào thời kỳ hạ chí. Luận án giải thích rằng vào các kỳ chí điểm, sự hình thành mầm mống LSWS do sóng trọng lực khí quyển (AGW) sinh ra từ hoạt động đối lưu nhiệt đới (ITCZ) và sự biến điệu của lực triều số sóng 2 đã cung cấp nguồn kích thích và gia tăng vận tốc nâng thẳng đứng, bù đắp cho sự thiếu hụt gradient độ dẫn hình học của STBA.

  4. Nguồn gốc vật lý của các dạng spread F quan trắc được tại trạm Phú Thụy (Việt Nam) là gì? Trạm Phú Thụy nằm ở khoảng $10^\circ$ vĩ độ từ, cách xa xích đạo từ. Các dạng spread F phân tầng phức hợp ghi nhận tại đây không sinh ra từ các bất ổn định trao đổi plasma tại chỗ ở đáy lớp F địa phương. Chúng có nguồn gốc từ các cấu trúc bong bóng plasma xích đạo (EPB) hình thành tại vùng xích đạo từ, sau đó phát triển vươn cao lên đỉnh lớp F và lan truyền dọc theo các đường sức từ trường Trái Đất xuống các vĩ độ thấp hơn. Tín hiệu thăm dò thẳng đứng từ trạm gặp các mặt nghiêng của bong bóng plasma này tạo thành các phản xạ xiên hiển thị trên điện ly đồ.

  5. Quy trình xác định thành phần nhiễu loạn PTEC từ vệ tinh C/NOFS trong luận án được thực hiện ra sao? Quy trình gồm hai bước: Đầu tiên, tín hiệu pha hai tần số 150 MHz và 400 MHz từ máy phát CERTO trên vệ tinh C/NOFS được máy thu GRBR mặt đất ghi nhận; tác giả áp dụng phương pháp một trạm của Yamamoto et al. cùng thuật toán khớp bình phương tối thiểu phi tuyến để tính toán giá trị nồng độ điện tử tổng cộng tuyệt đối theo phương thẳng đứng (VTEC). Tiếp đó, giá trị VTEC được làm mịn bằng hàm trung bình trượt trong khoảng thời gian 2,5 phút (tương đương quy mô không gian quỹ đạo khoảng 1000 km) để xác định biến thiên nền; hiệu số giữa VTEC thực đo và đường biến thiên nền này chính là thành phần nhiễu loạn PTEC, đại diện cho cấu trúc LSWS.

Kết luận

Luận án đã giải quyết thỏa đáng mục tiêu xác định cơ chế và các yếu tố điều khiển gây biến thiên hoạt động của hiện tượng spread F xích đạo thông qua việc phối hợp số liệu đa nguồn từ trạm mặt đất và vệ tinh chuyên dụng. Công trình khẳng định vết phản xạ nhiều lần trên điện ly đồ là dấu hiệu nhận biết gián tiếp tin cậy của cấu trúc dạng sóng quy mô lớn, đồng thời chứng minh vận tốc nâng lên cực đại của lớp F giữ vai trò chi phối quan trọng đối với biến thiên của hiện tượng ở quy mô hàng tháng. Ngoài ra, việc làm sáng tỏ nguồn gốc các cấu trúc spread F phức hợp tại trạm Phú Thụy liên quan đến sự lan truyền của bong bóng plasma từ xích đạo từ đã cung cấp luận cứ khoa học quan trọng cho việc khai thác chuỗi số liệu quan trắc điện ly tại Việt Nam.