# GR4J
**GR4J** (modèle du Génie Rural à 4 paramètres Journalier, French for "rural engineering model with 4 daily parameters") is a daily lumped conceptual rainfall-runoff model that relates catchment runoff to rainfall and potential evapotranspiration using four calibrated parameters and a fixed structure of two non-linear stores and two unit hydrographs.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) • [2](https://kalix.org/docs/gr4j-technical/) It was developed at Cemagref (now INRAE) in France and is widely used in hydrology for its parsimony, computational efficiency, and effectiveness across a range of climates.[2](https://kalix.org/docs/gr4j-technical/) • [3](https://doi.org/10.1016/j.jhydrol.2026.135366)
| Key fact | Detail |
|---|---|
| Type | Daily lumped conceptual rainfall-runoff model of the soil moisture accounting family[4](https://www.sciencedirect.com/science/article/pii/S0022169403002257) |
| Parameters | x1 production store capacity (mm), x2 water exchange coefficient, x3 routing store capacity (mm), x4 unit hydrograph time parameter (days)[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) |
| Default values | x1 = 350 mm, x2 = 0, x3 = 40 mm, x4 = 0.5 days; ranges 1–1500, −10.0 to 5.0, 1–500, and 0.5–4.0 respectively[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) |
| Reference version | Perrin, Michel, and Andréassian, *Journal of Hydrology*, 25 August 2003, vol. 279, pp. 275–289[4](https://www.sciencedirect.com/science/article/pii/S0022169403002257) |
| Routing split | 90% of effective rainfall routed through UH1 and a non-linear routing store, 10% through UH2[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) |
| Snow option | No snow routine of its own; commonly coupled with the two-parameter degree-day snow model CemaNeige[5](https://hess.copernicus.org/articles/29/2811/2025/index.html) |
| Software | Implemented in the R package airGR (Fortran model cores), the Source platform, and other tools[6](https://cran.r-project.org/web/packages/airGR/vignettes/V01_get_started.html) • [7](https://www.tandfonline.com/doi/abs/10.1080/09715010.2018.1556124) |
## History
Development of the GR family was initiated by Claude Michel at the beginning of the 1980s at Cemagref, a public research institute in France, starting from the nine-parameter CREC model of the University of Montpellier hydrology laboratory.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) • [8](https://webgr.inrae.fr/outils/modeles-hydrologiques/bref-historique) A simplified two-reservoir structure with two parameters (GR2) was applied to the Orgeval basin in 1983.[8](https://webgr.inrae.fr/outils/modeles-hydrologiques/bref-historique) Edijatno and Michel proposed the three-parameter daily version GR3 in 1989, tested on a sample of 110 French catchments.[8](https://webgr.inrae.fr/outils/modeles-hydrologiques/bref-historique) Nascimento (1995) added a fourth parameter governing groundwater exchanges, and Edijatno et al. (1999) proposed a further three-parameter version with a fixed production store capacity.[8](https://webgr.inrae.fr/outils/modeles-hydrologiques/bref-historique)
The reference four-parameter version, GR4J, was proposed by Charles Perrin (2000) and detailed by Perrin et al. (2003), introducing a percolation term from the production store that improved the simulation of low flows.[8](https://webgr.inrae.fr/outils/modeles-hydrologiques/bref-historique) • [1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) The modeling approach is mainly empirical: model structures are searched for and evaluated on large samples of catchments, with the objective of a general, efficient, and robust model.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) Successive extensions followed: the five-parameter GR5J with an improved exchange function (Le Moine, 2008) and the six-parameter GR6J with a second routing store for better low flows (Pushpalatha, 2013).[8](https://webgr.inrae.fr/outils/modeles-hydrologiques/bref-historique)
## Structure
GR4J operates at the catchment scale with a daily time step; all water quantities are expressed in millimeters by dividing volumes by catchment area.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) Inputs are catchment-average rainfall P and potential evapotranspiration E, which must be continuous and overlapping daily series.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG)
**Production store.** The first operation computes net rainfall Pn or net evapotranspiration En by subtracting E from P, as if through an interception storage of zero capacity.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) A part Ps of the net rainfall fills the production store, a soil moisture accounting reservoir of capacity x1 whose filling and evaporation functions follow parabolic forms with terms in (S/x1)².[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) A percolation leakage Perc, calculated as a power function of the store content, moves water to the routing functions; this term contributes little to total streamflow and matters mainly for low-flow simulation.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG)
**Routing.** The water reaching the routing functions, Pr, is split in a fixed proportion: 90% passes through a unit hydrograph UH1 and then a non-linear routing store of capacity x3, while the remaining 10% passes through a second unit hydrograph UH2.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) Both unit hydrographs share the time parameter x4 in days; UH1 has a time base of x4 days and UH2 a time base of 2·x4 days, and x4 must exceed 0.5 days.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) The routing store can simulate long streamflow recessions, and its capacity x3 is described as a "one day ahead maximum capacity" because the store level cannot exceed x3 at the end of a time step even when the input exceeds the saturation deficit.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG)
**Groundwater exchange.** A groundwater exchange term F acts on both flow components as a function of the routing store level R, its reference capacity x3, and the exchange coefficient x2.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) A positive x2 represents water imports from outside the catchment, a negative x2 exports to deeper aquifers, and zero means no exchange; in absolute value F cannot exceed x2.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) • [9](https://www.researchgate.net/publication/337773134_THE_IMPROVEMENT_OF_GR4J_MODEL_PARAMETER_TO_ESTIMATE_UNIT_HYDROGRAPH) Total streamflow is the sum of the routed component Qr and the direct component Qd.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG)
Several fixed coefficients appear in the equations, such as a power of 4 in the percolation and unit hydrograph computations, the 90%–10% split, a 2.5 exponent in the unit hydrographs, and a 2.25 coefficient in the percolation function. These were chosen as the values yielding the best results across many test conditions; leaving them free did not significantly improve results and added unhelpful complexity.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG)
## Calibration and software
Given the small number of parameters, simple optimization algorithms can generally identify parameter values yielding satisfactory results; the choice of objective function depends on the user's objectives, and a one-year warm-up is recommended at the start of each simulation to set appropriate initial states.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG) Median parameter values over large catchment samples, with approximate 80% confidence intervals, serve as initial parameter sets.[1](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG)
The R package airGR brings the GR models developed at INRAE's Catchment Hydrology Research Group into R, with each model core coded in Fortran for low computational time.[6](https://cran.r-project.org/web/packages/airGR/vignettes/V01_get_started.html) Its `RunModel_GR4J()` function implements the Perrin et al. (2003) model; by default the production and routing stores are initialized at 30% and 50% of capacity, and a warm-up period equal to or longer than one year is advised.[6](https://cran.r-project.org/web/packages/airGR/vignettes/V01_get_started.html) GR4J also runs in the Source modeling platform used in Australia.[7](https://www.tandfonline.com/doi/abs/10.1080/09715010.2018.1556124)
## Variants and temporal resolutions
The daily GR4J has an hourly counterpart, GR4H, and the model can be run at other temporal resolutions provided its time-dependent parameters are adjusted (Ficchì et al., 2016); a state-space formulation of the GR4 model has also been published (Santos et al., 2018).[10](https://unifhy-org.github.io/unifhy/science/surfacelayer/unifhycontrib.gr4.SurfaceLayerComponent.html) In snowy climates GR4J is coupled with CemaNeige, a two-parameter degree-day snow accumulation and melt model, as in the GR4J-CN configuration.[5](https://hess.copernicus.org/articles/29/2811/2025/index.html) • [11](https://www.sciencedirect.com/science/article/pii/S2214581823000460) In the Chile and Switzerland intercomparison of Muñoz-Castro et al., GR4J, GR5J, and GR6J coupled with CemaNeige carry 6, 7, and 8 parameters respectively.[12](https://hess.copernicus.org/articles/30/825/2026/hess-30-825-2026.pdf)
## Performance and applications
The 2003 reference paper assessed GR4J on a large sample of catchments and found its performance statistically among the best of the compared rainfall-runoff models, with clear gains over the previous three-parameter version in low-flow simulation; the tests indicated that a four-parameter structure corresponds to the maximum level of complexity the model could afford, since adding more free parameters did not bring significant improvements.[4](https://www.sciencedirect.com/science/article/pii/S0022169403002257) Because calibrated rainfall-runoff models match observed streamflow closely, GR4J and similar lumped models are used as benchmarks against land surface and global hydrological models; in one evaluation the four models were assessed against runoff from 644 unregulated catchments and actual evapotranspiration from 98 eddy covariance towers worldwide.[13](https://doi.org/10.1175/JHM-D-15-0107.1)
Operational and regionalization uses are extensive. GR4J was applied to generate continuous daily streamflow from 1993 onwards for 405 gauged and 771 ungauged coastal catchments across Australia, covering 2,549,000 km² and 835,000 km² respectively.[14](https://www.mdpi.com/2073-4441/16/4/527) In ungauged-basin studies, GR4J parameters are regionalized from gauged donor catchments; in the Nepal Himalaya, GR4J-CN achieved median Nash-Sutcliffe efficiency of 0.76 in calibration and 0.74 in validation, with physical similarity the most robust of four regionalization methods tested.[11](https://www.sciencedirect.com/science/article/pii/S2214581823000460) In Ireland, regionalized GR4J performed best for high flows while GR6J was best for low flows.[15](https://www.sciencedirect.com/science/article/pii/S2214581821000884) Hybridizations replace parts of the model with machine learning: DeepGR4J substitutes deep learning models (LSTM and CNN) for the routing storage processes, and QDeepGR4J builds a quantile-based ensemble for extreme flow prediction with uncertainty quantification.[16](https://doi.org/10.1016%2Fj.envsoft.2023.105831) • [17](https://arxiv.org/abs/2510.05453)
## Limitations
GR4J's main structural weaknesses follow from its lumped, conceptual design. It has no explicit groundwater module, and evaluations in the Sahelian Oti sub-basin found highly negative low-flow criteria (KGE on 1/Q of about −0.65 to −0.71), showing an inability to reproduce low-flow discharges; peak floods were also often underestimated, reflecting the limits of the single-path routing reservoir for rapid surface runoff.[18](https://www.mdpi.com/2073-4441/18/12/1501) The lumped formulation inadequately describes [runoff generation](https://www.edgepedia.org/runoff-generation) and flow concentration in large catchments, and conceptual lumped models are considered far from able to assess impacts of land use or forest variation.[19](https://hess.copernicus.org/preprints/hess-2021-646/hess-2021-646-ATC2.pdf) In spatial transferability tests across climates and underlying surface conditions, the more complex distributed models (XAJ and HBV) outperformed GR4J.[20](https://www.sciencedirect.com/science/article/abs/pii/S0022169420307368) On the theoretical side, a recent appraisal showed that GR4J's soil water model is mathematically identical to the Probability Distributed Model (PDM), employing a fixed probability density function resembling a log-normal distribution, and that its empirical unit hydrograph can be moment-matched with a diffusion wave model derived from the Saint-Venant equations.[21](https://www.sciencedirect.com/science/article/pii/S0022169425017330)
## References
1. [GR4J – SRG, Source User Guide 4.11, eWater](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG)
2. [GR4J technical reference, Kalix](https://kalix.org/docs/gr4j-technical/)
3. [Tracking rainfall contribution in a conceptual runoff model, Journal of Hydrology](https://doi.org/10.1016/j.jhydrol.2026.135366)
4. [Perrin, C., Michel, C., Andréassian, V. (2003). Improvement of a parsimonious model for streamflow simulation. Journal of Hydrology 279, 275–289](https://www.sciencedirect.com/science/article/pii/S0022169403002257)
5. [Assessing the adequacy of traditional hydrological models for climate change impact studies, HESS 29, 2811 (2025)](https://hess.copernicus.org/articles/29/2811/2025/index.html)
6. [Get Started with airGR, CRAN vignette](https://cran.r-project.org/web/packages/airGR/vignettes/V01_get_started.html)
7. [Performance evaluation of conceptual rainfall-runoff models GR4J and AWBM, ISH Journal of Hydraulic Engineering](https://www.tandfonline.com/doi/abs/10.1080/09715010.2018.1556124)
8. [Bref historique, webGR, INRAE](https://webgr.inrae.fr/outils/modeles-hydrologiques/bref-historique)
9. [The improvement of GR4J model parameter to estimate unit hydrograph, ResearchGate](https://www.researchgate.net/publication/337773134_THE_IMPROVEMENT_OF_GR4J_MODEL_PARAMETER_TO_ESTIMATE_UNIT_HYDROGRAPH)
10. [unifhycontrib.gr4 documentation](https://unifhy-org.github.io/unifhy/science/surfacelayer/unifhycontrib.gr4.SurfaceLayerComponent.html)
11. [Comparative performance of regionalization methods for model parameterization in ungauged Himalayan watersheds, Journal of Hydrology: Regional Studies](https://www.sciencedirect.com/science/article/pii/S2214581823000460)
12. [Muñoz-Castro, E. et al. Simulating compound streamflow extremes, HESS 30, 825 (2026)](https://hess.copernicus.org/articles/30/825/2026/hess-30-825-2026.pdf)
13. [Evaluating Regional and Global Hydrological Models against Streamflow and Evapotranspiration Measurements, J. Hydrometeorology](https://doi.org/10.1175/JHM-D-15-0107.1)
14. [Simulation of Gauged and Ungauged Streamflow of Coastal Catchments across Australia, Water 16(4), 527 (2024)](https://www.mdpi.com/2073-4441/16/4/527)
15. [Regionalization of hydrological models for flow estimation in ungauged catchments in Ireland, Journal of Hydrology: Regional Studies](https://www.sciencedirect.com/science/article/pii/S2214581821000884)
16. [DeepGR4J: A deep learning hybridization approach for conceptual rainfall-runoff modelling, Environmental Modelling & Software](https://doi.org/10.1016%2Fj.envsoft.2023.105831)
17. [QDeepGR4J: Quantile-based ensemble of deep learning and GR4J hybrid rainfall-runoff models, arXiv](https://arxiv.org/abs/2510.05453)
18. [Climate Simulation and Projection of Rainfall–Runoff Dynamics Using the GR4J Model in the Oti Sub-Basin, Water 18(12), 1501](https://www.mdpi.com/2073-4441/18/12/1501)
19. [HESS preprint discussion of GR4J limitations](https://hess.copernicus.org/preprints/hess-2021-646/hess-2021-646-ATC2.pdf)
20. [Temporal and spatial transferabilities of hydrological models under different climates and underlying surface conditions, Journal of Hydrology](https://www.sciencedirect.com/science/article/abs/pii/S0022169420307368)
21. [A theoretical appraisal of the GR4J rainfall-runoff modelling framework, Journal of Hydrology](https://www.sciencedirect.com/science/article/pii/S0022169425017330)
---
*Topic: Encyclopedia › Physical world and mathematics › Earth sciences › Hydrology and ocean science › Hydrology › Hydrological modeling and software*
*Initially written Oct 15, 2026*
*Copyright 2026 EdgeChat AI, a subsidiary of Biostate AI.*
License: Edgepedia Community License 1.0, https://www.edgepedia.org/edgepedia/license