All data are 2D matrixes (T*P), where T is the length of time series, P is the number of spatial grid points covers only continental regions. Landmask_50km.mat and Lat_Lon_50km.mat can be used to have a 360 * 720 grid of water stotage changes over the globe.