Hi Gaurav,
I am not sure that using runoff as a proxy for river and lake storage is a good strategy. In the GLDAS models, runoff, along with precipitation and evapotranspiration, is a flux out of (or into) a grid cell. To calculate groundwater storage change (dGW), you are using an equation like this:
dTWS = d(SM + SWE + CAN + GW + Rivers + Lakes)
Let’s say we ignore SWE and CAN (canopy storage) for now and combine rivers and lakes into one variable called (SURF).
dTWS = d(SM + GW + SURF)
GRACE measures dTWS and GLDAS gives you dSM. But there are still two unknowns (dGW and dSURF).
If you want to bring in runoff, that is a flux, not a component of water storage. It relates to total water storage as follows:
dTWS = P – ET – R
But this does not tell us any new information. We already know all four elements of this equation (dTWS from GRACE and P, ET, and R from GLDAS). In other words, runoff tells us about TWS change as a whole, but it does not tell us about the components of TWS change (GW vs. SWE vs. rivers and lakes, …).
Ideally, you’d have a land surface model that explicitly represents rivers and lakes with a runoff routing model, but none of the GLDAS models do this. That said, in some areas of the world, you may be able to assume surface water storage is small compared to TWS (see Getirana et al. 2017, GRL)
https://agupubs.onlinelibrary.wiley.com/doi/full/10.1002/2017GL074684 .
Also, keep in mind that soil moisture variations are sometimes overestimated by land surface models, so that when you subtract modeled soil moisture from GRACE TWS, the remaining (groundwater plus surface water) variations are underestimated.
You may be interested in the following:
• Monthly TWS anomalies from GLDAS published at PODAAC:
https://podaac.jpl.nasa.gov/dataset/TELLUS_GLDAS-NOAH-3.3_TWS-ANOMALY_MONTHLY# (also lists some papers citing the data, such as Lv et al. 2020, J. of Hydromet.).
• TWS calculation example from GLDAS Noah outputs:
https://grace.jpl.nasa.gov/data/get-data/land-water-content/
• Panda and Wahr, 2016, WRR, where they estimate groundwater storage by subtracting GLDAS NOAH SM from GRACE TWS and their justification citing Tewari et al. 2009 in GRL.
https://agupubs.onlinelibrary.wiley.com/doi/10.1002/2015WR017797
• The GLDAS Readme file, see (4) in Section 3.5 Data Interpretation. It sheds some light on how TWS is calculated in the CLSM and Noah models.
Hope this helps!
Jacob