Luận án Tiến sĩ Vật lý: Nghiên cứu EEJ từ CHAMP và số liệu mặt đất tại Việt Nam

Tổng quan về luận án

Nghiên cứu vật lý địa cầu hiện đại đòi hỏi sự kết hợp chặt chẽ giữa quan trắc vệ tinh quỹ đạo thấp (Low Earth Orbit - LEO) và mạng lưới trạm mặt đất nhằm giải mã các quá trình điện động lực học phức tạp trong hệ thống Trái Đất - Không gian. Luận án tiến sĩ vật lý với đề tài "Nghiên cứu dòng điện xích đạo (EEJ) từ số liệu vệ tinh CHAMP và từ số liệu mặt đất ở khu vực Việt Nam và các vùng lân cận" do nghiên cứu sinh Lê Trường Thanh thực hiện (chuyên ngành Vật lý địa cầu, mã số 62 44 01 11, dưới sự hướng dẫn của PGS. Hà Duyên Châu và TS. Lê Huy Minh tại Viện Hàn lâm Khoa học và Công nghệ Việt Nam) là công trình tiên phong giải quyết bài toán mô hình hóa trường từ dòng điện xích đạo (Equatorial Electrojet - EEJ) và trường từ bình thường (TTBT) tại khu vực Đông Nam Á.

Về mặt bối cảnh khoa học, EEJ là một dải dòng điện cường độ cao chạy theo hướng đông ở độ cao khoảng 105 km thuộc lớp E của tầng điện ly trong phạm vi hẹp $\pm 5^\circ$ quanh xích đạo từ (Chapman, 1951). Dù chỉ chiếm tỷ trọng nhỏ trong tổng trường địa từ ghi nhận được, biến thiên của EEJ có thể đạt biên độ từ hàng chục đến hàng trăm nanotesla (nT), gây ảnh hưởng nghiêm trọng đến các phép đo từ tellua, khảo sát địa vật lý và định vị vệ tinh. Luận án xác định ba khoảng trống nghiên cứu then chốt (research gaps):

1. Các nghiên cứu trước đây trên thế giới sử dụng số liệu vệ tinh LEO (như Doumouya et al., 2004; Lühr et al., 2004) chỉ phân tích các chuỗi dữ liệu rất ngắn (2–3 tháng) trong giai đoạn cực đại hoạt động Mặt Trời (2001), dẫn đến hiện tượng nhiễu từ quyển mạnh và chưa phản ánh được biến thiên theo mùa cũng như theo chu kỳ hoạt động Mặt Trời.
2. Tại khu vực kinh tuyến qua Việt Nam ($100^\circ\text{E} - 105^\circ\text{E}$), mặc dù các quan sát sơ khởi gợi mở sự tồn tại của một cực đại dòng EEJ toàn cầu, nhưng chưa có công trình nào sử dụng chuỗi số liệu vệ tinh dài hạn đồng bộ với số liệu trạm mặt đất để kiểm chứng định lượng.
3. Từ năm 2003 đến thời điểm thực hiện luận án, Việt Nam hoàn toàn chưa xây dựng bất kỳ bản đồ TTBT mới nào. Các mô hình trước đây (Hà Duyên Châu, 1997, 2003; Nguyễn Thị Kim Thoa et al., 1992) chỉ sử dụng phép xấp xỉ đa thức bậc hai trên mạng lưới 56–58 điểm đo lặp mặt đất, dẫn đến sai số biên lớn và bỏ trống toàn bộ vùng lãnh hải, thềm lục địa rộng lớn của Biển Đông.

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học tường minh:

• RQ1: Cấu trúc không gian - thời gian của hệ dòng EEJ trên quy mô toàn cầu và tại kinh tuyến Việt Nam biến đổi như thế nào qua các pha của chu kỳ hoạt động Mặt Trời?
• RQ2: Làm thế nào để loại bỏ nhiễu trường ngoài và tách biệt chính xác tín hiệu từ trường EEJ từ chuỗi số liệu từ trường tổng và vector của vệ tinh CHAMP mà không làm biến dạng tín hiệu cục bộ?
• RQ3: Phương pháp phân tích điều hòa chỏm cầu (Spherical Cap Harmonic Analysis - SCHA) có thể thay thế hoàn toàn các phương pháp xấp xỉ đa thức truyền thống để xây dựng bản đồ TTBT và dị thường từ độ chính xác cao cho cả vùng đất liền và vùng biển Việt Nam hay không?
• H1: Mật độ dòng điện cực đại tại tâm EEJ ($j_0$) đạt giá trị cao nhất toàn cầu tại dải kinh tuyến Đông Nam Á ($100^\circ\text{E} - 105^\circ\text{E}$) do sự cộng hưởng giữa độ dẫn Cowling và cấu trúc hình học của đường sức từ trường chính.
• H2: Tồn tại mối tương quan tuyến tính chặt chẽ giữa mật độ dòng $j_0$ với số vết đen Mặt Trời trung bình trượt 12 tháng ($R_z$) và chu kỳ biến thiên mùa với hai đỉnh cực đại tại các phân điểm.
• H3: Việc áp dụng phương pháp SCHA với số liệu vệ tinh CHAMP ban đêm cho phép mô hình hóa trường từ chính và trường vỏ khu vực với sai số bình phương trung bình (RMS) thấp hơn đáng kể so với mô hình IGRF toàn cầu và mô hình đa thức Taylor.

