Overview
This research investigates the long-term evolution, capacity, and sustainability of the Low Earth Orbit (LEO) environment through the lens of source-sink models, specifically addressing the implications of narrow orbital shells on model accuracy. The study formulates the multi-shell, multi-species LEO environment as a Markov jump process, deriving a stochastic differential equation (SDE) approximation from this framework. This SDE approximation, which accounts for fluctuations scaling as $V^{-1/2}$ in the per-shell volume $V$, recovers the conventional source-sink ordinary differential equation (ODE) as a large-volume limit.
Research Context
Deterministic source-sink models are commonly employed to assess the long-term dynamics of the LEO environment. These models propagate shell-averaged populations using ordinary differential equations, predicated on the assumption that individual collision, disposal, and decay events average out within sufficiently large altitude shells. However, with constellation traffic increasingly organized into kilometer- and sub-kilometer-scale shells, this averaging assumption faces strain.
Approach
The researchers formulated the LEO environment, characterized by multiple shells and species, as a Markov jump process. From this formulation, they derived a stochastic differential equation (SDE) approximation. The SDE introduces fluctuations that scale inversely with the square root of the per-shell volume ($V^{-1/2}$). The conventional source-sink ODE is identified as the large-volume limit of this Markov jump process formulation. To validate the diffusion approximation provided by the SDE, it was compared against an exact discrete-event simulation of the underlying jump process.
The validation process involved sweeping the shell volume $V$ while maintaining a fixed spatial density across the 450-800 km altitude band. This approach allowed for a direct comparison of the SDE and ODE descriptions under varying shell definitions.
Findings
- The conventional source-sink ODE is recovered as a large-volume limit of the multi-shell, multi-species LEO environment formulated as a Markov jump process.
- A stochastic differential equation (SDE) approximation was derived from this formulation, indicating fluctuations that scale as $V^{-1/2}$ in the per-shell volume $V$.
- Validation against an exact discrete-event simulation confirmed the SDE's diffusion approximation.
- At shell volumes comparable to those utilized in established source-sink models, the SDE and ODE descriptions exhibited agreement.
- As shells narrowed, the two descriptions diverged.
- At the finest shell volume considered, the mean debris population projected by the stochastic model reached approximately 4.5 times the prediction from the ODE.
- Several realizations of the stochastic model exhibited runaway growth, a phenomenon absent from the ODE trajectory.
- This departure is attributed to nonlinear collision terms, which lead to increased expected collision rates through population variance and covariance. This, in turn, generates further debris and reinforces the collision-debris feedback loop.
- Narrow-shell stochasticity both widens the distribution of outcomes and shifts the expected trajectory of the LEO environment evolution.
Why This Matters
The stochastic model shares its parameterization with the deterministic source-sink models, offering a scale-consistent extension. This allows for evaluation of shell configuration, collision risk, and long-term LEO sustainability, particularly as orbital traffic becomes increasingly concentrated in narrow shells. The identified divergence in predictions highlights the importance of incorporating stochastic effects for accurate assessment of the LEO environment's future state.