Ionospheric and Plasmaspheric Delay Characterization for Lunar Terrestrial GNSS Receivers with Global Core Plasma Model
تُوصِّف هذه الورقة التأخيرات الأيونوسفيرية والبلزموسفيرية لمستقبلات نظام الملاحة العالمي عبر الأقمار الصناعية (GNSS) القمرية من خلال استخدام نموذج البلازما الأساسي العالمي وخوارزمية تتبع الأشعة المخصصة لتحديد مقادير التأخير تحت ظروف شمسية وجيومغناطيسية متباينة، مما يساهم في تصميم خوارزميات تحديد المواقع والتوقيت القمرية القوية.
على مدى عقود، ظل القمر مكاناً صامتاً لأقمارنا الصناعية المخصصة للملاحة العالمية. فأنظمة مثل (GPS) و(Galileo) مصممة لبث إشارات نحو الأرض، حيث تصل إلى هواتفنا أو لوحات قيادة سياراتنا. ولكن مع توجه البشرية للعودة إلى القمر وبناء وجود دائم هناك، أدرك المهندسون أن هذه الأقمار الصناعية نفسها يمكن رؤيتها بالفعل من المدار القمري. إن الإشارات خافتة للغاية، تتدفق من حافة الكوكب مثل الضوء المنبعث من منارة يراها سفينة بعيدة، لكنها قوية بما يكفي لالتقاطها بواسطة أجهزة استقبال حساسة. وهذا يفتح الباب لاستخدام شبكة الملاحة الأرضية لتوجيه المركبات الفضائية ورواد الفضاء على سطح القمر، وهو مفهوم انتقل من حيز النظرية إلى الواقع مع الاختبارات الناجحة الأخيرة. ومع ذلك، هناك عقبة؛ فبينما تسافر هذه الإشارات الضعيفة مسافة 401,000 كيلومتر من الأرض إلى القمر، فإنها لا تسافر عبر فراغ مثالي. بل يجب أن تمر عبر طبقتين شاسعتين وغير مرئيتين من الغاز المشحون المحيط بكوكبنا: "الأيونوسفير" (الغلاف الأيوني)، الذي يقع بالقرب من الأرض، و"البلازموسفير" (الغلاف البلازمي)، وهو منطقة أكبر بكثير وأكثر تشتتاً تمتد بعيداً في الفضاء. تعمل هذه الطبقات مثل ضباب كثيف ومتحرك يبطئ موجات الراديو ويحرف مسارها، مما يخلق أخطاءً في حسابات التوقيت والمسافة الضرورية للملاحة.
قام فريق من الباحثين في جامعة ستانفورد الآن برسم خريطة توضح بدقة مدى تشويه هذه الطبقات غير المرئية للإشارات المتجهة إلى القمر. فقد بنوا محاكاة حاسوبية مفصلة لتتبع رحلة هذه الموجات الراديوية الخافتة أثناء انتقالها من أقمار (GPS) و(Galileo) الصناعية، عبر الغازات المشحونة في الغلاف الجوي العلوي للأرض، وصولاً إلى أجهزة الاستقبال في المدار القمري أو على السطح القمري. وبدلاً من افتراض أن الإشارات تسافر في خط مستقيم، استخدم الباحثون خوارزمية مخصصة لتتبع المسارات الفعلية المنحنية قليلاً التي تتخذها الموجات وهي تُدفع وتُسحب بفعل كثافة الإلكترونات في الفضاء المحيط بالأرض. وقد اختبروا هذه المسارات تحت مجموعة متنوعة من الظروف، حيث قاموا بمحاكاة أوقات النشاط الشمسي العالي عندما تقصف الشمس الأرض بالطاقة، وأوقات النشاط المنخفض، بالإضافة إلى مستويات مختلفة من العواصف المغناطيسية. كان هدفهم هو قياس التأخير الناتج عن هذه الرحلة، وحساب مقدار الوقت الإضافي الذي تستغرقه الإشارة للوصول مقارنة بما لو كانت قد سافرت عبر فضاء فارغ.
تكشف النتائج أن التأخير ليس رقماً ثابتاً، بل هو متغير يعتمد بشكل كبير على الزاوية التي تدخل بها الإشارة إلى الغلاف الجوي وعلى الحالة الراهنة للطقس الفضائي. فبالنسبل للإشارات التي تلامس حافة الغلاف الجوي بزوايا منخفضة، يمكن أن يكون التأخير هائلاً، حيث يمدد المسافة الظاهرية للقمر الصناعي بأكثر من مائة متر خلال فترات النشاط الشمسي المكثف. أما بالنسبة للإشارات التي تمر في مستويات أعلى، فإن التأخير يكون أصغر بكثير، وغالباً ما يكون متراً أو مترين فقط. ووجد الباحثون أن انحناء مسار الإشارة نفسه يضيف مسافة إضافية صغيرة ولكنها قابلة للقياس، منفصلة عن تأثير الإبطاء الناتج عن الغاز المشحون. كما اكتشفوا أن قوة الإشارة تلعب دوراً حاسماً في الخطأ الإجمالي؛ فالإشارات التي تمر عبر قدر أكبر من الغلاف الجوي تتأخر أكثر، لكنها تصل بإشارة أقوى وأوضح وأقل تأثراً بالضوضاء الإلكترونية في جهاز الاستقبال. وعلى العكس من ذلك، فإن الإشارات التي تمر عبر قدر أقل من الغلاف الجوي تصل بشكل أسرع ولكنها تكون أخفت وأكثر عرضة للضوضاء. وهذا يخلق مقايضة صعبة للمهندسين الذين يصممون أنظمة الملاحة القمرية: إذ يتعين عليهم الاختيار بين إشارة متأخرة ولكنها واضحة، أو إشارة سريعة ولكنها مليئة بالضوضاء.
كما أظهرت الدراسة مدى حساسية هذه التأخيرات تجاه الشمس والمجال المغناطيسي للأرض. فعندما يكون النشاط الشمسي مرتفعاً، تتمدد الطبقات المشحونة حول الأرض وتصبح أكثر كثافة، مما يسبب تأخيرات أكبر بكثير لجميع الإشارات. وعندما يضطرب المجال المغناطيسي للأرض بسبب العواصف، تتقلص الحدود الخارجية للغلاف البلازمي، وهو أمر يقلل بشكل مفاجئ من إجمالي التأخير لبعض المسارات لأن الإشارة تقضي وقتاً أقل في الجزء الأكثر كثافة من الغاز. وقد اختبر الباحثون ترددين شائعين يُستخدمان في نظام (GPS)، ووجدا أن إشارة التردد الأدنى تعاني من تأخير أكبر بكثير ولكنها تمتلك أيضاً ضوضاء أقل بكثير، بينما إشارة التردد الأعلى تعاني من تأخير أصغر ولكنها أكثر عرضة للضوض {الاضطرابات}. ومن خلال تحديد هذه التأثيرات كمياً، توفر الورقة البحثية دليلاً واضحاً لبناء أنظمة ملاحة قوية للقمر. فهي توضح أنه لتحقيق تحديد دقيق للمواقع، ستحتاج المهمات القمرية المستقبلية إلى مراعاة هذه التأخيرات الجوية باستخدام نماذج متطورة يمكنها التكيف مع الظروف المتغيرة للطقس الفضائي، لضمان سماع الهمسات الخافتة لأقمار الأرض الصناعية بوضوح عبر الفراغ.
ملخص تقني: توصيف تأخير الأيونوسفير والبلازموسفير لمستقبلات GNSS الأرضية مع نموذج البلازما الأساسي العالمي
بيان المشكلة أثبتت التطورات الأخيرة جدوى استخدام إشارات الأنظمة العالمية للملاحة عبر الأقمار الصناعية (GNSS) الأرضية، بما في ذلك عمليات الإرسال من الفصوص الجانبية الضعيفة، لأغراض تحديد الموقع والملاحة والتوقيت (PNT) عند القمر. ومع ذلك، تظل التأخيرات غير المبرمجة في الأيونوسفير والبلازموسفير مصدراً كبيراً للخطأ بالنسبة للمستخدمين القمريين. وخلافاً للمستخدمين الأرضيين الذين يواجهون مسارات إشارة شبه عمودية، يواجه المستقبلون القمريون إشارات تعبر الأيونوسفير والبلازموسفير بشكل أفقي عبر مسارات ممتدة. وبينما يمكن للمستقبلات ثنائية التردد تخفيف آثار الأيونوسفير من الدرجة الأولى، إلا أنها تؤدي إلى زيادة ضجيج المستقبل، وهو أمر بالغ الأهمية عند تتبع إشارات الفصوص الجانبية الضعيفة. علاوة على ذلك، فإن النماذج التجريبية القياسية (مثل Klobuchar وNTCM-G) لا يمكن تطبيقها مباشرة على الهندسات القمرية. ويتطلب التنبؤ الدقيق بهذه التأخيرات نمذجة صريحة لكثافة الإلكترونات في البلازموسفير، وهي منطقة تمتد من ارتفاع 1,000 كم تقريباً إلى ما بعد 20,000 كم، والتي تساهم بشكل كبير في المحتوى الكلي للإلكترونات (TEC) على طول المسارات المرتفعة أو الليلية.
المنهجية قام المؤلفون بتطوير إطار محاكاة لتوصيف تأخيرات المجموعة من الدرجة الأولى والثانية والثالثة، بالإضافة إلى تأخيرات المسار الزائدة الناتجة عن انحناء الأشعة، لإشارات GNSS التي يتم استقبالها في مدار قمري وعند القطب الجنوبي للقمر.
نمذجة البلازما: تستخدم الدراسة نموذج البلازما الأساسي العالمي (GCPM)، وهو نموذج تجريبي لكثافة الإلكترونات في الغلاف المغناطيسي الداخلي. تم اختيار GCPM لقدرته على دمج ظروف النشاط الجيومغناطيسي (مؤشر Kp) والنشاط الشمسي (مؤشر R12) المتغيرة. يدمج النموذج نماذج فرعية للأيونوسفير، والبلازموسفير، والحد الفاصل للبلازما (plasmapause)، والمنخفض المغناطيسي (magnetospheric trough)، والمنطقة القطبية.
خوارزمية تتبع الأشعة: تم تنفيذ خوارزمية مخصصة وفعالة حاسوبياً لتتبع الأشعة لحل مسارات الأشعة المنحنية بشكل تكراري. تستخدم الطريقة نهج "الإطلاق" حيث يتم تقسيم مسار الشعاع إلى قسم البلازموسفير (يتم نمذجته عبر معادلات تفاضلية) وقسم الفراغ (يُفترض أنه مستقيم لما بعد قطع 4 أنصاف أقطار أرضية). تستخدم الخوارزمية طريقة Nelder–Mead للتحسين لتحديد اتجاه الشعاع الأولي الذي يقلل المسافة بين المسار المنتشر وموقع المستقبل، مما يتجنب الحاجة إلى تقييمات التدرج أو استدعاءات كثافة الإلكترونات المتكررة أثناء التحسين.
حساب التأخير: تقوم المحاكاة بحساب معاملات الانكسار للطور والمجموعة لاستخلاص تأخيرات الدرجة الأولى والثانية والثالثة. وهي تفصل صراحةً المحتوى الكلي للإلكترونات (TEC) إلى مكونات خط البصر (LOS) ومكونات ناتجة عن الانحناء (ΔTECbend). يتم حساب تأخير المجموعة الكلي كمجموع تأخير مجموعة خط البصر (LOS)، والتأخيرات من درجات أعلى، وتأخيرات الدرجة الأولى الناتجة عن الانحناء، وطول المسار الزائد بسبب انحناء الشعاع.
تكوين المحاكاة: أُجريت عمليات المحاكاة لخمسة أقمار صناعية تابعة لنظام الاتصالات والملاحة والارتباط القمري (LCRNS) في مدار قمري، ولمستخدم عند القطب الجنوبي للقمر. حللت الدراسة إشارات GPS L1 وGalileo E1 تحت ظروف شمسية (R12 = 10 إلى 200) وجيومغناطيسية (Kp = 1 إلى 9) متغيرة. كما تم نمذجة ضجيج المستقبل بناءً على نسبة الموجة الحاملة إلى الضجيج (C/N0) وخصائص الضجيج الحراري لحلقة قفل التأخير (DLL).
المساهمات الرئيسية
تطوير الخوارزمية: تقترح الورقة خوارزمية تتبع أشعة منخفضة التكلفة باستخدام طريقة Nelder–Mead لتحديد مسارات الأشعة المنحنية بشكل تكراري دون الحاجة إلى تقييمات التدرج، مما يضمن تقارب التأخيرات المحسوبة.
توصيف شامل للتأخير: تقدم دراسة تفصيلية لتأخيرات المجموعة من الدرجة الأولى والثانية والثالثة وتأخيرات الانحناء لكل من المستخدمين في المدار القمري والمستخدمين على السطح، وهو نطاق لم يغطّه أحد من قبل لإشارات GNSS الأرضية المرصودة في المدار القمري.
التحليل البارامتري: تحدد الدراسة كمياً تأثير تردد الإشارة، والنشاط الجيومغناطيسي (مؤشر Kp)، والنشاط الشمسي (مؤشر R12) على مقادير التأخير وضجيج المستقبل.
النتائج
مقادير التأخير: عادة ما تكون متوسطات تأخيرات المجموعة في حدود 1 متر، ولكن يمكن أن تتجاوز 100 متر للمسارات ذات الارتفاعات المنخفضة (الارتفاعات المماسية < 500 كم) خلال فترات النشاط الشمسي العالي. بالنسبة للارتفاعات المماسية التي تقل عن 500 كم، يصل متوسط تأخير المجموعة من الدرجة الأولى إلى 50.59 م.
الاعتماد على الارتفاع: تتناقص التأخيرات بسرعة مع الارتفاع؛ ففوق 10,000 كم، تنخفض متوسطات تأخيرات المجموعة إلى أقل من 0.2 م. وتكون تأخيرات الانحناء أصغر عموماً من حدود الدرجة الأولى، ولكن يمكن أن تتجاوز 1 متر للمسارات ذات الارتفاعات المماسية تحت 1,000 كم. ومن الملاحظ أن التأخير الإضافي الناتج عن الانحناء غالباً ما يتجاوز الزيادة الفعلية في طول المسار الهندسي.
الآثار من درجات أعلى: تظل تأخيرات الدرجة الثانية والثالثة ضئيلة، حيث أنها أصغر بعدة مراتب من مساهمة الدرجة الأولى.
التأثير الجيومغناطيسي: خلافاً لبعض التوقعات الأرضية، تؤدي مؤشرات Kp الأعلى (التي تشير إلى نشاط جيومغناطيسي أقوى) عموماً إلى تأخيرات مجموع أصغر. ويعزى ذلك إلى تآكل البلازموسفير، حيث ينكمش الحد الفاصل للبلازما من 4-6 أنصاف أقطار أرضية (RE) في الظروف الهادئة إلى 2-3 RE أثناء النشاط القوي، مما يقلل من حجم الكثافة العالية للإلكترونات.
تأثير النشاط الشمسي: تنتج مؤشرات R12 الأعلى (النشاط الشمسي) تأخيرات مجموعة أكبر بسبب زيادة كثافة الإلكترونات في الأيونوسفير والبلازموسفير. على سبيل المثال، للحفاظ على تأخير المجموعة عند المئين التاسع والتسعين دون 5 أمتار، يجب زيادة قناع الارتفاع من ~2,000 كم (عند R12 ≤ 50) إلى ~10,000 كم (عند R12 = 200).
المفاضلات الترددية: توجد مفاضلة بين "الانحياز" و"التباين" بين إشارات L1 وL5. تظهر إشارات L1 تأخيرات مجموعة (انحياز) أصغر ولكن ضجيج مستقبل (تباين) أعلى، بينما تظهر إشارات L5 تأخيرات أكبر ولكن ضجيج حراري أقل بكثير. بالنسبة للارتفاعات المماسية حتى ~5,000 كم، تحقق L1 إجمالي خطأ مكافئ للمدى (UERE) أصغر لأن تقليل تأخير المجموعة يفوق زيادة ضجيج المستقبل. وعلى العكس من ذلك، تصبح L5 ميزة للمسارات ذات الفصوص الجانبية عالية الارتفاع حيث تكون تأخيرات المجموعة بطبيعتها صغيرة والضجيج الحراري هو المهيمن.
الأهمية تضع هذه الورقة أساساً لتصميم خوارزميات تحديد مواقع وتوقيت قوية تستخدم إشارات GNSS الأرضية في الفضاء القمري (cislunar). ومن خلال توصيف حجم وتغير آثار انتشار البلازموسفير تحت ظروف النشاط الشمسي والجيومغناطيسي المتغيرة، تساهم هذه الدراسة في توجيه عملية اختيار ترددات الإشارة، وتحديد أقنعة الارتفاع لاستبعاد الإشارات ذات التأخير العالي، وتطوير تقنيات التخفيف. وتسلط النتائج الضوء على أن النمذجة الدقيقة للتأخير ضرورية لتقليل الخطأ المكافئ للمدى (UERE) في ظل ظروف ضعف الملاحظة المتأصلة في الملاحة القمرية. وسيركز العمل المستقبلي، كما ذكر المؤلفون، على التحقق من صحة هذه المحاكاة باستخدام قياسات LuGRE القادمة ودمج تقدير التأخير مباشرة في مرشحات الملاحة.