Understanding celerity, velocity, and their connections is crucial for unraveling subsurface hydrology and improving numerical models. This observational study investigated leakage in soil columns subjected to diverse treatments under natural rainfall, aiming to characterize celerity, velocity, and their relationships during leakage, and assess treatment effects on this relationship. Over three years with 13 rainfall events, velocities (0.00-0.26 m/h) and celerities (0.09-442.45 m/h) displayed significant disparities. A fixed kinematic ratio (k = 200) partitioned the celerity-velocity relationship into High-Ratio Area and Low-Ratio Area. Assuming liquid-phase changes affect mechanical wave celerity in a three-phase medium, these areas are hypothesized to correspond to dominant leakage-stage flow patterns: film flow and macropore flow. Both the High-Ratio Area and the Low-Ratio Area showed a strong adherence to a power function model, indicating a consistent non-linear celerity-velocity relationship across leakage scenarios. Treatments applied to soil columns influenced the distribution of celerity-velocity pairs and their interaction, demonstrating the potential to modify subsurface hydrological dynamics through targeted interventions.