pith. sign in

arxiv: 2601.07299 · v2 · pith:2SALJS6Inew · submitted 2026-01-12 · 📊 stat.AP

Cauchy-Gaussian Overbound for Heavy-tailed GNSS Measurement Errors

Pith reviewed 2026-05-16 15:29 UTC · model grok-4.3

classification 📊 stat.AP
keywords GNSS measurement errorsheavy-tailed distributionCauchy-Gaussian overboundintegrity monitoringvertical protection levelunimodal errorsconvolution propertynavigation safety
0
0 comments X

The pith

A Cauchy core plus Gaussian tails tightly overbounds heavy-tailed GNSS errors while preserving the property through convolution.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper develops a procedure to choose parameters for a Cauchy distribution in the core region and Gaussian distributions in the tails so that the combined distribution overbounds heavy-tailed GNSS measurement errors. The method applies to both symmetric unimodal and non-symmetric unimodal errors and includes a proof that the overbounding property survives convolution. Experiments on simulated and real data show the bounds are sharp in both central and tail regions. In the position domain the approach lowers average vertical protection level by 15 percent for symmetric cases versus a single-CDF Gaussian overbound and by 21 to 47 percent for non-symmetric cases versus Navigation Discrete ENvelope and two-step Gaussian methods.

Core claim

The authors establish that heavy-tailed GNSS measurement errors can be overbounded by a Cauchy core combined with Gaussian tails, with a parameter-selection procedure that works for both symmetric and non-symmetric unimodal cases, and they prove the overbound property is preserved under convolution. Real and simulated data confirm the resulting distribution lies above the empirical distribution in both core and tails, yielding lower vertical protection levels in the position domain.

What carries the argument

The Cauchy-Gaussian overbound, formed by pairing the Cauchy distribution for the core with Gaussian tails, with parameters fixed by the paper's selection procedure so the cumulative distribution stays above the target error distribution.

If this is right

  • The overbounding property remains valid after convolution, so the method can be applied directly to position-domain calculations.
  • For symmetric unimodal errors the average vertical protection level drops by 15 percent relative to a single-CDF Gaussian overbound.
  • For non-symmetric unimodal errors the average vertical protection level drops by 21 to 47 percent relative to Navigation Discrete ENvelope and two-step Gaussian overbounds.
  • The same distribution provides a visibly tighter envelope than prior methods at both small and large error magnitudes on real data.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Navigation systems that rely on integrity monitoring could issue fewer conservative alerts when this bound replaces looser Gaussian-only envelopes.
  • The parameter procedure could be tested on other sensor types that produce heavy-tailed noise, such as inertial measurements or camera-based ranging.
  • Urban GNSS datasets with stronger multipath effects would supply a direct check on whether the unimodal premise continues to hold.

Load-bearing premise

Heavy-tailed GNSS measurement errors are unimodal and the chosen Cauchy core plus Gaussian tails stay above the actual error distribution in real data.

What would settle it

Any real GNSS error dataset in which the empirical cumulative distribution rises above the proposed Cauchy-Gaussian curve in the core or tail regions would show the overbound has been violated.

read the original abstract

Overbounds of heavy-tailed measurement errors are essential to meet stringent navigation requirements in integrity monitoring applications. This paper proposes to leverage the bounding sharpness of the Cauchy distribution in the core and the Overbounds of heavy-tailed measurement errors are essential for meeting stringent navigation requirements in integrity-monitoring applications. This paper proposes to leverage the bounding sharpness of the Cauchy distribution in the core and the Gaussian distribution in the tails to tightly bound heavy-tailedglobal navigation satellite system measurement errors. We develop a procedure to determine the overbounding parameters for both symmetric unimodal (SU)and non-symmetric unimodal (NSU) heavy-tailed errors and prove that the over-bounding property is preserved through convolution. Experiment results on both simulated and real-world data sets reveal that our method can sharply boundheavy-tailed errors in both the core and tail regions. In the position domain, the proposed method reduces the average vertical protection level by 15% for SU heavy-tailed errors compared with the single-cumulative-density-function Gaussian overbound and by 21%-47% for NSU heavy-tailed errors compared with the navigation discrete envelope and two-step Gaussian overbounds.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper proposes a Cauchy-Gaussian overbound for heavy-tailed GNSS measurement errors in integrity monitoring. It develops parameter-fitting procedures for both symmetric unimodal (s.u.) and non-symmetric unimodal (n.s.u.) cases, proves that the overbounding property is preserved under convolution, and reports experiments on simulated and real data showing tighter bounds than single-CDF Gaussian or other methods, with average vertical protection level (VPL) reductions of 15% for s.u. errors and 21-47% for n.s.u. errors.

Significance. If the fitting procedure reliably produces a valid overbound for the range of heavy-tailed distributions encountered in GNSS, the sharper core-and-tail bounding could meaningfully improve RAIM/ARAIM performance by lowering protection levels while maintaining integrity, which is a practical advance for high-precision navigation applications.

