目前時間分數階傳輸模式只被用來分析與地質材料不起吸脫附作用的非反應性示蹤劑實驗，因此本研究的目的是發展可用來分析與地質材料起吸脫附作用的反應性示蹤劑實驗的時間分數階傳輸模式。假設溶質與地質材料有非均衡吸脫附作用(nonequilibrium sorption)和溶質衰減，本模式修改Schumer et al.[2003]的冪律記憶函數(power-law memory function)來處理時間分數階的微分，但是加上時間尺度有效率係數F[Tγ-1]，用以使模式裡每項因次相同；γ為時間分數階階數，0<γ<1。此時間分數階傳輸模式對時間(包括分數階)作拉普拉斯逆轉換法求得數值解，並且使用有限差分隱示法求解作驗證，成功地分析Starr et al.[1985]的一組非反應性(溴離子)實驗得到γ=0.59，F=1.09[d-0.41]與四組反應性 (鍶-85)實驗得到γ=0.9，F值隨著有效縱向流速的變小而減少。反應性實驗γ值大於非反應性實驗γ值說明非均衡吸脫附作用抓住示蹤劑造成溶質質量損失，且不同類型的示蹤劑會有不同γ值。

For non-Fickian transport problems, the current fractional-in-time models can only deal with nonreactive tracers that do not chemically reacted with geological materials. This study develops a fractional-in-time model suitable for reactive tracers, assuming a non-equilibrium sorption process and biological/radioactive decay. The model development adopts the same modeling approach as given by Schumer et al. [2003], while allowing the power-law memory function involved to be modified by including a time-scaling effective rate coefficient, F in [Tγ-1], where 0<γ<1 is the fractional order of the time derivative in the governing equation. This modification renders a uniform dimensionality for the resultant fractional-in-time governing equation. Solutions to the model developed are determined using the Laplace transform with respect to time and a numerical Laplace inversion technique. The solutions so obtained are verified using a finite-difference model based on Grünwald series for the fractional-in-time term. There are four sets of reactive 85Sr tracer test data from a sand box consisted of a sand layer sandwiched between two silt layers. They cannot be analyzed by the conventional advective-dispersion equation models, but can be successfully analyzed using the fractional-in-time model developed herein. For the bromide breakthrough, it is found thatγ=0.59 and F=1.09 [d-0.41]. For the 85Sr breakthroughs, it is found that γ=0.9 which is larger than that of the bromide because of extra solute mass losses due to non-equilibrium sorption and radioactive decay, and F changes from 15.2~17.5 [d-0.1] under the influence of the flow velocity of the tracer tests. The fractional-in-time model developed appears to be useful for non-Fickian solute transport for both nonreactive and reactive tracers.

