Benchmarking Cylindrical Blast Wave Theory Against the OSIRIS-REx Sample Return Capsule Reentry

This study benchmarks cylindrical blast wave theory against the OSIRIS-REx Sample Return Capsule reentry using 39 infrasound stations, identifying the Sakurai formulation as the most accurate model for predicting signal characteristics of non-ablating hypersonic bodies while demonstrating that the signal period is a robust observable for constraining blast radius.

The Gist

Imagine a giant, invisible drum being struck by a speeding bullet as it slices through the sky. When an object moves faster than sound, it creates a shockwave—a sonic boom. As this boom travels far away, it changes shape and gets quieter. Scientists have a set of mathematical "recipes" (formulas) to predict exactly how loud that boom will be and how long the "thump" lasts when it finally reaches a listener on the ground.

For decades, these recipes were tested against meteors (falling space rocks). But meteors are messy: they burn up, break apart, and change size as they fall, making it hard to know if the recipe is wrong or if the rock just behaved unexpectedly.

This paper is like a "final exam" for those recipes, but instead of using a messy meteor, the scientists used a known, perfect object: the OSIRIS-REx Sample Return Capsule. This was a spacecraft that returned to Earth in 2023. Because it was a human-made machine, scientists knew its exact size, weight, speed, and path. It didn't burn up or break apart significantly. It was a "clean" test subject.

Here is what the study found, explained simply:

1. The "Ground Truth" Experiment

Think of the 39 microphones (infrasound stations) scattered across the desert like a giant net catching the sound of the capsule's reentry. Because the capsule's path was perfectly known, the scientists could calculate exactly what the sound should have been at each microphone. They then compared the "should have been" math against the "actually heard" data.

2. The Six Recipes vs. The Three Rules

The scientists tested six different mathematical recipes for calculating the "blast radius" (how big the initial shockwave is). They also tested three different "transition rules" (mathematical switches that decide when the shockwave stops acting like a violent explosion and starts acting like a normal sound wave).

• The Winner: One specific recipe, called the Sakurai formulation, was the clear champion. It predicted the "thump" duration (the signal period) with incredible accuracy—within about 9% of what was actually heard.
• The Runner-Up: Another recipe (Jones/Plooster) was almost as good, provided the scientists used the right "transition rule."
• The Losers: Three other recipes, which are commonly used for meteors, failed miserably. They predicted the sound would last much longer than it actually did.
  • The Analogy: Imagine trying to predict how far a rubber band snaps. The "meteor" recipes assume the rubber band is sticky and leaves a trail of goo that makes it snap further. But the capsule was a clean, rigid metal ball. Using the "sticky" recipes for the "clean" ball made the prediction way too big (overestimating the blast radius by more than 3 times).

3. The "Thump" vs. The "Volume"

The study made a crucial discovery about what to measure:

• The "Thump" (Period): This is how long the sound wave lasts. The paper found that measuring the duration of the sound is a very reliable way to figure out the source's energy. It's like judging the size of a drum by how long the vibration lasts; it's stable and hard to mess up.
• The "Volume" (Amplitude): This is how loud the sound is. The study found that predicting the loudness was a disaster. No recipe could get the volume right.
  • The Analogy: Imagine trying to guess how hard a drum was hit by listening to it in a windy, echoey canyon. The length of the sound might still be clear, but the volume gets messed up by the wind, the rocks, and the echo. The paper concludes that for these types of events, you should trust the "thump" (duration) and ignore the "volume" (loudness) because the volume is too easily distorted by the atmosphere.

4. The Altitude Problem

The study also found a pattern based on height.

• When the capsule was low (thick air), the recipes slightly underestimated the sound.
• When the capsule was high (thin air), the recipes slightly overestimated the sound.
• The Analogy: It's like a map that is slightly too small for the bottom of a mountain and slightly too big for the top. The map works okay in the middle, but it drifts as you go up or down. The scientists found that the "Sakurai" recipe works best between 46 and 58 km altitude, but it starts to drift outside that range.

5. Why This Matters (According to the Paper)

The paper doesn't claim this will change how we build spaceships or treat diseases. Instead, it establishes a baseline of truth.

• It proves that for rigid, non-burning objects (like spacecraft returning to Earth), we can now use the "Sakurai" recipe to accurately estimate the energy of the event just by listening to the sound duration.
• It confirms that we should stop using the "meteor" recipes for these clean spacecraft, as they give wildly incorrect results.
• It tells future scientists: "If you want to know what happened during a reentry, measure the time the sound lasts, not how loud it is, and use the Sakurai math."

In short, the paper took a messy, complicated problem (predicting space sounds) and used a perfect, known object to figure out which math tools actually work and which ones are broken. The result is a much clearer, more accurate way to listen to the sky.

Technical Summary

Technical Summary: Benchmarking Cylindrical Blast Wave Theory Against the OSIRIS-REx Sample Return Capsule Reentry

