Reconciling muted NH4 seasonality with observed concentrations in 2-D riparian transect model with low-K soils (silty clay loam)

[Bisesh Joshi]

Dear Pflotran experts,

I am modeling reactive transport in a fine-grained riparian sediment (silty clay loam) using PFLOTRAN.

Using Kx = 1–5 × 10⁻¹⁴ m² and Kz = 1–5 × 10⁻¹⁵ m² produces steady-state average ammonium concentrations that are reasonably consistent with observations. However, the simulated Darcy velocities are quite low (~0.02–0.1 m/yr), corresponding to porewater velocities of roughly ~0.04–0.2 m/yr assuming porosity ≈ 0.5. This results in long residence times and muted seasonal variability in transient simulations.

When I increase permeability by one to two orders of magnitude to improve connectivity and advective transport, ammonium is flushed out and no longer matches observations.

I am trying to understand the most defensible way to reconcile realistic solute concentrations with more realistic transport dynamics. Has anyone encountered a similar tradeoff?

Would you recommend:

• Using an effective/mobile porosity instead of total porosity for transport calculations?

• Representing preferential flow paths using a dual-domain or multiple-continuum approach?

• Keeping the lower permeability and interpreting the system as transport-limited?

• Another approach altogether?

Any suggestions or experiences with similar riparian reactive transport systems would be greatly appreciated.

Thank you!

Best regards,

Bisesh

--

Bisesh Joshi

PhD student

Water Science and Policy

Graduate College

University of Delaware

901-297-6761

[Hammond, Glenn E]

Bisesh,

I would try using sorption through a linear KD to mimic mobile/immobile porosity behavior.

Glenn

From: pflotra...@googlegroups.com <pflotra...@googlegroups.com> on behalf of Bisesh Joshi <bjo...@udel.edu>
Date: Monday, June 22, 2026 at 4:30 PM
To: pflotra...@googlegroups.com <pflotra...@googlegroups.com>
Subject: [pflotran-users: 8805] Reconciling muted NH4 seasonality with observed concentrations in 2-D riparian transect model with low-K soils (silty clay loam)

Check twice before you click! This email originated from outside PNNL.

--
You received this message because you are subscribed to the Google Groups "pflotran-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to pflotran-user...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/pflotran-users/CAEvi1yQnRNud97QebURZJQ17C5M6v6%2BAybgOkbBx2arSVhPZtQ%40mail.gmail.com.