Model Idea

Lagging Theory revisits one groundwater-flow assumption. Classical Darcy’s law relates water flux to the hydraulic gradient at the same instant at the continuum scale. Flux means water moving through the aquifer; the gradient describes the spatial change in head or drawdown that drives that flow. Lin and Yeh (2017) allowed flux and drawdown gradient to adjust on different macroscopic response times in a constant-rate pumping model.

In this context, “lag” can change the shape of a response, rather than simply shift a curve along the time axis. The two lag parameters summarize a delayed relation at the model’s observation scale. They do not identify a unique microscopic cause.

The two lag times in the 2017 formulation have distinct roles, using the same notation as the demonstration below:

  • flux-response lag, τq, belongs to the water-flux side of the relation;
  • drawdown-gradient lag, τs, belongs to the drawdown-gradient side.

For the initial conditions used in this demo, equal lags cancel in the transfer relation; setting both to zero gives the classical Darcy limit directly. Try the Classical Darcy preset in the pumping demo, then separate the lags to see how the early response changes. This is a teaching calculation, not a calibrated field test; it omits wellbore storage and finite well radius.

Published Cases

  • Constant-rate pumping, 2017. Lin, Y.-C., and Yeh, H.-D., Water Resources Research formulated lagging flow in a leaky confined aquifer and compared its solution with a fractured-aquifer pumping record from South Dakota. The reported fit was particularly good at early pumping time, and fitted lag behavior varied with observation distance. The NYCU publication record provides the abstract and published citation.
  • Periodic hydraulic head, 2026. Lin and Kurylyk, Journal of Hydrology examined why diffusivity inferred from amplitude damping can disagree with diffusivity inferred from phase lag. Their two-lag formulation reconciled these signatures in reanalyses of existing Tuolumne Meadows and Meghna River records. See the publisher record for the study and its assumptions.

These comparisons support the model’s ability to describe the reported responses under the stated assumptions. Improved calibration fit alone does not establish prediction at a new well, frequency, or site, nor prove which physical mechanism caused a fitted lag. Flow-path adjustment, domain exchange, drainage, leakage, storage, and unresolved heterogeneity remain possible interpretations to distinguish with additional evidence. The minimum tests below describe validation goals, rather than results established for every application.

Concise Definition

Lagging Theory tests whether asynchronous hydraulic response changes interpretation. Boundary movement and thermal response are explored as extensions; each requires its own model assumptions and validation.

The testable claim is narrow. If asynchronous response matters, a lagging formulation should improve residual structure, parameter transfer, held-out prediction, and at least one engineering decision variable after complexity and identifiability checks.

Extensions

Subsequent work applies the same test to different response settings:

  • in unconfined aquifers, lag times enter the free-surface condition and represent capillary-fringe release and capillary-suction drainage;
  • in periodic head signals, flux lag and head lag separate amplitude damping from phase offset, addressing phase-amplitude diffusivity mismatch;
  • in engineering interpretation, the lagging equation becomes one candidate analytical model for transforming measured response into inferred properties and decision variables.

Lagging Theory complements Darcy, Theis, Neuman, delayed-yield, leakage, and dual-porosity models by testing one specific possibility: the hydraulic response may contain flux-gradient asynchrony.

Minimum Tests

A lagging model should pass more than calibration fit:

  1. Residual structure improves in a meaningful way.
  2. Complexity penalties do not erase the gain.
  3. Parameters are identifiable enough for the intended decision.
  4. Synthetic known-truth coverage is acceptable.
  5. Field prediction improves on held-out time, recovery, or wells.
  6. The difference propagates to an engineering decision variable.

For pumping-test applications, see When Does a Pumping Test Need Lagging Darcy Law?.

Relation to Other Non-Equilibrium Models

Mathematical overlap with dual-porosity, delayed-yield, and other non-equilibrium models is informative. It provides a compact way to test whether flow-path adjustment, inter-domain exchange, capillary drainage, hydro-mechanical coupling, or field-scale delayed response affects the interpretation.

Discuss a Research Question

For an academic discussion, use the research collaboration route. A useful starting point is a specific question about transfer across frequencies or sites, comparison with a dual-domain model, or reanalysis of pumping and recovery records. Describe the forcing, the observed mismatch, the competing explanations, and what would count as a successful test. Contact Ying-Fan Lin at yflin1110@cycu.edu.tw.