Propagation of Non-Gaussian Uncertainties from Storm-Time Neutral Density to Orbital Position

Propagation of Non-Gaussian Uncertainties from Storm-Time Neutral Density to Orbital Position | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Short Report Propagation of Non-Gaussian Uncertainties from Storm-Time Neutral Density to Orbital Position Emmett Peters, Nicholas Dietrich, Tomoko Matsuo This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-8042026/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Satellites in low Earth orbit (LEO) are subject to atmospheric drag driven by neutral density, which is highly variable and amplified during geomagnetic storms. Accurate drag modeling during storms is critical for collision avoidance and space situational awareness (SSA). While orbit prediction errors are often modeled under Gaussian assumptions and represented with covariance ellipsoids, prior work suggests that density distributions are non-Gaussian especially during geomagnetic storms. The impact of this variability on orbital position uncertainty has not been quantified.In this brief report, we propagate circular LEO orbits through Monte Carlo density fields, obtained from a particle filter neutral density estimation framework that employs the Thermosphere Ionosphere Electrodynamics General Circulation Model (TIEGCM), during both quiet and storm periods. Reference orbits without drag are compared against drag-included trajectories, and position differences are analyzed in the in-track direction. The statistics of density and in-track position differences are presented.The results show that storm-time density exhibits larger mean and variance compared to quiet-time conditions, as well as elevated skewness and kurtosis, confirming its non-Gaussian structure. These properties drive in-track position differences to be larger with a higher variance and to deviate from Gaussian distributions. The magnitude of storm-time in-track position differences resulting from these density uncertainties is comparable to position uncertainties caused by initial orbital state parameter errors typically considered in SSA orbital determination applications.This study suggests the importance of accounting for realistic neutral density variability in drag modeling to improve the reliability of orbit predictions during geomagnetic storms. neutral mass density atmospheric drag orbital errors geomagnetic storms non-Gaussian position errors Full Text Additional Declarations No competing interests reported. Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-8042026","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Short Report","associatedPublications":[],"authors":[{"id":540988710,"identity":"a7dd53ba-008a-47c9-8d7a-0c58b0db4895","order_by":0,"name":"Emmett Peters","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABBklEQVRIiWNgGAWjYBAC+wYIXQ+hKhgY+EA0D4g4gF2LAVgugSGxgQ3EOsPAwEaaFsY2YrSInU58/PEHQ2K/fPOxjz/n2cmxsfc+/PCmgkGO70YCdr9I5242ANpSL9nGljybd1uyMRvPcWPJOWcYjCVxaDGQzt0mAdSSYHCMx5iZcduBxDaJNAZp3jaGxA24tWz/AdJif4z/M+PPOWAtzL+BWurxaNnGALaFjYeZgbcBrIUNZEuCAW4tmyXOpEkkSBxLM2bmOQbyyzE2yzlnJAxnnnmA3fuzczd+qLCxSeBvPvyY8UeNnRw/exvzjTcVNvJ8x7HbAgUSRIiMglEwCkbBKCAeAAAuP1sEdwMOXQAAAABJRU5ErkJggg==","orcid":"","institution":"University of Colorado Boulder","correspondingAuthor":true,"prefix":"","firstName":"Emmett","middleName":"","lastName":"Peters","suffix":""},{"id":540988711,"identity":"e1ae96ad-e2c7-40fd-b4db-12c7010f089e","order_by":1,"name":"Nicholas Dietrich","email":"","orcid":"","institution":"United States Naval Research Laboratory","correspondingAuthor":false,"prefix":"","firstName":"Nicholas","middleName":"","lastName":"Dietrich","suffix":""},{"id":540988712,"identity":"ae9b6eed-2d66-4e2d-b914-ced25f2b3e07","order_by":2,"name":"Tomoko Matsuo","email":"","orcid":"","institution":"University of Colorado Boulder","correspondingAuthor":false,"prefix":"","firstName":"Tomoko","middleName":"","lastName":"Matsuo","suffix":""}],"badges":[],"createdAt":"2025-11-05 22:38:13","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-8042026/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-8042026/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":95527681,"identity":"48b85122-5fb8-4f08-923d-1aaa76594598","added_by":"auto","created_at":"2025-11-10 10:14:29","extension":"json","order_by":0,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":5547,"visible":true,"origin":"","legend":"","description":"","filename":"ec9ed2b3fa4743e78976cf9998adc507.json","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/583e3f7ac793033d407d3175.json"},{"id":95493086,"identity":"148070ad-84b9-42a1-9956-01ceadea9011","added_by":"auto","created_at":"2025-11-10 03:17:45","extension":"xml","order_by":1,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":54183,"visible":true,"origin":"","legend":"","description":"","filename":"ec9ed2b3fa4743e78976cf9998adc5071enriched.xml","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/0eff99e736656f4f887d61c3.xml"},{"id":95493093,"identity":"21c6d5b1-9273-4f0e-aa53-99d73b9611ff","added_by":"auto","created_at":"2025-11-10 03:17:45","extension":"pdf","order_by":2,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":86009,"visible":true,"origin":"","legend":"","description":"","filename":"EPetersCoverLetter.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/1eac3e9e2f2f7cdf6de4e282.pdf"},{"id":95527878,"identity":"fb770a37-a565-493c-906d-eba35c19d25d","added_by":"auto","created_at":"2025-11-10 10:15:07","extension":"pdf","order_by":3,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":932513,"visible":true,"origin":"","legend":"","description":"","filename":"EPetersPropagationNGStormDensityOrbitalPosition.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/f338d43ee349823890b0b1ab.pdf"},{"id":95493087,"identity":"6bbbce50-8bbb-4ed8-95ce-2127074fbba3","added_by":"auto","created_at":"2025-11-10 03:17:45","extension":"png","order_by":4,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":46165,"visible":true,"origin":"","legend":"","description":"","filename":"Kp.png","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/689b413c12ab461c59f59182.png"},{"id":95528122,"identity":"6c887aca-32d8-4451-aac1-23a55da863e9","added_by":"auto","created_at":"2025-11-10 10:15:38","extension":"png","order_by":5,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":276984,"visible":true,"origin":"","legend":"","description":"","filename":"density.png","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/f6f9c4f7abe514264cd8fbb7.png"},{"id":95493090,"identity":"7041e81f-5282-459a-b5b5-96a8f6c0ae9b","added_by":"auto","created_at":"2025-11-10 03:17:45","extension":"png","order_by":6,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":154640,"visible":true,"origin":"","legend":"","description":"","filename":"densitystats.png","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/0b33d1da951f6e12c4e63ca6.png"},{"id":95493088,"identity":"aca86a9a-e3ce-42c5-a868-b4b3e4cea4c9","added_by":"auto","created_at":"2025-11-10 03:17:45","extension":"eps","order_by":7,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":2890,"visible":true,"origin":"","legend":"","description":"","filename":"empty.eps","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/af15105c1b2ab4c4e6c09b4d.eps"},{"id":95528141,"identity":"7c1a5c06-c2e8-41b6-b2c2-b307925af29e","added_by":"auto","created_at":"2025-11-10 10:15:39","extension":"eps","order_by":8,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":91593,"visible":true,"origin":"","legend":"","description":"","filename":"fig.eps","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/0b2a29d04dc3da8622a69211.eps"},{"id":95528535,"identity":"39d4814c-e403-4daa-ad7d-5fa884e4ce24","added_by":"auto","created_at":"2025-11-10 10:16:15","extension":"png","order_by":9,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":219422,"visible":true,"origin":"","legend":"","description":"","filename":"intrack.png","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/46c700850576e7e6c2b533c7.png"},{"id":95493095,"identity":"ca6fc349-b091-49f6-a25b-a3f0baca742d","added_by":"auto","created_at":"2025-11-10 03:17:45","extension":"png","order_by":10,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":171546,"visible":true,"origin":"","legend":"","description":"","filename":"intrack1.png","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/9a56a9d45041794d5b29e311.png"},{"id":95528467,"identity":"6572cb52-eb33-42bd-923a-9a57f63a6690","added_by":"auto","created_at":"2025-11-10 10:16:10","extension":"png","order_by":11,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":100794,"visible":true,"origin":"","legend":"","description":"","filename":"intrackstats.png","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/c8a9c5c79b983be8853e0b90.png"},{"id":95527754,"identity":"9864bdf7-b81e-4dd4-9aba-36f5d01651bd","added_by":"auto","created_at":"2025-11-10 10:14:41","extension":"bst","order_by":12,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":146013,"visible":true,"origin":"","legend":"","description":"","filename":"snapacite.bst","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/f46aa23f6d9901c7439633b6.bst"},{"id":95493106,"identity":"085509a8-443d-41e6-9772-b85058707105","added_by":"auto","created_at":"2025-11-10 03:17:45","extension":"bst","order_by":13,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":29828,"visible":true,"origin":"","legend":"","description":"","filename":"snaps.bst","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/b80d64792472c4e7d24cecb3.bst"},{"id":95528184,"identity":"de0cac3c-d699-44b5-9eb5-f6f7ebb04c41","added_by":"auto","created_at":"2025-11-10 10:15:41","extension":"bst","order_by":14,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":35515,"visible":true,"origin":"","legend":"","description":"","filename":"snbasic.bst","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/1d5dcb965d1d6f3ee888609a.bst"},{"id":95493107,"identity":"389ff70b-79a7-4796-80eb-64e766bd0a92","added_by":"auto","created_at":"2025-11-10 03:17:46","extension":"bst","order_by":15,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":33968,"visible":true,"origin":"","legend":"","description":"","filename":"snchicago.bst","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/2e10bc3eb88303b00e007697.bst"},{"id":95527525,"identity":"bcb17da2-f1da-4e97-baae-d1e34af0c5e0","added_by":"auto","created_at":"2025-11-10 10:14:01","extension":"cls","order_by":16,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":55857,"visible":true,"origin":"","legend":"","description":"","filename":"snjnl.cls","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/271a6121d0adc9e3ad7ad3e9.cls"},{"id":95493097,"identity":"1647762f-bef0-4fb8-b06c-e640a0575728","added_by":"auto","created_at":"2025-11-10 03:17:45","extension":"bst","order_by":17,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":64023,"visible":true,"origin":"","legend":"","description":"","filename":"snmathphysay.bst","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/f623f7bcd745b4d7ffda9fc2.bst"},{"id":95493110,"identity":"29646195-aa67-4f6d-bfaf-d137e07a3dfc","added_by":"auto","created_at":"2025-11-10 03:17:46","extension":"bst","order_by":18,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":64166,"visible":true,"origin":"","legend":"","description":"","filename":"snmathphysnum.bst","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/6053df381d78d14777332faa.bst"},{"id":95528481,"identity":"f2fff774-c34f-4cd3-9b84-99c913449ef3","added_by":"auto","created_at":"2025-11-10 10:16:11","extension":"bst","order_by":19,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":37333,"visible":true,"origin":"","legend":"","description":"","filename":"snnature.bst","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/929022c4dd8e5469703c6c1b.bst"},{"id":95493104,"identity":"a13b2beb-b67d-4063-aea1-932bc8aff551","added_by":"auto","created_at":"2025-11-10 03:17:45","extension":"bst","order_by":20,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":39951,"visible":true,"origin":"","legend":"","description":"","filename":"snvancouveray.bst","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/efae6b97208079f4300e58c8.bst"},{"id":95493099,"identity":"2144cf75-f648-476f-bf0c-7e400074146b","added_by":"auto","created_at":"2025-11-10 03:17:45","extension":"bst","order_by":21,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":40758,"visible":true,"origin":"","legend":"","description":"","filename":"snvancouvernum.bst","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/4f891a37d2bfce994ad91dc7.bst"},{"id":95493111,"identity":"46901276-e2ce-4509-b160-7df19681235c","added_by":"auto","created_at":"2025-11-10 03:17:46","extension":"png","order_by":26,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":146696,"visible":true,"origin":"","legend":"","description":"","filename":"Onlineintrack1.png","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/3d84460da8d43a32db957c04.png"},{"id":95493103,"identity":"5ce89051-0cd9-4299-a7c5-a70ebef149cb","added_by":"auto","created_at":"2025-11-10 03:17:45","extension":"xml","order_by":28,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":59010,"visible":true,"origin":"","legend":"","description":"","filename":"ec9ed2b3fa4743e78976cf9998adc5071structuring.xml","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/343621d20e125a4edbb770ce.xml"},{"id":95493108,"identity":"4ea318b8-f32b-441d-b7e7-c77f755df979","added_by":"auto","created_at":"2025-11-10 03:17:46","extension":"html","order_by":29,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":66193,"visible":true,"origin":"","legend":"","description":"","filename":"earlyproof.html","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1/566ea93f89ab85fc30da0007.html"},{"id":99641356,"identity":"8ece0273-987b-4554-b3f1-8dfc9222d55b","added_by":"auto","created_at":"2026-01-06 18:39:13","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":870897,"visible":true,"origin":"","legend":"","description":"","filename":"EPetersPropagationNGStormDensityOrbitalPosition.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8042026/v1_covered_5307b994-6dd5-477b-b9c5-e697352d8544.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"Propagation of Non-Gaussian Uncertainties from Storm-Time Neutral Density to Orbital Position","fulltext":[],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":false,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":true,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":true,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"neutral mass density, atmospheric drag, orbital errors, geomagnetic storms, non-Gaussian position errors","lastPublishedDoi":"10.21203/rs.3.rs-8042026/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-8042026/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eSatellites in low Earth orbit (LEO) are subject to atmospheric drag driven by neutral density, which is highly variable and amplified during geomagnetic storms. Accurate drag modeling during storms is critical for collision avoidance and space situational awareness (SSA). While orbit prediction errors are often modeled under Gaussian assumptions and represented with covariance ellipsoids, prior work suggests that density distributions are non-Gaussian especially during geomagnetic storms. The impact of this variability on orbital position uncertainty has not been quantified.In this brief report, we propagate circular LEO orbits through Monte Carlo density fields, obtained from a particle filter neutral density estimation framework that employs the Thermosphere Ionosphere Electrodynamics General Circulation Model (TIEGCM), during both quiet and storm periods. Reference orbits without drag are compared against drag-included trajectories, and position differences are analyzed in the in-track direction. The statistics of density and in-track position differences are presented.The results show that storm-time density exhibits larger mean and variance compared to quiet-time conditions, as well as elevated skewness and kurtosis, confirming its non-Gaussian structure. These properties drive in-track position differences to be larger with a higher variance and to deviate from Gaussian distributions. The magnitude of storm-time in-track position differences resulting from these density uncertainties is comparable to position uncertainties caused by initial orbital state parameter errors typically considered in SSA orbital determination applications.This study suggests the importance of accounting for realistic neutral density variability in drag modeling to improve the reliability of orbit predictions during geomagnetic storms.\u003c/p\u003e","manuscriptTitle":"Propagation of Non-Gaussian Uncertainties from Storm-Time Neutral Density to Orbital Position","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-11-10 03:17:40","doi":"10.21203/rs.3.rs-8042026/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"6706ce57-424f-4a10-a7d1-f92262323f77","owner":[],"postedDate":"November 10th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[],"tags":[],"updatedAt":"2026-01-06T18:39:02+00:00","versionOfRecord":[],"versionCreatedAt":"2025-11-10 03:17:40","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-8042026","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-8042026","identity":"rs-8042026","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}