Về quy mô và phạm vi dữ liệu, luận án xử lý chuỗi dữ liệu 6 năm liên tục (2002–2007) từ vệ tinh CHAMP với hơn 31 triệu điểm đo mỗi năm ở tần số lấy mẫu 1 Hz (từ kế proton OVM) và 50 Hz (từ kế fluxgate FGM), kết hợp dữ liệu đối sánh từ 6 đài địa từ quốc tế và quốc gia (Bạc Liêu - BCL, Phú Thụy - PHU, Huancayo - HUA, Fuquene - FUQ, Addis Ababa - AAE, Qsaybeh - QSB). Khung lý thuyết tích hợp thuyết dynamo tầng điện ly, mô hình dòng 3EM (3-Dimensional Empirical Model) và giải tích điều hòa chỏm cầu Haines (1985).

Literature Review và Positioning

Nghiên cứu về EEJ và TTBT đã trải qua quá trình phát triển sâu rộng trong địa vật lý quốc tế và Việt Nam. Về mặt lịch sử, hiện tượng tăng đột biến của biên độ biến thiên ngày đêm thành phần nằm ngang $H$ gần xích đạo từ được phát hiện lần đầu bởi McNish (1934) tại trạm Huancayo (Peru) và Egedal (1947). Chapman (1951) đã giải thích cơ chế vật lý này bằng sự tồn tại của hệ dòng điện chảy về phía đông ở độ cao xích đạo, sinh ra bởi sự bất đồng nhất độ dẫn điện ly dưới tác động của bức xạ ion hóa Mặt Trời.

                    TIẾN TRÌNH NGHIÊN CỨU DÒNG ĐIỆN XÍCH ĐẠO (EEJ)
┌──────────────────────────────────────────────────────────────────────────────────┐
│ Quan sát mặt đất & Tên lửa (1934 - 1980)                                         │
│ • McNish (1934), Egedal (1947), Chapman (1951): Phát hiện EEJ & cơ chế độ dẫn    │
│ • Fambitakoye & Mayaud (1976): Mô hình dải dòng mặt đất Châu Phi (a, j0)         │
│ • Sampath et al. (1979): Đo đạc tên lửa tại Thumba xác định đỉnh dòng ở 106 km   │
└────────────────────────────────────────┬─────────────────────────────────────────┘
                                         │
                                         ▼
┌──────────────────────────────────────────────────────────────────────────────────┐
│ Vệ tinh giai đoạn đầu (1970 - 2000)                                              │
│ • POGO (Cain & Sweeney, 1973; Onwumechili, 1980): Phát hiện đỉnh EEJ tại 100°E   │
│ • MAGSAT (Cohen, 1985; Langel, 1985): Phát hiện dòng ngược ban đêm/sáng sớm      │
│ • Ørsted (Jadhav et al., 2002): Đỉnh dòng tại 135°E nhưng quỹ đạo bình minh/hoàng hôn │
└────────────────────────────────────────┬─────────────────────────────────────────┘
                                         │
                                         ▼
┌──────────────────────────────────────────────────────────────────────────────────┐
│ Kỷ nguyên Vệ tinh Quỹ đạo Thấp CHAMP (2000 - 2010)                               │
│ • Doumouya et al. (2004): 4 đỉnh EEJ toàn cầu (chuỗi 2 tháng, solar max)        │
│ • Lühr et al. (2004): Mật độ dòng trung bình 150 A/km (chưa xét biến thiên mùa)  │
│ • Alken & Maus (2007): Đối chiếu EEJ với mô hình gió GSWM-02                     │
└────────────────────────────────────────┬─────────────────────────────────────────┘
                                         │
                                         ▼
┌──────────────────────────────────────────────────────────────────────────────────┐
│ Đóng góp đột phá của Luận án (Lê Trường Thanh, 2015)                             │
│ ➔ Chuỗi dữ liệu đồng bộ 6 năm (2002-2007) CHAMP + 6 đài địa từ toàn cầu          │
│ ➔ Lọc đa thức thích nghi (Adaptive Polynomials bậc 6-18) tách chính xác ΔF       │
│ ➔ Xác nhận đỉnh cực đại j0 tại kinh tuyến Việt Nam (100°E-105°E) xuyên suốt chu kỳ│
│ ➔ Lần đầu ứng dụng SCHA xây dựng bản đồ TTBT 2007.0 phủ kín đất liền & Biển Đông │
└──────────────────────────────────────────────────────────────────────────────────┘

Trong dòng nghiên cứu quan sát mặt đất, Fambitakoye và Mayaud (1976) đã đặt nền móng định lượng khi sử dụng mạng lưới 9 trạm dọc kinh tuyến $17^\circ\text{E}$ (Châu Phi) để xác định nửa bề rộng dòng điện ($a$) đạt $486\text{ km}$ vào giữa trưa và thu hẹp còn $291\text{ km}$ lúc chiều tối, với mật độ dòng điện tại tâm ($j_0$) đạt cực đại $145\text{ A/km}$. Doumouya et al. (1998) tiếp tục chuẩn hóa mô hình tại Tây Phi ($5^\circ\text{W}$), ghi nhận $j_0 = 232\text{ A/km}$ và $a = 303 \pm 29\text{ km}$. Tuy nhiên, các công trình mặt đất thuần túy bộc lộ hạn chế cố hữu: không thể bao quát các vùng đại dương và phân bố trạm không đối xứng qua xích đạo từ.

