Posts by Nicolae Suciu


Reactive transport in saturated/unsaturated porous media is numerically upscaled to the spacetime scale of a hypothetical measurement through coarse grained space-time (CGST) averages. The one-dimensional reactive transport is modeled at the fine-grained Darcy scale by the actual number of molecules involved in reactions which undergo advective and diffusive movements described by global random walk (GRW) simulations. The CGST averages verify identities similar to a local scale balance equation which allow us to derive expressions for the flow velocity and the intrinsic diffusion coefficient in terms of averaged microscopic quantities. The latter are further used to verify the CGST-GRW numerical approach. The upscaling approach is applied to biodegradation processes in saturated aquifers and variably saturated soils and the CGST averages are compared to classical volume averages. One finds that if the process is characterized by slow variations in time, as in homogeneous systems or in case of observations of reactive transport in heterogeneous aquifers made at large times or far away from the contaminant source, the differences between the two averages are negligible. Instead, the differences are significant if the averages are computed close to thesource at early times, in case of aquifer simulations, and can be extremely large in simulations of biodegradation in soils. In the latter case, the volume average is totally inappropriate as model for experimental measurements, leading for instance to overestimations by 100% of the CGST average.


Nicolae Suciu
Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy,  Cluj-Napoca, Romania

Florin A. Radu
Department of Mathematics, University of Bergen, Allegaten 41, 5007 Bergen, Norway


Space-time upscaling, Global random walk, Reactive transport, Richards equation MSC: 76S05, 35K57, 86A05, 65C35


Cite this paper as:

N. Suciu and F.A. Radu, Space-time upscaling of reactive transport in porous media, Arxiv:

About this paper


Arxiv (preprint)

Publisher Name

Not available yet.

Print ISSN

Not available yet.

Online ISSN

Not available yet.

Google Scholar Profile