major comments (3)
  1. [§4] §4 (parameter selection procedure): the quantile- or moment-matching steps for choosing Cauchy scale/location and Gaussian variance are data-dependent and lack a general proof or worst-case guarantee that the composite CDF strictly dominates the true error CDF for arbitrary unimodal heavy-tailed distributions beyond the calibration sets; this directly affects the validity of the claimed VPL reductions.
  2. [§3] Theorem on convolution preservation (likely §3): while the mathematical preservation under convolution is shown assuming pointwise dominance, the theorem does not address whether the fitting procedure can produce a violation for distributions with heavier tails or different core shape than those tested, rendering the preservation result conditional on an unproven fitting guarantee.
  3. [§5] Experiments (§5, real-data results): the reported 15% and 21-47% VPL reductions presuppose that the fitted overbound holds on the actual GNSS datasets, but no cross-validation or out-of-sample tail-quantile checks are described to confirm dominance at all quantiles.
minor comments (2)
  1. [Introduction] The definitions of s.u. and n.s.u. should be stated explicitly in the introduction rather than deferred to §2.
  2. [Figures] Figure captions for the CDF plots should include the specific quantile ranges used for visual comparison to allow readers to assess core vs. tail behavior.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive comments, which help clarify the scope and limitations of our Cauchy-Gaussian overbound. We address each major point below with clarifications and proposed revisions. The fitting procedure is empirical and tailored to GNSS-like unimodal heavy-tailed errors; we do not claim a universal guarantee for arbitrary distributions.

read point-by-point responses
  1. Referee: [§4] §4 (parameter selection procedure): the quantile- or moment-matching steps for choosing Cauchy scale/location and Gaussian variance are data-dependent and lack a general proof or worst-case guarantee that the composite CDF strictly dominates the true error CDF for arbitrary unimodal heavy-tailed distributions beyond the calibration sets; this directly affects the validity of the claimed VPL reductions.

    Authors: We agree that the parameter selection relies on quantile/moment matching calibrated to observed GNSS error characteristics and lacks a general worst-case proof for all possible unimodal heavy-tailed distributions. The procedure is designed to ensure pointwise dominance for the symmetric and asymmetric unimodal cases encountered in GNSS, as verified through extensive simulation and real-data checks in the manuscript. We will revise §4 to explicitly state the assumptions (unimodality and heavy tails consistent with GNSS multipath/ionospheric errors), add a limitations paragraph, and include additional synthetic tests with varied tail indices to illustrate robustness. The reported VPL reductions remain valid under the verified dominance on the datasets used. revision: partial

  2. Referee: [§3] Theorem on convolution preservation (likely §3): while the mathematical preservation under convolution is shown assuming pointwise dominance, the theorem does not address whether the fitting procedure can produce a violation for distributions with heavier tails or different core shape than those tested, rendering the preservation result conditional on an unproven fitting guarantee.

    Authors: The theorem in §3 establishes that pointwise CDF dominance is preserved under convolution, which is the key integrity property needed for RAIM/ARAIM. We will revise the theorem statement and surrounding text to clarify that preservation holds conditional on the fitted Cauchy-Gaussian being a valid overbound for the input distribution. The fitting procedure is shown empirically to produce valid overbounds for the GNSS-relevant class of distributions tested; we do not claim it succeeds for arbitrary heavier tails or core shapes outside this class. A new remark will be added noting this scope. revision: partial

  3. Referee: [§5] Experiments (§5, real-data results): the reported 15% and 21-47% VPL reductions presuppose that the fitted overbound holds on the actual GNSS datasets, but no cross-validation or out-of-sample tail-quantile checks are described to confirm dominance at all quantiles.

    Authors: The manuscript verifies dominance by direct CDF comparison on the full real and simulated datasets at core and tail quantiles. To address the concern, we will add a cross-validation subsection in §5: the real GNSS datasets will be split into calibration and test portions, with out-of-sample tail-quantile checks (e.g., 1-10^{-6} levels) reported to confirm the overbound continues to hold. This will strengthen the empirical support for the VPL reductions without altering the numerical results. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper's central derivation consists of a procedure for selecting Cauchy-Gaussian overbound parameters for unimodal heavy-tailed errors plus an explicit mathematical proof that the overbounding property is preserved under convolution. Neither step reduces to a self-definition, a fitted input renamed as prediction, or a load-bearing self-citation chain; the proof is presented as an independent mathematical argument, and the experimental VPL reductions are reported as empirical outcomes rather than forced by construction. The derivation therefore remains self-contained against external benchmarks and does not exhibit any of the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The approach depends on fitting parameters to achieve the overbound and assumes the error distributions are unimodal, with the convolution property proven from those assumptions.

free parameters (2)
  • Cauchy scale and location parameters
    Chosen via the procedure to tightly bound the core region of the heavy-tailed error distribution.
  • Gaussian variance parameter
    Selected to ensure the tails are overbounded while maintaining the overall bound.
axioms (1)
  • domain assumption Measurement errors follow a unimodal distribution (symmetric or non-symmetric).
    Invoked to develop separate procedures for s.u. and n.s.u. heavy-tailed errors.

pith-pipeline@v0.9.0 · 5462 in / 1423 out tokens · 47077 ms · 2026-05-16T15:29:00.990326+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.