Sự xuất hiện của các vệ tinh địa từ đã mở ra kỷ nguyên mới. Chuỗi vệ tinh POGO (Cain & Sweeney, 1973; Onwumechili & Agu, 1980) lần đầu tiên cung cấp dữ liệu toàn cầu nhưng chịu hạn chế do quỹ đạo bay quá cao và độ chính xác thiết bị thấp. Vệ tinh MAGSAT (Cohen, 1985; Langel et al., 1985) và Ørsted (Jadhav et al., 2002) phát hiện hiện tượng dòng điện ngược xích đạo (Counter Electrojet - CEJ), song Ørsted bị giới hạn bởi quỹ đạo đồng bộ Mặt Trời lúc bình minh/hoàng hôn.

Khi vệ tinh CHAMP đi vào hoạt động với quỹ đạo bay thấp ($300 - 454\text{ km}$), Doumouya et al. (2004) phát hiện cấu trúc 4 đỉnh cực đại của EEJ trên toàn cầu, trong đó đỉnh tại kinh tuyến $105^\circ\text{E}$ đạt giá trị lớn nhất. Lühr et al. (2004) ghi nhận mật độ dòng trung bình toàn cầu khoảng $150\text{ A/km}$ nhưng chưa giải quyết được biến thiên theo mùa. Alken và Maus (2007) so sánh EEJ từ CHAMP với mô hình vận tốc gió GSWM-02 nhưng ghi nhận sự sai lệch lớn, chứng minh rằng động lực học gió đơn thuần không thể giải thích toàn bộ sự biến thiên theo kinh độ của EEJ.

Tại Việt Nam, các nghiên cứu địa từ ban đầu của Đặng Văn Hưng (1985), Trương Quang Hảo và Lê Huy Minh (1987) đã xác định tốc độ dịch chuyển của xích đạo từ về phía nam là khoảng $4\text{ km/năm}$. Nguyễn Văn Giảng (1988) dùng số liệu MAGSAT phân tích điều hòa cầu (Spherical Harmonic Analysis - SHA) bậc $n=13$ nhưng sai số mô hình còn lớn ($\pm 60\text{ nT}$). Nguyễn Thị Kim Thoa et al. (1990, 1992) và Hà Duyên Châu (1997, 2003) xây dựng TTBT bằng đa thức bậc hai cho các niên đại 1991.5, 1997.5 và 2003.5 nhưng chỉ giới hạn trong phạm vi đất liền.

Trong bối cảnh quốc tế về mô hình hóa khu vực, Haines (1985) phát triển phương pháp SCHA, mở ra khả năng mô hình hóa trường từ trên chỏm cầu với số liệu ở các độ cao khác nhau. Phương pháp này được áp dụng thành công tại Canada (Haines, 1985), Ý (De Santis, 1990; Duka et al., 1998), Nam Mỹ (Kotzé & Haak, 2002), Trung Quốc (Gu et al., 2006) và Philippines (Zhi et al., 2011). Luận án của Lê Trường Thanh định vị chính xác tại điểm giao thoa: lấp đầy khoảng trống dữ liệu dài hạn về EEJ tại Đông Nam Á bằng chuỗi dữ liệu CHAMP 6 năm, đồng thời là công trình đầu tiên ứng dụng SCHA tại Việt Nam để giải bài toán TTBT cho cả vùng đất liền và vùng biển.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng và kiểm chứng thực nghiệm sâu sắc thuyết Dynamo tầng điện ly lớp E và cơ chế độ dẫn Cowling do Chapman (1951) và Richmond (1973) thiết lập. Theo lý thuyết cổ điển, độ dẫn song song ($\sigma_\parallel$), độ dẫn Pedersen ($\sigma_P$) và độ dẫn Hall ($\sigma_H$) trong tầng điện ly biến đổi mạnh theo độ cao:

$$\sigma_\parallel = \frac{n_e e^2}{m_e \nu_{en}}, \quad \sigma_P = \frac{n_e e}{B} \left( \frac{\kappa_i}{1 + \kappa_i^2} + \frac{\kappa_e}{1 + \kappa_e^2} \right), \quad \sigma_H = \frac{n_e e}{B} \left( \frac{1}{1 + \kappa_e^2} - \frac{1}{1 + \kappa_i^2} \right)$$

Tại mặt phẳng xích đạo từ, nơi đường sức từ trường Trái Đất nằm ngang hoàn toàn ($I = 0^\circ$), sự tích tụ điện tích ở ranh giới trên và dưới của lớp E tạo ra điện trường phân cực thẳng đứng $E_z = (\sigma_H / \sigma_P) E_y$. Điện trường này khuếch đại dòng điện Hall theo phương ngang, làm xuất hiện độ dẫn tổng hợp Cowling:

$$\sigma_C = \sigma_P + \frac{\sigma_H^2}{\sigma_P}$$

                CƠ CHẾ VẬT LÝ HÌNH THÀNH ĐỘ DẪN COWLING VÀ EEJ
┌─────────────────────────────────────────────────────────────────────────────┐
│ 1. Bức xạ ion hóa Mặt Trời tác động lên tầng điện ly lớp E (độ cao ~105 km) │
│    ➔ Tạo nồng độ ion và điện tử ne                                         │
└──────────────────────────────────────┬──────────────────────────────────────┘
                                       │
                                       ▼
┌─────────────────────────────────────────────────────────────────────────────┐
│ 2. Điện trường Đông - Tây ban ngày Ey kết hợp Từ trường chính ngang B       │
│    ➔ Dòng Hall thẳng đứng làm tích tụ điện tích tại biên lớp E             │
└──────────────────────────────────────┬──────────────────────────────────────┘
                                       │
                                       ▼
┌─────────────────────────────────────────────────────────────────────────────┐
│ 3. Thiết lập Điện trường phân cực thẳng đứng: Ez = (σH / σP) * Ey           │
│    ➔ Khuếch đại mạnh mẽ dòng hạt tải điện theo hướng Đông                   │
└──────────────────────────────────────┬──────────────────────────────────────┘
                                       │
                                       ▼
┌─────────────────────────────────────────────────────────────────────────────┐
│ 4. Xuất hiện Độ dẫn Cowling: σC = σP + (σH^2 / σP)                          │
│    ➔ Hình thành Dải dòng điện xích đạo (EEJ) với mật độ j0 cực đại          │
└─────────────────────────────────────────────────────────────────────────────┘

Đóng góp lý thuyết của luận án thể hiện qua 3 mệnh đề khoa học (Propositions):

• Mệnh đề P1: Độ lớn của độ dẫn Cowling $\sigma_C$ và cường độ EEJ không thuần nhất theo kinh độ mà chịu sự điều chế sâu sắc bởi độ lớn của từ trường chính Trái Đất $B_0$. Tại kinh tuyến Đông Nam Á ($105^\circ\text{E}$), do giá trị $B_0$ nhỏ hơn so với nhiều khu vực cùng vĩ độ, tỷ số $\sigma_C \propto 1/B_0$ tăng cao, dẫn đến mật độ dòng $j_0$ đạt cực đại toàn cầu.
• Mệnh đề P2: Cấu trúc không gian của dải dòng EEJ tuân theo phân bố parabol bậc hai mở rộng $j(x) = j_0 [1 - (x/a)^2]^m$ với $m=2$, trong đó bề rộng dải dòng tỷ lệ nghịch với độ dốc của từ khuynh $\partial I / \partial \phi$.
• Mệnh đề P3: Biến thiên theo mùa của EEJ chịu sự kiểm soát đồng thời của góc phương vị Mặt Trời và hệ thống sóng triều khí quyển di chuyển (diurnal migrating tides), tạo nên hai đỉnh đối xứng tại phân điểm xuân và phân điểm thu.

Luận án đánh dấu bước chuyển biến nhận thức luận (epistemological shift): từ mô hình quan sát tĩnh cục bộ tại từng trạm mặt đất đơn lẻ chuyển sang mô hình động lực học không gian đa chiều, liên kết tương tác trường ngoài (tầng điện ly) và trường trong (vỏ và nhân Trái Đất).

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đồng thời ba cấu trúc lý thuyết và phương pháp toán học phức hợp:

1. Mô hình thực nghiệm ba chiều 3EM (3-Dimensional Empirical Model): Mô hình hóa mật độ dòng diện tích $j(x, \lambda, t) = j(x) \cdot G(t) \cdot f(\lambda)$, kết hợp hàm phân bố vĩ độ $j(x)$, hàm phân bố thời gian địa phương $G(t) = \cos^d \left( \frac{\pi (t - t_m)}{T} \right)$ và hàm điều chế kinh độ $f(\lambda)$.
2. Phương pháp phân tích điều hòa chỏm cầu (SCHA): Giải phương trình Laplace cho thế từ vị $\nabla^2 V = 0$ trên miền chỏm cầu giới hạn bởi góc mở $\theta_0$:

$$V(r, \theta, \lambda) = a \sum_{k=0}^{K_{max}} \sum_{m=0}^{k} \left( \frac{a}{r} \right)^{n_k(m) + 1} P_{n_k(m)}^m (\cos \theta) \left( g_k^m \cos m\lambda + h_k^m \sin m\lambda \right)$$

Điểm đột phá là bậc $n_k(m)$ của đa thức Legendre liên kết $P_{n_k(m)}^m (\cos \theta)$ là các số thực không nguyên, được xác định từ điều kiện biên tại biên chỏm cầu $\theta = \theta_0$: hoặc $P_{n_k(m)}^m (\cos \theta_0) = 0$ (cho đạo hàm tiếp tuyến), hoặc $dP_{n_k(m)}^m (\cos \theta)/d\theta |_{\theta_0} = 0$ (cho đạo hàm pháp tuyến). 3. Kỹ thuật nghịch đảo ma trận thích nghi: Tách biệt thành phần trường nội và trường ngoại thông qua hệ số Gauss mở rộng trên chỏm cầu.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án được thiết kế theo trường phái nhận thức thực chứng định lượng (quantitative positivism), sử dụng quy trình xử lý dữ liệu thực nghiệm quy mô lớn kết hợp mô phỏng giải tích. Thiết kế nghiên cứu đa tầng (multi-level design) bao gồm:

• Cấp độ 1 (Toàn cầu): Xử lý số liệu quỹ đạo cực của vệ tinh CHAMP trên toàn bộ 360 kinh tuyến nhằm xác lập bức tranh tổng thể về EEJ và phát hiện các dị thường kinh độ.
• Cấp độ 2 (Khu vực): Ứng dụng SCHA trên chỏm cầu có góc mở $\theta_0 = 20^\circ$ (tâm chỏm cầu tại $\phi_0 = 16^\circ\text{N}, \lambda_0 = 106^\circ\text{E}$) bao trùm toàn bộ lãnh thổ Việt Nam và Biển Đông.
• Cấp độ 3 (Điểm trạm chuẩn): Kiểm chứng mô hình đối sánh với 6 đài địa từ mặt đất đạt chuẩn INTERMAGNET và MAGDAS.

┌─────────────────────────────────────────────────────────────────────────────┐
│                    DỮ LIỆU ĐẦU VÀO VỆ TINH VÀ MẶT ĐẤT                       │
│ • Vệ tinh CHAMP Mức 3 (2002-2007): Từ kế OVM (1 Hz) & FGM (50 Hz), NEC      │
│ • 6 Đài địa từ: BCL, PHU (Việt Nam); HUA, FUQ (Mỹ); AAE, QSB (Châu Phi)     │
└──────────────────────────────────────┬──────────────────────────────────────┘
                                       │
                   ┌───────────────────┴───────────────────┐
                   ▼                                       ▼
┌──────────────────────────────────────┐┌──────────────────────────────────────┐
│  DỮ LIỆU BAN NGÀY (10:00 - 14:00 LT) ││  DỮ LIỆU BAN ĐÊM (22:00 - 04:00 LT)  │
│  Tách và mô hình hóa dòng EEJ        ││  Xây dựng Từ trường bình thường SCHA │
└──────────────────┬───────────────────┘└──────────────────┬───────────────────┘
                   │                                       │
                   ▼                                       ▼
┌──────────────────────────────────────┐┌──────────────────────────────────────┐
│ 1. Trừ trường chính POMME 3.0 / IGRF ││ 1. Lọc ngày yên tĩnh từ (Kp <= 2+)   │
│ 2. Lọc đa thức thích nghi (bậc 6-18) ││ 2. Thiết lập chỏm cầu: θ0=20°, Kmax=8│
│ 3. Tính mật độ dòng cực đại j0       ││ 3. Khai triển Legendre bậc phi nguyên│
│ 4. Mô hình hóa thực nghiệm 3EM       ││ 4. Nghịch đảo bình phương tối thiểu  │
└──────────────────┬───────────────────┘└──────────────────┬───────────────────┘
                   │                                       │
                   └───────────────────┬───────────────────┘
                                       ▼
┌─────────────────────────────────────────────────────────────────────────────┐
│                    KẾT QUẢ ĐẦU RA ĐỘ PHÂN GIẢI CAO                          │
│ • Quy luật không gian - thời gian của EEJ (Đỉnh cực đại tại 100°E - 105°E)  │
│ • Bộ bản đồ TTBT 7 thành phần (X, Y, Z, F, H, D, I) niên đại 2007.0         │
│ • Bản đồ dị thường từ vỏ Trái Đất (ΔXa, ΔYa, ΔZa, ΔFa) khu vực Biển Đông   │
└─────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình xử lý dữ liệu tuân thủ các bước nghiêm ngặt:

1. Lựa chọn và tiền xử lý dữ liệu vệ tinh CHAMP: Luận án sử dụng dữ liệu Mức 3 (Level 3 data) lưu trữ dưới định dạng CDF (Common Data Format). Các lát cắt quỹ đạo được chọn lọc theo điều kiện từ trường yên tĩnh ($Kp \le 2+$ hoặc $Kp \le 3$).
   • Dữ liệu ban ngày (10h00 đến 14h00 LT) dùng để trích xuất tín hiệu EEJ.
   • Dữ liệu ban đêm (22h00 đến 04h00 LT) khi tầng điện ly không còn bị ion hóa được dùng để tính TTBT nhằm loại bỏ hoàn toàn ảnh hưởng của dòng EEJ và dòng $Sq$.
2. Thuật toán tách trường EEJ bằng đa thức thích nghi: Trừ trường từ chính bên trong Trái Đất bằng mô hình POMME 3.0 để thu được trường dư $F_{res}$. Để loại bỏ hiệu ứng trường ngoài quy mô lớn mà không làm triệt tiêu tín hiệu EEJ cục bộ, luận án áp dụng phép lọc đa thức với bậc biến thiên từ 6 đến 18 tùy theo cấu hình lát cắt, khắc phục triệt để nhược điểm của việc cố định bậc đa thức (như bậc 12 của Doumouya et al., 2004).
3. Hiệu chuẩn và loại bỏ hiệu ứng Backus: Backus (1970) đã chỉ ra rằng việc chỉ sử dụng cường độ trường tổng $F$ để nghịch đảo sẽ gây ra sai số lớn ở vùng vĩ độ thấp ($I < 20^\circ$). Luận án sử dụng trọn vẹn cả 3 thành phần vector $X, Y, Z$ trong hệ tọa độ NEC (North-East-Center) kết hợp với $F$ để triệt tiêu hoàn toàn hiệu ứng này.

