{ "id": "epen24w97w", "slug": "gr4j", "title": "GR4J", "updated": "2026-10-15", "topic_path": [ { "id": "physical", "label": "Physical world and mathematics", "api_url": "https://www.edgepedia.org/api/v1/topics/physical" }, { "id": "physical.earth", "label": "Earth sciences", "api_url": "https://www.edgepedia.org/api/v1/topics/physical.earth" }, { "id": "physical.earth.water", "label": "Hydrology and ocean science", "api_url": "https://www.edgepedia.org/api/v1/topics/physical.earth.water" }, { "id": "physical.earth.water.hydrology", "label": "Hydrology", "api_url": "https://www.edgepedia.org/api/v1/topics/physical.earth.water.hydrology" }, { "id": "physical.earth.water.hydrology.modeling-software", "label": "Hydrological modeling and software", "api_url": "https://www.edgepedia.org/api/v1/topics/physical.earth.water.hydrology.modeling-software" } ], "geo": [ { "id": "geo.weu.t1946.physical.earth.water", "label": "Western Europe · 1946 to 2000: Hydrology and ocean science", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.weu.t1946.physical.earth.water", "path": [ { "id": "geo.weu", "label": "Western Europe", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.weu" }, { "id": "geo.weu.t1946", "label": "Western Europe · 1946 to 2000", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.weu.t1946" }, { "id": "geo.weu.t1946.physical", "label": "Physical world and mathematics", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.weu.t1946.physical" }, { "id": "geo.weu.t1946.physical.earth", "label": "Earth sciences", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.weu.t1946.physical.earth" }, { "id": "geo.weu.t1946.physical.earth.water", "label": "Hydrology and ocean science", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.weu.t1946.physical.earth.water" } ] }, { "id": "geo.weu.t2001.physical.earth", "label": "Western Europe · 2001 to 2020: Earth sciences", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.weu.t2001.physical.earth", "path": [ { "id": "geo.weu", "label": "Western Europe", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.weu" }, { "id": "geo.weu.t2001", "label": "Western Europe · 2001 to 2020", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.weu.t2001" }, { "id": "geo.weu.t2001.physical", "label": "Physical world and mathematics", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.weu.t2001.physical" }, { "id": "geo.weu.t2001.physical.earth", "label": "Earth sciences", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.weu.t2001.physical.earth" } ] } ], "excerpt": "GR4J, French for rural engineering model with 4 daily parameters, is a daily lumped rainfall-runoff model developed at Cemagref in France and widely used across climates.", "snippet": "GR4J, French for rural engineering model with 4 daily parameters, is a daily lumped rainfall-runoff model developed at Cemagref in France and widely used across climates.", "node": "physical.earth.water.hydrology.modeling-software", "markdown": "# GR4J\n\n**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)\n\n| Key fact | Detail |\n|---|---|\n| Type | Daily lumped conceptual rainfall-runoff model of the soil moisture accounting family[4](https://www.sciencedirect.com/science/article/pii/S0022169403002257) |\n| 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) |\n| 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) |\n| 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) |\n| 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) |\n| 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) |\n| 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) |\n\n## History\n\nDevelopment 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)\n\nThe 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)\n\n## Structure\n\nGR4J 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)\n\n**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)\n\n**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)\n\n**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)\n\nSeveral 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)\n\n## Calibration and software\n\nGiven 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)\n\nThe 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)\n\n## Variants and temporal resolutions\n\nThe 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)\n\n## Performance and applications\n\nThe 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)\n\nOperational 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)\n\n## Limitations\n\nGR4J'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)\n\n## References\n\n1. [GR4J – SRG, Source User Guide 4.11, eWater](https://ewater.atlassian.net/wiki/spaces/SD411/pages/53281097/GR4J%2B-%2BSRG)\n2. [GR4J technical reference, Kalix](https://kalix.org/docs/gr4j-technical/)\n3. [Tracking rainfall contribution in a conceptual runoff model, Journal of Hydrology](https://doi.org/10.1016/j.jhydrol.2026.135366)\n4. [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)\n5. [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)\n6. [Get Started with airGR, CRAN vignette](https://cran.r-project.org/web/packages/airGR/vignettes/V01_get_started.html)\n7. [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)\n8. [Bref historique, webGR, INRAE](https://webgr.inrae.fr/outils/modeles-hydrologiques/bref-historique)\n9. [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)\n10. [unifhycontrib.gr4 documentation](https://unifhy-org.github.io/unifhy/science/surfacelayer/unifhycontrib.gr4.SurfaceLayerComponent.html)\n11. [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)\n12. [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)\n13. [Evaluating Regional and Global Hydrological Models against Streamflow and Evapotranspiration Measurements, J. Hydrometeorology](https://doi.org/10.1175/JHM-D-15-0107.1)\n14. [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)\n15. [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)\n16. [DeepGR4J: A deep learning hybridization approach for conceptual rainfall-runoff modelling, Environmental Modelling & Software](https://doi.org/10.1016%2Fj.envsoft.2023.105831)\n17. [QDeepGR4J: Quantile-based ensemble of deep learning and GR4J hybrid rainfall-runoff models, arXiv](https://arxiv.org/abs/2510.05453)\n18. [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)\n19. [HESS preprint discussion of GR4J limitations](https://hess.copernicus.org/preprints/hess-2021-646/hess-2021-646-ATC2.pdf)\n20. [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)\n21. [A theoretical appraisal of the GR4J rainfall-runoff modelling framework, Journal of Hydrology](https://www.sciencedirect.com/science/article/pii/S0022169425017330)\n\n---\n*Topic: Encyclopedia › Physical world and mathematics › Earth sciences › Hydrology and ocean science › Hydrology › Hydrological modeling and software*\n\n*Initially written Oct 15, 2026*\n\n*Copyright 2026 EdgeChat AI, a subsidiary of Biostate AI.*\n\nLicense: Edgepedia Community License 1.0, https://www.edgepedia.org/edgepedia/license\n", "same_as": [], "url": "https://www.edgepedia.org/gr4j", "markdown_url": "https://www.edgepedia.org/gr4j.md", "license": { "name": "Edgepedia Community License 1.0", "url": "https://www.edgepedia.org/edgepedia/license", "summary": "Free with credit, commercial use included. AI training is open to everyone. For other uses, organizations over USD 100M in revenue or 100M monthly users license separately.", "spdx": "LicenseRef-Edgepedia-Community-1.0" }, "credit": "\"GR4J\", Edgepedia (EdgeChat), https://www.edgepedia.org/gr4j. Edgepedia Community License 1.0.", "credit_md": "\"[GR4J](https://www.edgepedia.org/gr4j)\", Edgepedia (EdgeChat), [https://www.edgepedia.org/gr4j](https://www.edgepedia.org/gr4j). [Edgepedia Community License 1.0](https://www.edgepedia.org/edgepedia/license).", "credit_html": "\"GR4J\", Edgepedia (EdgeChat), https://www.edgepedia.org/gr4j. Edgepedia Community License 1.0.", "speakable": "GR4J, French for rural engineering model with 4 daily parameters, is a daily lumped rainfall-runoff model developed at Cemagref in France and widely used across climates." }