Layer-Stripping Velocity Analysis Method for GPR/LPR Data

Diffraction-based velocity analysis is a key data interpretation technique in geophysical exploration, typically relying on the geometric characteristics, energy distribution, or propagation paths of diffraction waves. The hyperbola-based method is a classical strategy in this category, which extracts depth-dependent velocity (or dielectric properties) by correlating the hyperbolic shape of diffraction events with subsurface parameters for characterizing subsurface structures and material compositions. In this study, we propose a layer-stripping velocity analysis method applicable to ground-penetrating radar (GPR) and lunar-penetrating radar (LPR) data, with two main innovations: (1) replacing traditional local optimization algorithms with an intuitive parallelism check scheme, eliminating the need for complex nonlinear iterations; (2) performing depth-progressive velocity scanning of radargram diffraction signals, where shallow-layer velocity analysis constrains deeper-layer calculations. This strategy avoids misinterpretations of deep geological objects’ burial depth, morphology, and physical properties caused by a single average velocity or independent deep-layer velocity assumptions. The workflow of the proposed method is first demonstrated using a synthetic rock-fragment layered model, then applied to derive the near-surface dielectric constant distribution (down to 27 m) at the Chang’e-4 landing site. The estimated values range from 2.55 to 6, with the depth-dependent profile revealing lunar regolith stratification and interlayer material property variations. Consistent with previously reported results for the Chang’e-4 region, our findings confirm the method’s applicability to LPR data, providing a new technical framework for high-resolution subsurface structure reconstruction.