Data và phân tích

Toàn bộ hệ thống thuật toán nghịch đảo ma trận SCHA, khai triển đa thức Legendre liên kết bậc phi nguyên, chuyển đổi tọa độ cầu - đề các và lọc thống kê được lập trình trực tiếp bằng ngôn ngữ khoa học Visual Fortran 5.0. Đồ họa và trực quan hóa bản đồ đồng mức được thực hiện qua phần mềm mã nguồn mở Gnuplot 4.0 kết hợp công cụ trích xuất dữ liệu của NASA.

Đánh giá độ tin cậy và sai số mô hình được kiểm soát thông qua chỉ số sai số bình phương trung bình:

$$RMS = \sqrt{\frac{1}{N} \sum_{i=1}^N (F_{obs, i} - F_{model, i})^2}$$

Trong quá trình tối ưu hóa bậc khai triển SCHA ($K_{int}$), luận án khảo sát $K_{int}$ từ 1 đến 10 và xác định $K_{int} = 8$ là bậc tối ưu, đảm bảo giá trị RMS của các thành phần $X, Y, Z$ đạt mức cực tiểu ổn định ($RMS_X \approx 12.4\text{ nT}, RMS_Y \approx 11.8\text{ nT}, RMS_Z \approx 14.2\text{ nT}$), tránh hiện tượng phân kỳ do bậc quá cao (over-fitting).

Phát hiện đột phá và implications

Những phát hiện then chốt

1. Xác nhận sự tồn tại của đỉnh cực đại EEJ toàn cầu tại kinh tuyến Việt Nam: Phân tích 6 năm số liệu CHAMP (2002–2007) khẳng định mật độ dòng điện tại tâm EEJ ($j_0$) trên toàn cầu có cấu trúc 4 đỉnh cực đại dạng sóng 4 (wave-4 structure). Trong đó, đỉnh tại kinh tuyến qua Việt Nam ($100^\circ\text{E} - 105^\circ\text{E}$) luôn đạt giá trị cao nhất thế giới trong suốt 6 năm, dao động từ $150\text{ A/km}$ đến $232\text{ A/km}$, vượt trội so với đỉnh tại Nam Mỹ ($60^\circ\text{W}$) và Châu Phi ($10^\circ\text{E}$–$40^\circ\text{E}$).
2. Định lượng quy luật biến thiên theo mùa và hoạt động Mặt Trời: Mật độ dòng $j_0$ đạt cực đại vào các tháng phân điểm (Equinoxes: tháng 3–4 và tháng 9–10) và đạt cực tiểu vào mùa đông (Solstices). Luận án phát hiện mối tương quan đồng biến chặt chẽ giữa $j_0$ với số vết đen Mặt Trời trung bình tháng $R_z$: trong giai đoạn 2002 (hoạt động Mặt Trời cao, $R_z \approx 110$), $j_0$ tại trạm Bạc Liêu đạt trên $160\text{ A/km}$; đến năm 2007 (hoạt động Mặt Trời tối thiểu, $R_z \approx 15$), $j_0$ suy giảm xuống còn khoảng $85\text{ A/km}$.

3. Cấu trúc dải dòng và sự trôi dạt của xích đạo từ: Nửa bề rộng dòng điện $a$ được xác định chính xác nằm trong khoảng $235 - 413\text{ km}$, tâm dòng điện nằm ở độ cao $105 \pm 5\text{ km}$. Tại khu vực Nam Việt Nam, vị trí thành phần thẳng đứng $Z = 0$ xác nhận xích đạo từ tiếp tục dịch chuyển về phía Nam với tốc độ trung bình $4\text{ km/năm}$, hiện đã nằm cách mũi Cà Mau chưa đầy $100\text{ km}$.
4. Xây dựng thành công bộ bản đồ TTBT 7 thành phần niên đại 2007.0 bằng SCHA: Lần đầu tiên thiết lập tập bản đồ số độ chính xác cao cho 7 thành phần địa từ ($X, Y, Z, F, H, D, I$) bao phủ toàn vẹn lãnh thổ và lãnh hải Việt Nam. Độ lệch giữa mô hình SCHA và số liệu quan trắc vệ tinh ban đêm đạt mức tuyệt đối nhỏ ($\Delta F_{\text{ĐL}} < 15\text{ nT}$), thể hiện ưu thế vượt trội so với mô hình IGRF toàn cầu (vốn có sai số lên tới $40 - 60\text{ nT}$ tại vùng rìa).
5. Bản đồ dị thường từ lớp vỏ Biển Đông: Bằng cách loại trừ TTBT mô hình hóa từ SCHA, luận án trích xuất thành công các dị thường từ cục bộ ($\Delta X_a, \Delta Y_a, \Delta Z_a, \Delta F_a$). Các dải dị thường từ phản ánh rõ nét ranh giới phân định mảng kiến tạo, đứt gãy sâu kinh tuyến $109^\circ\text{E}$ và các bồn trũng trầm tích Đệ Tam ở thềm lục địa Biển Đông.

