A Singular Value Decomposition Framework for Jovian Radio Emissions from Parker Solar Probe

ملخص تقني: إطار عمل تفكيك القيم المفردة لانبعاثات راديو المشتري من مسبار باركر الشمسي

بيان المشكلة
منذ إطلاقه في عام 2018، عمل مسبار باركر الشمسي (PSP) كمراقب بعيد طويل الأمد للنظام المشتري، حيث سجل بيانات راديوية وبلازما عالية الدقة خلال 21 عملية اقتراب من الشمس. ومع ذلك، فإن استخراج انبعاثات المشتري من نوع "الديكامتر" (DAM) و"الهكتومتر" (HOM) الضعيفة من هذه المجموعة من البيانات يمثل تحديًا تحليليًا كبيرًا. فمن المعتاد أن تتعرض البيانات للتشبع بسبب ضجيج البلازما شبه الحراري، والتداخل الأدواتي للمركبة الفضائية، والانفجارات الشمسية واسعة الطيف (النوعين الثاني والثالث). إن طرق الاستخراج التقليدية، التي تعتمد على قيود مكانية قريبة أو تحديد اتجاه محلي، غير فعالة وعرضة للتحيز الذاتي عند تطبيقها على منصة غير مخصصة وبعيدة مثل PSP. علاوة على ذلك، لا تتوفر بيانات استقطاب ستوكس $V$ القياسية إلا خلال نوافذ الإرسال عالي المعدل (حوالي يومين شهريًا)، مما يحول دون المراقبة المستمرة طويلة الأمد.

المنهجية
يقدم المؤلفون خط معالجة تجريبيًا عامًا مصممًا للفصل الأعمى بين الإشارات الكوكبية المهيكلة والخلفية هليوسفيرية العشوائية. تمر المنهجية عبر أربع مراحل رئيسية:

1. طي الطور وتنسيق البيانات:
   يتم تحويل الأطياف الديناميكية الخام إلى إطار مرجعي دوري مطوي الطور بالنسبة لطول نظام المشتري الثالث ($\lambda_{III}$) وطور مدار قمر "آيو" ($\Phi_{Io}$). باستخدام تقويمات ناسا (SPICE)، يتم تصحيح الطوابع الزمنية لزمن انتقال الضوء وإسقاطها على هذه الإحداثيات الدوارة. هذا إعادة الهيكلة تحول البيانات الزمنية إلى مصفوفات دورية متكررة، مما يؤدي إلى إنشاء خط أساس لتقليل الضجيج دون اللجوء إلى المتوسط الزمني.

2. تصفية الضجيج الحتمي:

• تبييض الضجيج شبه الحراري: تقوم دالة نسبة مئوية ديناميكية للضجيج بتقييم الوسيط والنسبة المئوية الـ 75 لجانب البيانات الموجب لكل تردد على حدة. ينتج عن ذلك مضاعف يعتمد على التردد لتسطيح مستمر الضجيج الخلفي، مما يضمن عدم تصنيف حد تردد بلازما الإلكترونات المحلي خطأً كنبضات عابرة.
• تخفيف الانفجارات الشمسية: تكتشف خوارزمية تحديد العتبة الآلية الانفجارات الطيفية الواسعة (الانفجارات الشمسية) عن طريق تحديد الزيادات المستمرة في الانحراف المعياري عبر قنوات التردد بينما يظل متوسط السعة قريبًا من الصفر. يتم إزالة هذه الفترات الفاسدة واستبدالها عبر "ترميم الوسيط" (median inpainting) بمقياس يتناسب مع أرضية الضجيج المحلية.

3. استخراج تفكيك القيم المفردة (SVD):
   يتم تفكيك الأطياف المرشحة باستخدام SVD ($M = U\Sigma V^T$). ويحكم فصل الإشارة عن الضجيج توزيع القيم المفردة، التي تظهر انخفاضًا حادًا يليه "كتف ضجيج".

• اختيار الرتبة: يتم تحديد معلم رتبة حتمي $k$ من خلال الاحتفاظ بالقيم المفردة حتى نقطة الانعطاف.
• تحديد عتبة الإنتروبيا: لاختيار الأنماط بشكل موضوعي، يتم حساب إنتروبيا شانون لمتجهات القاعدة الطيفية الموحدة. تتركز انبعاثات المشتري المهيكلة في مجموعة فرعية من القنوات (إنتروبيا منخفضة)، بينما توزع الضجيج العشوائي الطاقة بشكل موحد (إنتروبيا عالية). يتم الاحتفاظ بالأنماط فقط إذا كانت إنتروبيا $S_i$ أقل من عتبة $S_{lim} = 3.0$ نات (nats).
• إعادة البناء: يتم إعادة بناء المصفوفة باستخدام الأنماط ذات الإنتروبيا المنخفضة المحتفظ بها فقط، مما يعزل فعليًا إشارة المشتري.

