Issue in creating LISFLOOD input files for Tapi basin in India

[diksha922]

Sir/ Madam,

I want to use LISFLOOD model for hydrological modelling over Tapi basin in India, using the rainfall data from Indian Meteorological Department at 0.25 degree and temperature at 1 degree, on a daily time step. However, I am facing issues in creating the meteorological inputs files using these inputs.

Further, I am not getting how to create the other input files such as topography, landuse, soil map etc. required for running the hydrological model over Tapi basin.

Please help me in creating the input files for the study area for running the hydrological model.

[StefaniaGrimaldi]

Dear @diksha922,

thank you for your interest in using OS LISFLOOD for hydrological modelling in the Tapi basin.

OS LISFLOOD requires the following meteorological inputs: total precipitation, average temperature, reference evapotranspiration.

All the input meteorological forcings and static maps (topography, land use, soil maps, ...) must have the same spatial extent, the same reference system, the same pixel size. This is a strict requirement for model implementation: more details can be found here https://ec-jrc.github.io/lisflood-code/4_Static-Maps-introduction/

Therefore, the first step for the preparation of the input meteorological forcings and of static maps consists in defining spatial extent, reference system, and pixel size for your application.

Your rainfall data and temperature data have different pixel size, and they require re-gridding prior to their use within OS LISFLOOD. Depending on the format and structure (e.g. GRIB) of your rainfall and temperature dataset, you might be interested in checking the functionalities of the open-source tool pyg2p: https://github.com/ec-jrc/pyg2p

As a follow up of this previous conversation https://github.com/ec-jrc/lisflood-lisvap/issues/56, we shared with you the mask map and the ERA5 meteorological variables for the computation of reference evapotranspiration values for the Tapi basin using OS LISVAP. The mask map and all the meteorological variables had 3arcmin (or 0.05 degrees) resolution.

In case you decided to use 3arcmin (or 0.05 degrees) resolution for your OS LISFLOOD model set up, you might be interested in the openly available LISFLOOD static maps dataset at 3arcmin (or 0.05 degrees) resolution, for the global domain: https://data.jrc.ec.europa.eu/dataset/68050d73-9c06-499c-a441-dc5053cb0c86. This manuscript https://egusphere.copernicus.org/preprints/2023/egusphere-2023-1306/ provides a detailed description of the raw datasets and of the protocols used for the generation of each map of the openly available dataset.

Clearly, users can generate their own input maps dataset: the technical requirements and the use of each static map within LISFLOOD are explained here: https://ec-jrc.github.io/lisflood-code/4_Static-Maps-introduction/.

We hope that our answer helps, please do not hesitate to get in touch for further questions, kind regards, on behalf of the developers team, Stefania

[diksha922]

Mam,

I have generated the potential evapotranspiration using the input file provided at: https://github.com/ec-jrc/lisflood-lisvap/issues/56, at 3 arc resolution. I have rainfall available at 0.25 degree and temperature at 1 degree.

I am not getting how to convert all the input files to same resolution. Please help me for getting the same, may be some reference examples may help.

Also, it will be helpful if some guide/ reference may be provided to generate LISFLOOD compatible inputs from .mat file of rainfall and temperature.

[StefaniaGrimaldi]

Dear @diksha922,

what is the format of your rainfall and temperature dataset? Different solutions are possible depending on the format of your datasets.

Kind regards, on behalf of the developers team, Stefania

[diksha922]

Madam,

The rainfall data is available as gridded dataset having three dimensions (X, Y, value) with 0.25 degree for rainfall and 1 degree for temperature, in .mat format.

[StefaniaGrimaldi]

Dear @diksha922,

thanks for providing more details.

Unfortunately, we cannot provide support for the conversion and re-gridding of matlab files.

However, online blogs propose various solutions. A quick internet search allowed to find a few examples: https://it.mathworks.com/matlabcentral/answers/565163-write-a-mat-file-into-netcdf-file and https://gdal.org/index.html. Please note that these are only examples. These are neither recommendations nor guidelines. We do not work with .mat files in our operational protocols.

When defining your working protocol, you should use the template map used in LISVAP to guarantee the consistency between LISFLOOD forcings files.

Kind regards, on behalf of the developers team, Stefania