我们曾经在如下博文中探讨过关于WGHM数据的读取和绘制
https://blog.csdn.net/WH_G_Y/article/details/130514813?spm=1001.2014.3001.5501
最近在进行水文研究时发现,每个水文模型结合水平衡方程计算得到的地下水储量异常(GWSA)时各不相同的,归根结底还是因为每个模型输入水平衡方程的分量种类不一致导致的。简单点说即,A模型TWSA=GWSA+土壤水异常+雪水异常+树冠含水异常,B模型TWSA=GWSA+土壤水异常+雪水异常+树冠含水异常+水库储量异常
查看H. Müller Schmied 2014年题为“Sensitivity of simulated global-scale freshwater fluxes and storages to input data, hydrological model structure, human water use and calibration”论文可知,
TWSA是由绿色框内的分量过程的,是与不是,让我们进行试验验证。
试验思路:
(1)利用TWSA -(绿色框内除黄色Groundwater以外的分量), 绘制全球格网,并绘制区域GWSA时间序列(我们把水平衡方程求得GWSA命名为“GWSA_WB”)
(2)直接读取GWS分量,绘制全球格网和GWSA时间序列
(3)前面两个步骤作差
贴上部分代码
clearvars -except
addpath('E:\RSE\data\WGHM');
filetws='watergap_22d_WFDEI-GPCC_histsoc_tws_monthly_1901_2016.nc4';
info = ncinfo('watergap_22d_WFDEI-GPCC_histsoc_reservoirstor_monthly_1901_2016.nc4');
lon=ncread(filetws,'lon');lat=ncread(filetws,'lat');
tws=rot90(ncread(filetws,'tws'));time=ncread(filetws,'time');
filesoil='watergap_22d_WFDEI-GPCC_histsoc_soilmoist_monthly_1901_2016.nc4';
soil=rot90(ncread(filesoil,'soilmoist'));
filereser='watergap_22d_WFDEI-GPCC_histsoc_reservoirstor_monthly_1901_2016.nc4';
reser=rot90(ncread(filereser,'reservoirstor'));
filecan='watergap_22d_WFDEI-GPCC_histsoc_canopystor_monthly_1901_2016.nc4';
can=rot90(ncread(filecan,'canopystor'));
fileriv='watergap_22d_WFDEI-GPCC_histsoc_riverstor_monthly_1901_2016.nc4';
riv=rot90(ncread(fileriv,'riverstor'));
fileswe='watergap_22d_WFDEI-GPCC_histsoc_swe_monthly_1901_2016.nc4';
swe=rot90(ncread(fileswe,'swe'));
fileglowet='watergap_22d_WFDEI-GPCC_histsoc_glowetlandstor_monthly_1901_2016.nc4';
glowet=rot90(ncread(fileglowet,'glowetlandstor'));
filelowet='watergap_22d_WFDEI-GPCC_histsoc_locwetlandstor_monthly_1901_2016.nc4';
lowet=rot90(ncread(filelowet,'locwetlandstor'));
fileglla='watergap_22d_WFDEI-GPCC_histsoc_glolakestor_monthly_1901_2016.nc4';
glla=rot90(ncread(fileglla,'glolakestor'));
filelola='watergap_22d_WFDEI-GPCC_histsoc_loclakestor_monthly_1901_2016.nc4';
lola=rot90(ncread(filelola,'loclakestor'));
gws_WB=tWS-(soil+reser+can+riv+swe+...
glowet+lowet+glla+lola);
filegroundwstor='E:\RSE\data\WGHM\watergap_22d_WFDEI-GPCC_histsoc_groundwstor_monthly_1901_2016.nc4';
groundwstor=rot90(ncread(filegroundwstor,'groundwstor'));
结果
对比区域的时间序列图
欢迎各位做GRACE时变重力场和水文数据处理的友友私信与我交流,b站频道同步更新哦