Problem Statement
Weak shock theory, based on cylindrical blast wave models, has long been used to interpret meteor infrasound. However, previous applications have relied on ablating meteoroids where source parameters (mass, density, ablation rate) must be inferred indirectly, introducing significant uncertainty. Consequently, the operational validity of existing blast radius ($R_0$) formulations and weak-shock transition coefficients ($C$) had not been systematically benchmarked against a non-ablating hypersonic source with independently known parameters. This study addresses the need to isolate source physics from propagation effects by evaluating these formulations against a "ground truth" dataset where the source geometry, trajectory, and emission points are precisely constrained.

Methodology
The study utilizes infrasound observations from the OSIRIS-REx Sample Return Capsule (SRC) reentry on September 24, 2023. The SRC provided a unique benchmark as a rigid, effectively non-ablating body with known dimensions ($d_m = 0.81$ m), mass, and trajectory. Data were collected from 39 infrasound stations across Nevada and Utah, each sampling a distinct emission point along the trajectory (altitudes ~44–62 km).

The methodology involved a systematic forward and inverse modeling approach:

1. Forward Modeling: The study evaluated six published $R_0$ formulations (Few, Jones/Plooster, Sakurai, Tsikulin standard, Tsikulin modified, and the Mach-diameter approximation) and three weak-shock transition (WST) coefficients ($C = 5.38, 34.3, 57.2$). These were applied within a stratified atmospheric propagation model to predict signal period ($\tau$) and peak overpressure ($\Delta p$) at the receiver stations.
2. Inverse Modeling: Observed periods and amplitudes were used to invert for $R_0$ under each formulation and $C$ assumption.
3. Statistical Analysis: Model performance was quantified using signed percent residuals and median absolute percentage residuals (MAPR) across the 39 stations. Statistical significance was assessed using bootstrap confidence intervals, paired Wilcoxon signed-rank tests, and Friedman repeated-measures tests.

Key Contributions

• First Controlled Benchmark: This work provides the first systematic benchmarking of cylindrical blast wave formulations against a non-ablating hypersonic source with independently known source conditions, removing the ambiguity inherent in meteor studies.
• Formulation Ranking: It establishes a performance hierarchy for $R_0$ formulations specifically for rigid bodies, distinguishing between energy-normalized and scaling-type formulations.
• Observable Reliability: The study rigorously demonstrates that signal period is a robust observable for constraining $R_0$, whereas peak overpressure amplitude is not, due to its high sensitivity to propagation effects and transition assumptions.
• Source-Condition Dependence: It identifies a systematic, monotonic trend in residuals relative to source altitude, revealing structural limitations in constant-normalization formulations when applied to decelerating hypersonic sources.

Results

• Best-Performing Formulation: The Sakurai (1965) formulation is identified as the best-performing model for non-ablating bodies. When paired with the Towne ($C=34.3$) transition coefficient, it yields a forward-model period MAPR of 9% with a near-zero median signed residual (+1%). The Jones/Plooster formulation performs comparably (MAPR ~11%) when a physically appropriate $C$ is adopted.
• Inadequacy of Other Formulations: The three scaling-type formulations (Tsikulin standard/modified and Mach-diameter) systematically overestimate $R_0$ by factors of 2.5 to 3.5, leading to period MAPRs exceeding 80%. Specifically, the Mach-diameter approximation ($R_0 \approx M d_m$), commonly used for meteors, overestimates $R_0$ by more than a factor of 3 for the non-ablating SRC.
• Period vs. Amplitude: The signal period is only weakly sensitive to the WST coefficient $C$ at the propagation distances sampled (~7,000–10,000 $R_0$). In contrast, amplitude predictions are highly sensitive to $C$ and show large discrepancies (MAPRs of 28–578%). No single formulation successfully reconciles both period and amplitude constraints simultaneously.
• Along-Trajectory Bias: All energy-normalized formulations exhibit a strong monotonic dependence on source altitude (Spearman $r_s = 0.85$). They tend to underpredict $R_0$ at low altitudes (high density, low Mach) and overpredict at high altitudes (low density, high Mach). The Sakurai formulation achieves residuals within ±25% for source altitudes between approximately 46 and 58 km.

Significance and Claims
The paper claims to establish a quantitative performance baseline for applying weak-shock theory to non-ablating hypersonic reentries. It demonstrates that for rigid bodies, the signal period is the operationally reliable observable for constraining source-level blast radii and energy deposition per unit path length, while amplitude remains poorly constrained by current semi-analytical models.

The study concludes that the Sakurai formulation, potentially with altitude-dependent corrections to address the observed bias, provides the most accurate framework for interpreting infrasound from sample return capsules and similar non-ablating sources. The authors note that while the Mach-diameter approximation is unsuitable for rigid bodies, the period-based inversion offers a practical pathway for source characterization in scenarios where only infrasound data are available, such as reentries over remote regions. The work does not propose new theoretical frameworks but rather validates and refines the operational application of existing cylindrical blast wave theory.