Root solute uptake in dual permeability model

A discussion forum for old Hydrus-1D users.
Post Reply
YOON
Posts: 14
Joined: Sun Jul 23, 2023 1:51 pm
Location: Korea

Root solute uptake in dual permeability model

Post by YOON » Tue Jan 09, 2024 7:40 am

Dear all,
I am currently working on a modeling study using dual-permeability model(physical non-equilibrium) and have a question.
I am attempting to estimate the solute concentration of the root water uptake in the dual-perm model. However, I have noticed some discrepancies between the results of cvRoot and sum(cvRoot), where the values in cvRoot seem greater than what they should be.
I have tried to estimate sum(cvRoot) using cvRoot as the cumulative sum of the average cvRoot between two-time steps multiplied by the length of the time step, and it seems that the estimate is larger than sum(cvRoot) in general.
I also found a relationship between the {estimated sum(cvRoot) for the final time step} / {the reported sum(cvRoot) in solute1.out file for the final time step} and w (the ratio of fracture) based on the results gathered from numerous model runs using different parameters. It appears that these two variables exhibit a certain clear relationship as depicted in the attached figure. I would appreciate any insights or explanations regarding this phenomenon.
Image
<img src="https://drive.google.com/file/d/1CzZYL- ... drive_link">

Best regards,
Yoon.

YOON
Posts: 14
Joined: Sun Jul 23, 2023 1:51 pm
Location: Korea

Re: Root solute uptake in dual permeability model

Post by YOON » Tue Mar 05, 2024 12:49 am

I'm writing again since I haven't received a response yet.
If anyone could offer assistance, it would be greatly appreciated.
Thank you.

Jirka
Posts: 6258
Joined: Sat Mar 16, 2002 3:47 pm
Location: USA
Location: Riverside, CA

Re: Root solute uptake in dual permeability model

Post by Jirka » Tue Mar 05, 2024 3:43 am

One of the assumptions in the dual permeability model is that root water uptake occurs only from the matrix (by setting qTopF=0; a fraction of uptake from fractures). At each time step, the code calculates RWU from both domains (even though it is zero from the fracture domain), and integrates cumulative RWU.
For water:
qTopF=0. ! =0: only matrix
TPotM=TPot*(1.-qTopF)/(1.-wF) ! potential transpiration
TPotF=TPot*qTopF/wF

vRootT=vRoot*(1.-wF)+vRootF*wF ! actual total RWU=RWU from matrix + RWU from fractures
CumQ(4)=CumQ(4)+vRootT*dt ! cumulative total RWU
For solute :
CumCh(5,jS)=CumCh(5,jS)+cvChR(jS)*dt*(1.-wF) ! uptake from matrix
CumChF(5,jS)=CumChF(5,jS)+cvChRF(jS)*dt*wF ! uptake from fractures

J.

Post Reply