I noticed in phreatic methods that the elements with negative pressure and residual saturation show variation of storage (release/capture).

How is this variation computed on the DOI balances?

And how could these volumes released by the residual saturation elements change the entire volume? 

When comparing the volumes released from the conceptual storage with the modeled one, I found a huge discrepancy.