4. الوجوه الذاتية (Eigenfaces) للاتجاهات العالمية:
   يتم تطبيق تحليل مصفوفة ثانوي، مقتبس من تقنية "الوجوه الذاتية" المستخدمة في التعرف على الوجوه، على الأرشيف بأكمله. يتم تسطيح مصفوفات الأطياف الفردية إلى متجهات وتكديسها لاستخلاص نمط هيكلي عالمي مهيمن ("الوجه الذاتي"). يتم حساب معامل ارتباط معياري $C$ (تشابه جيب التمام) لكل طيف مقابل هذا الأساس، ليعمل كمقياس قياسي للكشف عن الاتجاهات المورفولوجية طويلة الأمد والمسح الأرشيفي.

المساهمات والنتائج الرئيسية

• تحسين نسبة الإشارة إلى الضجيج (SNR): أظهر خط المعالجة تحسينات تجريبية في SNR بمعدل $1.54 \pm 0.11$ ديسيبل، مع وصول التحسينات في الأحداث الفردية إلى $3.32$ ديسيبل.
• التحقق من فصل الإشارة:
  • الارتباط الذاتي: تظهر مصفوفة SVD المعاد بناؤها ارتباطًا ذاتيًا بفارق زمن قدره 1 ($lag-1$) بقيمة$r_{svd} = 0.797 \pm 0.003$، مما يؤكد وجود بنية مترابطة مكانيًا. وعلى العكس من ذلك، تعطي المصفوفة المتبقية$r_{res} = 0.040 \pm 0.017$، وهو ما يتوافق مع الضجيج الأبيض غير المترابط.
  • التحقق من الاستقطاب المتقاطع: تظهر بيانات الارتباط المتقاطع $C_{XY}$ المستخرجة (والتي تعمل كوكيل لستوكس $V$) ارتباط بيرسون $r = 0.783$ مع أطياف ستوكس $V$ المعالجة بشكل مستقل خلال نوافذ معدل البيانات العالي. أظهرت التجارلات الضابطة التي طبقت SVD مباشرة على بيانات ستوكس $V$ أداءً أقل بكثير ($r_{svj,V} = 0.32$ وتحسن في SNR قدره $0.64$ديسيبل)، مما يؤكد أن ميزة الاستخراج تنبع من سلسلة إشارة$C_{XY}$ المحددة وليس من خوارزمية SVD نفسها.
• الاستقلال الهندسي: نجح خط المعالجة في استعادة السمات الفيزيائية المعروفة دون فرض قيود هندسية أثناء الاستخراج. ويشمل ذلك:
  • ذروات طور "آيو" المركزة عند $\Phi_{Io} \approx 100^\circ$ و $260^\circ$.
  • منخفض انخفاض بالقرب من $\Phi_{Io} \approx 150^\circ$.
  • عدم تناظر نصف كروي منهجي في توزيعات خط الطول (ثنائي النمط للانبعاثات ذات الاستقطاب الدائري الأيسر، وأكثر تسطحًا لليمين).
• التكامل مع مسبار جونو: نجح الإطار في عزل عابرات HOM في نطاقات التردد (150 كيلوهرتز–3.5 ميجاهرتز) حيث تكون بيانات Juno-Waves ملوثة بالتداخل الأدواتي. وتتوافق عمليات الكشف المتزامنة بين PSP وJuno مع تأخير زمن انتقال الضوء المتوقع هندسيًا ($\sim 36.84$ دقيقة)، مما يوفر تأكيدًا مستقلًا للأصل المشتري للإشارات المستخرجة.

الأهمية والادعاءات
يؤسس البحث جدوى اعتبار PSP مرصدًا راديويًا للمشتري طويل الأمد وغير مخصص. ومن خلال استخدام إطار عمل رياضي آلي وقوي يعتمد على SVD وطي الطور، يوضح المؤلفون أنه يمكن استخراج إشارات كوكبية باهتة من بيئات هليوسفيرية عالية الضجيج دون الاعتماد على قيود مكانية محلية أو تنسيق يدوي.

يصيغ المؤلفون مقياس "الوجه الذاتي" بتواضع ليس كنموذج تصنيف نهائي على مستوى البكسل، بل كأداة فحص أولي احتمالية وفعالة حاسبيًا تُثري مجموعات المرشحين بنحو 30 ضعفًا. تكمن الأهمية الأساسية في توفير منهجية عامة مفتوحة المصدر يمكن تطبيقها بأثر رجعي على مجموعات بيانات هليوسفيرية أرشيفية أخرى أو دمجها في مهمات كوكبية مستقبلية لتحديد الإشارات العابرة بشكل موضوعي في بيئات ضجيج معقدة. ويؤكد العمل أن المركبات الفضائية البعيدة، عندما تُجهز بأدوات حساسة وتُعالج عبر أطر رياضية صارمة، يمكن أن تقدم بيانات علمية قيمة في علم الفلك الراديوي الكوكبي.