Berkowitz, B., and H. Scher (1995), On characterization of anomalous dispersion in porous and fractured media, Water Resour. Res., 31(6), 1461-1466.
Berkowitz, B., S. Emmanuel, and H. Scher (2008), Non-Fickian transport and multiple-rate mass transfer in porous media, Water Resour. Res., 44, W03402, doi: 10.1029/2007WR0059068.
Brusseau, M. L. and P. S. C. Rao (1989), Modeling the transport of solutes influenced by multiprocess nonequilibrium, Water Resour. Res., 25(9), 1971-1988, 1989.
Chen, W., and R. J. Wagenet (1995), Solute transport in porous media with sorption-site heterogeneity, Environ. Sci. Technol., 29(11), 2725-2734.
Connaughton, D. F., J. R. Stedinger, L. W. Lion, and M. L. Schuler (1993), Desorption of time-varying desorption kinetics: Release of naphthalene from contaminated soils, Environ. Sci. Technol., 27(12), 2397-2403.
Cortis, A., and B. Berkowitz (2004), Anomalous Transport in "Classical" Soil and Sand Columns, Soil Sci. Soc. Am. J. 68:1539-1548.
Culver, T. B., S. P. Hallisey, D. Sahoo, J. J. Deitsch, and J. A., Smith (1997), Modeling the desorption of organic contaminants from long-term contaminated soil using distributed mass transfer rates, Environ. Sci. Technol., 31(6), 1581-1588.
Dentz, M., and B. Berkowitz (2003), Transport behavior of a passive solute in continuous time random walks and multirate mass transfer, Water Resour. Res., 39 (5), 1111, doi: 10.1029/2001WR001163.
Farrell, J. and M. Reinhard (1994), Desorption of halogenated organics from model solids, sediments, and soil under unsaturated conditions, 2. Kinetics. Environ. Sci. Technol., 28(1), 63-72.
Fomin, S., V. Chugunov, and T. Hashida (2005), The effect of non-Fickian diffusion into surrounding rocks on contaminant transport in a fractured porous aquifer, Proc. R. Soc. A, 461, 2923-2939, doi: 10.1098/rapa.2005.1487.
Gorenflo, R., and F. Mainardi (1998), Random walk models for space-fractional diffusion processes, Fract. Calculus. Appl. Anal. 1, 167-191.
Haggerty, R., and S. M. Gorelick (1995), Multiple-rate mass transfer for modeling diffusion and surface reactions in media with pore-scale heterogeneity, Water Resour. Res., 31(10), 2383-2400.
Haggerty, R., S. W. Fleming, L. C. Meigs, and S A. McKenna (2001), Tracer tests in a fractured dolomite 2. Analysis of mass transfer in single-well injection-withdrawal tests, Water Resour. Res., 37(5), 1129 – 1142.
Hollenbeck, K. J., C. F. Harvey, R. Haggerty, and C. J. Werth (1999), A method for estimating distributions of mass transfer rate coefficients with application to purging and batch experiments, J Contam Hydrol., 37(3-4):367-88.
Jespersen, S., R. Metzler, and H. C. Fogedby (1999), Lévy flights in external force fields: Langevin and fractional Fokker-Planck equations and their solutions, Phys. Rev. 59(3):2736-2745.
Jiang, X. Y., and M. Y. Xu (2004), Fractional model and solution of pollution sources concentration distribution, J. Shan. Univ. 39(3):37-41.
Liu, F., P. Zhuang, V. Anh, and I. Turner (2006), A fractional-order implicit difference approximation for the space-time fractional diffusion equation, ANZIAM J. 47 C48-C68. EMAC 2005.
Mainardi, F., G. Pagnini, and R. K. Saxena (2005), Fox H function in fractional diffusion, J. Comput. Appl. Math. 178:321-331.
McKenna, S. A., L. C. Meigs, and R. Haggerty (2001), Tracer tests in a fractured dolomite: 3. Double-porosity, multiple-rate mass transfer processes in convergent flow tracer tests, Water Resour. Res., 37(5), 1143-1154, 2001.
Meerschaert, M. M., and C. Tadjeran (2004), Finite difference approximations for fractional advection-dispersion flow equations, J. Comput. Appl. Math. 172:65-77.
Neretnieks, I. (1980), Diffusion in the rock matrix: an important factor in radionuclide retardation, J. Geophys. Res., 85, 4379-4397.
Pedit, J. A., and C. T. Miller (1994), Heterogeneous sorption processes in subsurface systems, 1, Model formulations and applications, Environ. Sci. Technol., 28(12), 2094-2104.
Pignatello, J., and B. Xing (1995), Mechanisms of slow sorption of organic chemicals to natural particles, Environ. Sci. Technol., 30(1), 1-11.
Rahman, M. M., R. Liedl, and P. Grathwohl (2004), Sorption kinetics during macropore transport of organic contaminants in soils: laboratory experiments and analytical modeling, Water Resour. Res., 40, W01503. (doi:10.1029/2002WR001946)
Scher, H., and M. Lax (1973), Stochastic transport in a disordered solid. I. theory, Phys. Rev. B 7, 4491-4502.
Schumer R., D. A. Benson, M. M., Meerschaert, and B. Baeumer (2003), Fractal mobile/immobile solute transport, Water Resour. Res., 39(10), 1296, doi: 10.1029/2003WR002141.
Starr, R. C., R. W. Gillham, and E.A. Sudicky (1985), Experimental investigation of solute transport in stratified porous media: 2. The reactive case, Water Resour. Res., 21(7), 1043-1050.
Wylie, C. R. and L. C. Barrett (1982), Advanced Engineering Mathematics, 5th edition, McGraw-Hill Inc.
Xiong , Y., G. Huang, and Q. Huang (2006), Modeling solute transport in one-dimensional homogeneous and heterogeneous soil columns with continuous time random walk, J. Contam Hydrol., 86:163-175.