Feature request: use of non gauged stations in the network

The idea is to be able to integrate GR models in the network with parameters dependant of neighbour sub-basins.

Ungauged stations

Todo:

• Add SubGraph containing ungauged stations in plot.GRiwrm
• Adapt CreateInputsModel.GRiwrm for handling ungauged nodes
• Create a function that creates a RunModel function including the modelling of ungauged nodes for run in Calibration
• Modify Calibration.GRiwrmInputsModel for calibrate the model above and set the calibration parameters for each sub-basin

changed milestone to %v0.7.0

added suggestion label

Use a calibration based on:

Lobligeois, F., 2014. Mieux connaître la distribution spatiale des pluies améliore-t-il la modélisation des crues ? Diagnostic sur 181 bassins versants français (phdthesis). AgroParisTech. https://pastel.archives-ouvertes.fr/tel-01134990/document

added feature label

Specification

I propose to use the solution M2 above which consists in having 2 separated GR models. The latter have the same parameters except for the unit hydrogram parameter which is denoted as:

$$x_{4i} = \left( \dfrac{S_i}{S_{BV}} \right) ^ {0.3} X_4$$

With $S_i$ the area of the intermediate basin delimited by the i^th non gauged station and $S_{BV}$ the area of the gauged intermediate basin (remember that other gauged stations can be located upstream).

We allow multiple upstream non-gauged stations and these ones can eventually be queued as below:

graph LR
    subgraph BV gauged with X4 parameter
    BVD(BV 0 <br/> x4_0)
    BVU1(BV up 1 <br/> x4_1)
    dots(BV up ...<br/> x4_...)
    dots2(BV up ...<br/> x4_...)
    BVUI(BV up n <br/> x4_n)
    BVU1-->BVD
    dots2-->dots-->BVD
    BVUI-->BVD
    end

Outstanding issue

Question for @guillaume.thirel: the document above only mentions GR5J model. Does the formula applicable for GR4J and GR6J models ? What about annual, monthly, and hourly time step models ?

Implementation

The use of a non-gauged station only implies modification of the Calibration process.

Configuration of ungauged stations

First, we need a parameter indicating the a station is non gauged. I think that this one should be defined in CreateCalibOptions.GRiwrmInputsModel through a parameter UngaugedIds containing a character vector of the ungauged station identifiers.

Adaptation of the Calibration process

1. the loop over gauged stations should skip the ungauged stations
2. for each stations, we search for direct link with upstream ungauged stations
3. a custom RunModel function running all the model is injected into the Calibration function
4. the calibrated parameters are copied and converted to all the gauged station and upstream ungauged ones

Outstanding issue

• Use a featured RunModel for a given ensemble of one gauged and ungauged stations and make sure that Calibration process it as an original GR model (for instance "RunModelGR5J_Cemaneige").

Question for @guillaume.thirel: the document above only mentions GR5J model. Does the formula applicable for GR4J and GR6J models ? What about annual, monthly, and hourly time step models ?

The formula is also applied to GR4J (see Param_Sets_GR4J) and to hourly models in GRSD (Fortran edition). By extension, we can surely apply it to any daily or hourly GR model. GR2M and GR1A do not present any unit hydrograph, so we don't need this formula for them.

the loop over gauged stations should skip the ungauged stations

Except if useUpstreamQsim, right? In that case, we need to perform a RunModel, to give upstream streamflows for the next downstream station Well, maybe we can have a call, I'm not sure what your plans are and what is the best solution (but the good thing is that I think this is not very difficult to implement).

created merge request !40 (merged) to address this issue

mentioned in merge request !40 (merged)

created merge request !44 (closed) to address this issue

mentioned in merge request !44 (closed)

changed the description

mentioned in commit 2458ca57

marked the checklist item Add SubGraph containing ungauged stations in plot.GRiwrm as completed

mentioned in commit 1ef804d6

marked the checklist item Adapt CreateInputsModel.GRiwrm for handling ungauged nodes as completed

mentioned in commit 27229352

mentioned in issue #87

mentioned in commit 763adf21

mentioned in commit e4d5e384

mentioned in commit d72a0b62

mentioned in commit ece31a8d

mentioned in commit b9a36d46

mentioned in commit 35825f1b

mentioned in commit ea15a215

mentioned in commit 9009e720

marked the checklist item Create a function that creates a RunModel function including the modelling of ungauged nodes for run in Calibration as completed

marked the checklist item Modify Calibration.GRiwrmInputsModel for calibrate the model above and set the calibration parameters for each sub-basin as completed

mentioned in commit 0b3a0d26

mentioned in commit fbb87832

closed via merge request !40 (merged)

mentioned in commit c2931db8

mentioned in issue #92 (closed)

mentioned in commit 6ea126b4

mentioned in commit 5f56aa0e