Implications đa chiều

• Về mặt lý thuyết: Cung cấp bằng chứng thực nghiệm củng cố mô hình tương tác sóng triều khí quyển - điện ly tầng cao, giải quyết tranh luận khoa học về nguyên nhân tạo nên cấu trúc 4 đỉnh của EEJ.
• Về phương pháp luận: Chuẩn hóa quy trình xử lý số liệu vệ tinh LEO kết hợp trạm mặt đất, thiết lập tiền đề phương pháp luận để khai thác các thế hệ vệ tinh tiếp theo như SWARM (ESA) và CSES.
• Về thực tiễn và chính sách: Cung cấp cơ sở dữ liệu trường từ chuẩn quốc gia phục vụ hiệu chỉnh sai số từ biến cho các công tác thăm dò dầu khí, khoáng sản biển sâu, an ninh quốc phòng (dẫn đường tàu ngầm, hàng không) và dự báo thời tiết không gian (Space Weather).

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, luận án cũng thẳng thắn chỉ ra các giới hạn nội tại:

1. Sự suy giảm độ cao quỹ đạo của vệ tinh CHAMP: Trong giai đoạn 2001–2007, độ cao của CHAMP suy giảm liên tục từ $454\text{ km}$ xuống $348\text{ km}$. Mặc dù đã sử dụng mô hình chuẩn hóa POMME 3.0, sự thay đổi độ cao này vẫn tạo ra độ dốc nhất định trong chuỗi số liệu theo thời gian.
2. Mạng lưới trạm mặt đất khu vực còn thưa thớt: Tại Đông Nam Á, nghiên cứu chủ yếu dựa vào hai đài từ Bạc Liêu và Phú Thụy. Việc thiếu các đài trạm quan trắc liên tục ở phía nam đường xích đạo từ trong cùng kinh tuyến làm hạn chế khả năng kiểm chứng tính bất đối xứng bắc - nam của dải dòng.
3. Chưa tích hợp số liệu gió trung hòa thực đo: Việc giải thích cơ chế 4 đỉnh kinh độ mới chỉ dừng ở mức đối sánh với mô hình lý thuyết gió GSWM-02 mà chưa có số liệu quan trắc gió tầng cao trực tiếp tại trạm xích đạo.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) bao gồm:

• Hướng 1: Khai thác chuỗi dữ liệu đa vệ tinh của chùm vệ tinh SWARM (gồm 3 vệ tinh Alpha, Bravo, Charlie) nhằm nâng cao độ phân giải không gian và tách biệt tức thời biến thiên theo thời gian thực.
• Hướng 2: Phát triển phương pháp SCHA phụ thuộc thời gian (Time-dependent SCHA - T-SCHA) để mô hình hóa liên tục biến thiên thế kỷ của từ trường khu vực.
• Hướng 3: Nghiên cứu hiện tượng dòng điện ngược xích đạo (CEJ) và hiện tượng F-rải rác (Equatorial Spread-F) bằng hệ thống thu sóng vô tuyến GRBR (GNU Radio Beacon Receiver) kết hợp quan sát địa từ.

Tác động và ảnh hưởng

Luận án của Lê Trường Thanh tạo ra những tác động lan tỏa sâu rộng trên nhiều phương diện:

• Tác động học thuật: Đặt nền tảng cho trường phái nghiên cứu địa từ vệ tinh tại Việt Nam; các bài báo trích xuất từ luận án đã được công bố trên các tạp chí chuyên ngành uy tín như Tạp chí Các khoa học về Trái Đất (2010, 2011) và các kỷ yếu hội nghị quốc tế của Hiệp hội Địa từ và Khí quyển Quốc tế (IAGA, 2009), Liên đoàn Trắc địa và Địa vật lý Quốc tế (IUGG, 2011).
• Chuyển đổi kỹ thuật công nghiệp: Cung cấp lưới dữ liệu TTBT số hóa phục vụ trực tiếp cho Tổng cục Địa chất và Khoáng sản, Tập đoàn Dầu khí Việt Nam (PVN) trong việc xử lý tài liệu từ hàng không và từ biển, nâng cao độ chính xác khi xác định ranh giới các khối magma và cấu trúc móng chứa dầu.
• Tác động an ninh - chủ quyền biển đảo: Lần đầu tiên cung cấp hệ thống bản đồ 7 yếu tố địa từ bao phủ toàn bộ vùng đặc quyền kinh tế và thềm lục địa Việt Nam, phục vụ công tác dẫn đường hàng hải và các nhiệm vụ quốc phòng trên biển.

Đối tượng hưởng lợi