[1] Alecsa, C.D., Boros, I., Frank, F., Knabner, P., Nechita, M., Prechtel, A., Rupp, A., Suciu, N., 2019. Numerical benchmark study for fow in heterogeneous aquifers. Adv. Water Resour., 138, 103558.
[2] Andricevic R., 1998. Effects of local dispersion and sampling volume on the evolution of concentration fluctuations in aquifers. Water Resour. Res., 34(5):1115–29.
[3] de Barros, F.P.J., Fiori, A., 2014. First-order based cumulative distribution function for solute concentration in heterogeneous aquifers: Theoretical analysis and implications for human health   risk assessment. Water Resour. Res., 50(5), 4018–4037.
[4] Battiato, I. Multiscale dynamics of reactive fronts in heterogeneous media with fluctuating forcings. SIAM Conference on Mathematical and Computational Issues in the Geosciences September 11–14, 2017, Erlangen, Germany.
[5] Bayer-Raich, M., Jarsjo, J., Liedl, R., Ptak, T., Teutsch, G.2004. Average contaminant concentration and mass flow in aquifers from time-dependent pumping well data: Analytical framework. Water Resour. Res., 40(8), W08303.
[6] Bause, M., Knabner, P., 2004. Numerical simulation of contaminant biodegradation by higher order methods and adaptive time stepping. Comput. Visual. Sci., 7(2), 61–78.
[7] Berkowitz, B., 2021. HESS Opinions: Chemical transport modeling in subsurface hydrological systems – Space, time, and the holy grail of “upscaling”. Hydrol. Earth. Syst. Sci. Discuss., [preprint].
[8] Brunner, F., Radu, F.A., Bause, M., Knabner, P., 2012. Optimal order convergence of a modified BDM1 mixed finite element scheme for reactive transport in porous media. Adv. Water Resour., 35, 163–171.
[9] Bringedal, C., Berre, I., Pop, I.S., Radu, F.A., 2015. A model for non-isothermal flow and mineral precipitation and dissolution in a thin strip. J. Comput. Appl. Math., 289, 346–355.
[10] Bringedal, C., Berre, I., Pop, I. S., Radu, F. A., 2016. Upscaling of nonisothermal reactive porous media flow under dominant Peclet number: the effect of changing porosity. Multiscale Model. Simul. 14(1), 502–533.
[11] Caroni, E., Fiorotto, V., 2005. Analysis of concentration as sampled in natural aquifers. Transport Porous Media, 59(1),19–45.
[12] Cirpka, O.A., Frind, E.O., Helmig, R., 1999. Numerical simulation of biodegradation controlled by transverse mixing. J. Contam. Hydrol., 40(2), 159-182.
[13] Cushman, J.H., 1983. Multiphase transport equations: I-general equation for macroscopic statistical, local space-time homogeneity. Transport Theor. Stat. Phys., 12(1), 35–71.
[14] Destouni, G., Graham, W., 1997. The influence of observation method on local concentration statistics in the subsurface. Water Resour. Res., 33(4), 663–676.
[15] Eberhard, J.P., Suciu, N., Vamos, C., 2007. On the self-averaging of dispersion for transport in quasi-periodic random media. J. Phys. Math. Theor., 40(4), 597–610.
[16] Gray, W.G., Miller, C.T., 2013. A generalization of averaging theorems for porous medium analysis. Adv. Water Resour. 62, 227–37.
[17] Gray, W.G., Dye, A.L., McClure, J.E., Pyrak-Nolte, L.J., Miller, C.T., 2015. On the dynamics and kinematics of two-fluid-phase flow in porous media. Water Resour. Res., 51(7), 5365-81.
[18] He, Y., Sykes, J.F., 1996. On the spatial-temporal averaging method for modeling transport in porous media. Transport Porous Media, 22(1), 1–51.
[19] Hornung, U., 1996. Homogenization and porous media (Vol. 6). Springer Science & Business Media.
[20] Illiano, D., Pop, I.S., Radu, F.A., 2020. Iterative schemes for surfactant transport in porous media. Comput. Geosci.
[21] Klofkorn, R., Kr¨oner, D., Ohlberger, M., 2002. Local adaptive methods for convection dominated problems. Int. J. Numer. Meth. Fluid., 40(1-2), 79–91.
[22] Kumar, K., Van Noorden, T.L., Pop, I.S., 2011. Effective dispersion equations for reactive flows involving free boundaries at the microscale. Multiscale Model. Simul., 9(1), 29–58.
[23] Kumar, K., Neuss-Radu, M., Pop, I.S., 2016. Homogenization of a pore scale model for precipitation and dissolution in porous media. IMA J. Appl. Math., 81(5), 877–897.
[24] List, F. and Radu, F.A., 2016. A study on iterative methods for solving Richards’ equation. Comput. Geosci. 20 (2), 341–353.
[25] Mailloux B.J., Fuller, M.E., Rose, G.F., Onstott, T.C., DeFlaun, M.F., Alvarez, E., Hemingway, C., Hallet, R.B., Phelps, T.J., Griffin, T., 2003. A modular injection system, multilevel sampler, and manifold for tracer tests. Ground Water, 41(6), 816–27.
[26] McClure, J.E., Dye, A.L., Miller, C.T., Gray, W.G., 2017. On the consistency of scale among experiments, theory, and simulation. Hydrol. Earth Syst. Sci., 21(2),
[27] Nolen, J., Papanicolaou, G., Pironneau, O., 2008. A framework for adaptive multiscale methods for elliptic problems. Multiscale Model. Simul., 7, 171–196.
[28] Nordbotten, J.M., Celia, M.A., Dahle, H.K., Hassanizadeh, S.M., 2007. Interpretation of macroscale variables in Darcy’s law. Water Resour. Res., 43(8), W08430.
[29] Pop, I.S., Radu, F.A., Knabner, P., 2004. Mixed finite elements for the Richards’ equation: linearization procedure. Journal of computational and applied mathematics 168 (1-2), 365–373.
[30] Porta, G.M., Riva, M., Guadagnini, A., 2012. Upscaling solute transport in porous media in the presence of an irreversible bimolecular reaction. Adv. Water Resour., 35, 151–162.
[31] Schuler, L., Suciu, N., Knabner, P., Attinger, S. 2016. A time dependent mixing model to close PDF equations for transport in heterogeneous aquifers. Adv. Water Resour., 96:55–67.
[32] Schwarze, H., Jaekel, U., Vereecken, H., 2001. Estimation of macrodispersion by different approximation methods for flow and transport in randomly heterogeneous media. Transport Porous Media, 43(2), 265-287.
[33] Schwede, R.L., Cirpka, O.A., Nowak, W., Neuweiler, I., 2008. Impact of sampling volume on the probability density function of steady state concentration. Water Resour. Res., 44(12), W12433.
[34] Skøien, J.O., Bloschl, G., Western, A.W., 2003. Characteristic space scales and timescales in hydrology. Water Resour. Res., 39(10).
[35] Srzic, V., Cvetkovic, V., Andricevic, R., Gotovac, H., 2013. Impact of aquifer heterogeneity structure and local-scale dispersion on solute concentration uncertainty. Water Resour. Res., 49(6), 3712–3728.
[36] Suciu N., 2001. On the connection between the microscopic and macoscopic modeling of the thermodynamic processes. Ed. Univ. Pitesti (in Romanian). English abstract:
[37] Suciu, .N, Radu, F.A., Attinger, S., Sch¨uler, L., Knabner, P., 2015. A Fokker-Planck approach for probability distributions of species concentrations transported in heterogeneous media. J. Comput. Appl. Math., 289, 241–252.
[38] Suciu, N., Schuler, L., Attinger, S., Knabner, P., 2016. Towards a filtered density function approach for reactive transport in groundwater. Adv. Water Resour., 90, 83–98.
[39] Suciu, N., 2019. Diffusion in Random Fields. Applications to Transport in Groundwater. Birkhauser, Cham.
[40] Suciu, N., Illiano, D., Prechtel, A., Radu, F.A., 2021. Global random walk solvers for fully coupled flow and transport in saturated/unsaturated porous media. Adv. Water Resour., 152, 103935.
[41] Suciu, N., Radu, F.A., 2021. Global random walk solvers for reactive transport and biodegradation processes in heterogeneous porous media. arXiv: 2107.08745
[42] Suciu, N., Radu, F.A., 2021. Git repository.
[43] Suciu, N., Radu, F.A., 2021. Git repository.
[44] Teutsch, G., Ptak, T., Schwarz, R., Holder, T., 2000. Ein neues Verfahren zur Quantifizierung der Grundwasserimmission: I. Theoretische Grundlagen. Grundwasser, 4(5):170–175.
[45] Tonina, D., Bellin, A., 2008. Effects of pore-scale dispersion, degree of heterogeneity, sampling size, and source volume on the concentration moments of conservative solutes in heterogeneous formations. Adv. Water Resour., 31(2), 339–354.
[46] Vamos C., Georgescu A., Suciu N., Turcu I., 1996. Balance equations for physical systems with corpuscular structure. Phys. Stat. Mech. Appl., 227(1), 81–92.
[47] Vamos C., Suciu N,, Georgescu, A., 1997. Hydrodynamic equations for one-dimensional systems of inelastic particles. Phys. Rev. E, 55(5), 6277–6280.
[48] Vamos, C., Suciu, N., Peculea, M., 1997. Numerical modelling of the one-dimensional diffusion by random walkers. Rev. Anal. Numer. Theor. Approx., 26(1), 237–247.
[49] Vamos, C., Suciu, N., Blaj, W., 2000. Derivation of one-dimensional hydrodynamic model for stock price evolution. Phys. Stat. Mech. Appl., 287(3), 461–467.
[50] Vamos, C., Suciu, N., Vereecken, H., 2003. Generalized random walk algorithm for the numerical modeling of complex diffusion processes. J. Comput. Phys., 186, 527–544.
[51] Vamos C., 2007. Continuous fields associated to corpuscular systems. Risoprint, Cluj-Napoca (in Romanian). English version:
[52] Wiedemeier, T.H., Rifai, H.S., Newell, C.J., Wilson, J.T., 1999. Natural attenuation of fuels and chlorinated solvents in the subsurface. John Wiley & Sons.
[53] Wright, E. E., Sund, N. L., Richter, D. H., Porta, G. M., Bolster, D., 2021. Upscaling bimolecular reactive transport in highly heterogeneous porous media with the LAgrangian Transport Eulerian Reaction Spatial (LATERS) Markov model. Stoch. Environ. Res. Risk. Assess., 35, 1529–1547.

Related Posts