• Nghiên cứu sinh và Giới học thuật Địa vật lý: Tiếp cận phương pháp toán học SCHA tiên tiến, thuật toán lọc đa thức thích nghi và cơ sở dữ liệu đối sánh vệ tinh - mặt đất chuẩn mực.
• Các đơn vị Thăm dò Khoáng sản & Địa chất Biển: Sử dụng mô hình TTBT niên đại 2007.0 làm trường chuẩn để bóc tách chính xác các dị thường từ quặng hóa và cấu trúc địa chất dầu khí.
• Cơ quan Quản lý Không gian và Định vị Dẫn đường: Ứng dụng mô hình EEJ để hiệu chỉnh tầng điện ly cho các trạm định vị vi sai DGPS, GNSS thời gian thực và quản lý liên lạc vô tuyến sóng ngắn.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì? Trả lời: Luận án đã mở rộng thuyết Dynamo và độ dẫn Cowling bằng việc chứng minh và định lượng hóa quy luật: cường độ EEJ tỷ lệ nghịch với độ lớn từ trường chính $B_0$ dọc theo kinh tuyến từ, thiết lập cơ sở vật lý giải thích tại sao khu vực Đông Nam Á ($100^\circ\text{E} - 105^\circ\text{E}$) có mật độ dòng $j_0$ đạt đỉnh cực đại toàn cầu.
2. Đột phá phương pháp luận so với các nghiên cứu quốc tế trước đây? Trả lời: So với nghiên cứu của Doumouya et al. (2004) dùng bậc đa thức cố định 12 trên chuỗi 2 tháng và Jadhav et al. (2002) dùng vệ tinh Ørsted ở quỹ đạo bình minh/hoàng hôn, luận án đã phát triển thuật toán lọc đa thức thích nghi bậc 6–18 trên chuỗi 6 năm CHAMP, kết hợp lần đầu tiên thuật toán SCHA ($K_{max}=8, \theta_0=20^\circ$) để mô hình hóa toàn diện cả đất liền và lãnh hải.
3. Phát hiện gây bất ngờ nhất có số liệu thực nghiệm chứng minh? Trả lời: Đó là sự tương quan đồng pha tuyệt đối giữa mật độ dòng $j_0$ và chu kỳ hoạt động Mặt Trời kéo dài 6 năm (2002–2007) xuyên suốt từ cực đại xuống cực tiểu, với giá trị $j_0$ tại trạm Bạc Liêu suy giảm tuyến tính từ $162.4\text{ A/km}$ (năm 2002) xuống $86.1\text{ A/km}$ (năm 2007), sai lệch so với quan trắc vệ tinh luôn duy trì dưới $5\text{ A/km}$.
4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không? Trả lời: Có. Toàn bộ cấu trúc thuật toán giải đa thức Legendre phi nguyên, quy chuẩn xử lý file CDF Mức 3 của NASA, mã nguồn chương trình Visual Fortran 5.0 và các tham số điều kiện biên chỏm cầu đều được mô tả chi tiết, cho phép tái lập hoàn toàn quy trình trên các tập dữ liệu vệ tinh mới.
5. Định hướng nghiên cứu 10 năm tiếp theo được vạch ra như thế nào? Trả lời: Luận án định hình lộ trình phát triển mô hình T-SCHA 4 chiều không gian - thời gian, tích hợp dữ liệu chùm vệ tinh SWARM để theo dõi liên tục biến thiên thế kỷ và các cấu trúc dị thường sâu trong lòng vỏ Trái Đất khu vực Biển Đông.

Kết luận

1. Khẳng định vững chắc bằng thực nghiệm chuỗi 6 năm vệ tinh CHAMP rằng mật độ dòng điện xích đạo EEJ tại dải kinh tuyến Đông Nam Á ($100^\circ\text{E} - 105^\circ\text{E}$) đạt giá trị cực đại toàn cầu ($150 - 232\text{ A/km}$).
2. Phát triển thành công thuật toán tách trường EEJ thích nghi bằng đa thức bậc biến thiên (bậc 6 đến 18), loại bỏ triệt để hiện tượng méo tín hiệu trường cục bộ.
3. Định lượng hóa toàn diện quy luật biến thiên không gian - thời gian của EEJ: chu kỳ ngày đêm cực đại lúc 11h30–12h30 LT, chu kỳ mùa cực đại tại các phân điểm và tương quan tuyến tính chặt chẽ với số vết đen Mặt Trời $R_z$.
4. Lần đầu tiên tại Việt Nam ứng dụng thành công phương pháp phân tích điều hòa chỏm cầu (SCHA) với $K_{int}=8, \theta_0=20^\circ$, giải quyết trọn vẹn bài toán mô hình hóa TTBT cho khu vực đất liền và toàn bộ Biển Đông niên đại 2007.0 với độ chính xác cao ($RMS < 15\text{ nT}$).
5. Thành lập tập bản đồ dị thường từ vỏ Trái Đất khu vực Việt Nam và lân cận, làm sáng tỏ các cấu trúc kiến tạo địa chất và ranh giới bồn trũng trầm tích thềm lục địa.
6. Mở ra hướng nghiên cứu mới trong việc kết hợp dữ liệu vệ tinh LEO hiện đại với mạng lưới trạm mặt đất, đặt nền móng kỹ thuật vững chắc cho ngành vật lý địa cầu và thăm dò tài nguyên biển Việt Nam trong kỷ nguyên không gian.