hydrological model
{{Short description|Predicting and managing water resources}}
A hydrologic model is a simplification of a real-world system (e.g., surface water, soil water, wetland, groundwater, estuary) that aids in understanding, predicting, and managing water resources. Both the flow and quality of water are commonly studied using hydrologic models.
Analog models
Prior to the advent of computer models, hydrologic modeling used analog models to simulate flow and transport systems. Unlike mathematical models that use equations to describe, predict, and manage hydrologic systems, analog models use non-mathematical approaches to simulate hydrology.
Two general categories of analog models are common; scale analogs that use miniaturized versions of the physical system and process analogs that use comparable physics (e.g., electricity, heat, diffusion) to mimic the system of interest.
= Scale analogs =
File:Miss-Basin-Model2-s.jpg (US Army Corps of Engineers, 2006)]]
Scale models offer a useful approximation of physical or chemical processes at a size that allows for greater ease of visualization.{{Cite journal|last=Rodhe|first=A.|date=2012-09-03|title=Physical models for classroom teaching in hydrology|journal=Hydrol. Earth Syst. Sci.|volume=16|issue=9|pages=3075–3082|doi=10.5194/hess-16-3075-2012|issn=1607-7938|bibcode=2012HESS...16.3075R|doi-access=free}} The model may be created in one (core, column), two (plan, profile), or three dimensions, and can be designed to represent a variety of specific initial and boundary conditions as needed to answer a question.
Scale models commonly use physical properties that are similar to their natural counterparts (e.g., gravity, temperature). Yet, maintaining some properties at their natural values can lead to erroneous predictions.{{Cite journal|last=Beven|first=Keith|title=Changing ideas in hydrology – The case of physically-based models|journal=Journal of Hydrology|language=en|volume=105|issue=1–2|pages=157–172|doi=10.1016/0022-1694(89)90101-7|bibcode=1989JHyd..105..157B|year=1989}} Properties such as viscosity, friction, and surface area must be adjusted to maintain appropriate flow and transport behavior. This usually involves matching dimensionless ratios (e.g., Reynolds number, Froude number).
File:Physical Aquifer Model.jpg
Groundwater flow can be visualized using a scale model built of acrylic and filled with sand, silt, and clay.Humphrey, M.D., 1992. Experimental design of physical aquifer models for evaluation of groundwater remediation strategies (Doctoral dissertation). Water and tracer dye may be pumped through this system to represent the flow of the simulated groundwater. Some physical aquifer models are between two and three dimensions, with simplified boundary conditions simulated using pumps and barriers.{{cite journal | last1 = Lee | first1 = S.S. | last2 = Kim | first2 = J.S. | last3 = Kim | first3 = D.J. | year = 2001 | title = Monitoring of drawdown pattern during pumping in an unconfined physical aquifer model | journal = Hydrological Processes | volume = 15 | issue = 3| pages = 479–492 | doi=10.1002/hyp.162| bibcode = 2001HyPr...15..479L | s2cid = 129919508 }}
= Process analogs =
Process analogs are used in hydrology to represent fluid flow using the similarity between Darcy's law, Ohm's law, Fourier's law, and Fick's law. The analogs to fluid flow are the flux of electricity, heat, and solutes, respectively.Principles of Soil and Plant Water Relations
https://books.google.com/books?isbn=0124200788 The corresponding analogs to fluid potential are voltage, temperature, and solute concentration (or chemical potential). The analogs to hydraulic conductivity are electrical conductivity, thermal conductivity, and the solute diffusion coefficient.
An early process analog model was an electrical network model of an aquifer composed of resistors in a grid.{{Cite web |url=http://www.isws.illinois.edu/hilites/achieve/images/gwmodded06.jpg |title=Archived copy |access-date=2017-05-01 |archive-date=2015-12-29 |archive-url=https://web.archive.org/web/20151229032549/http://www.isws.illinois.edu/hilites/achieve/images/gwmodded06.jpg |url-status=dead }} Voltages were assigned along the outer boundary, and then measured within the domain. Electrical conductivity paper{{Cite web | url=https://www.pasco.com/prodGroups/conductive-paper-and-pen/index.cfm | title=Conductive Paper and Pen: PASCO}} can also be used instead of resistors.
Statistical models
Statistical models are a type of mathematical model that are commonly used in hydrology to describe data, as well as relationships between data.Beard, Leo R. Statistical methods in hydrology. HYDROLOGIC ENGINEERING CENTER DAVIS CA, 1962. Using statistical methods, hydrologists develop empirical relationships between observed variables,{{Cite journal|last=Wallis|first=James R.|date=1965-12-01|title=Multivariate statistical methods in hydrology—A comparison using data of known functional relationship|journal=Water Resources Research|language=en|volume=1|issue=4|pages=447–461|doi=10.1029/WR001i004p00447|issn=1944-7973|bibcode=1965WRR.....1..447W}} find trends in historical data,{{Cite journal|last=Hamed|first=Khaled H.|date=2008-02-01|title=Trend detection in hydrologic data: The Mann–Kendall trend test under the scaling hypothesis|journal=Journal of Hydrology|volume=349|issue=3–4|pages=350–363|doi=10.1016/j.jhydrol.2007.11.009|bibcode=2008JHyd..349..350H}} or forecast probable storm or drought events.Yevjevich, Vujica. Probability and statistics in hydrology. Fort Collins, CO: Water resources publications, 1972.
= Moments =
Statistical moments (e.g., mean, standard deviation, skewness, kurtosis) are used to describe the information content of data. These moments can then be used to determine an appropriate frequency distribution,Zaharia, L. "L-MOMENTS AND THEIR USE IN MAXIMUM DISCHARGES'ANALYSIS IN CURVATURE CARPATHIANS REGION." Aerul si Apa. Componente ale Mediului (2013): 119. which can then be used as a probability model.{{Cite book|title=Computational Probability Applications|last1=Vargo|first1=Erik|last2=Pasupathy|first2=Raghu|last3=Leemis|first3=Lawrence M.|date=2017-01-01|publisher=Springer International Publishing|isbn=9783319433158|editor-last=Glen|editor-first=Andrew G.|series=International Series in Operations Research & Management Science|pages=149–164|language=en|doi=10.1007/978-3-319-43317-2_12|editor-last2=Leemis|editor-first2=Lawrence M.|citeseerx = 10.1.1.295.9820}} Two common techniques include L-moment ratios{{Cite journal|last1=PEEL|first1=MURRAY C.|last2=WANG|first2=Q. J.|last3=VOGEL|first3=RICHARD M.|last4=McMAHON|first4=THOMAS A.|title=The utility of L-moment ratio diagrams for selecting a regional probability distribution|journal=Hydrological Sciences Journal|volume=46|issue=1|pages=147–155|doi=10.1080/02626660109492806|year=2001|s2cid=14783093|doi-access=free|bibcode=2001HydSJ..46..147P }} and Moment-Ratio Diagrams.{{Cite journal|last1=Bobee|first1=B.|last2=Perreault|first2=L.|last3=Ashkar|first3=F.|date=1993-03-01|title=Two kinds of moment ratio diagrams and their applications in hydrology|journal=Stochastic Hydrology and Hydraulics|language=en|volume=7|issue=1|pages=41–65|doi=10.1007/BF01581566|issn=0931-1955|bibcode=1993SHH.....7...41B|s2cid=122128745}}
The frequency of extremal events, such as severe droughts and storms, often requires the use of distributions that focus on the tail of the distribution, rather than the data nearest the mean. These techniques, collectively known as extreme value analysis, provide a methodology for identifying the likelihood and uncertainty of extreme events.{{Cite journal|last=Sharma|first=T. C.|date=1998-03-30|title=An analysis of non-normal Markovian extremal droughts|journal=Hydrological Processes|language=en|volume=12|issue=4|pages=597–611|doi=10.1002/(sici)1099-1085(19980330)12:4<597::aid-hyp596>3.0.co;2-n|bibcode=1998HyPr...12..597S|issn=1099-1085}}{{Cite journal|last1=Katz|first1=Richard W|last2=Parlange|first2=Marc B|last3=Naveau|first3=Philippe|date=2002-08-01|title=Statistics of extremes in hydrology|journal=Advances in Water Resources|volume=25|issue=8–12|pages=1287–1304|doi=10.1016/S0309-1708(02)00056-8|bibcode=2002AdWR...25.1287K|url=http://infoscience.epfl.ch/record/53928}} Examples of extreme value distributions include the Gumbel, Pearson, and generalized extreme value. The standard method for determining peak discharge uses the log-Pearson Type III (log-gamma) distribution and observed annual flow peaks.https://water.usgs.gov/osw/bulletin17b/dl_flow.pdf {{Bare URL PDF|date=March 2022}}
= Correlation analysis =
The degree and nature of correlation may be quantified, by using a method such as the Pearson correlation coefficient, autocorrelation, or the T-test.Helsel, Dennis R., and Robert M. Hirsch. Statistical methods in water resources. Vol. 49. Elsevier, 1992 The degree of randomness or uncertainty in the model may also be estimated using stochastics,{{Cite journal|author1-link=Lynn Gelhar|last=Gelhar|first=Lynn W.|date=1986-08-01|title=Stochastic subsurface hydrology from theory to applications|journal=Water Resources Research|language=en|volume=22|issue=9S|pages=135S–145S|doi=10.1029/WR022i09Sp0135S|issn=1944-7973|bibcode=1986WRR....22R.135G}} or residual analysis.{{Cite journal|last1=Gupta|first1=Hoshin Vijai|last2=Sorooshian|first2=Soroosh|last3=Yapo|first3=Patrice Ogou|date=1998-04-01|title=Toward improved calibration of hydrologic models: Multiple and noncommensurable measures of information|journal=Water Resources Research|language=en|volume=34|issue=4|pages=751–763|doi=10.1029/97WR03495|issn=1944-7973|bibcode=1998WRR....34..751G|url=https://escholarship.org/uc/item/3zg7f6c9|doi-access=free}} These techniques may be used in the identification of flood dynamics,{{Cite journal|last1=Ouarda|first1=Taha B. M. J.|last2=Girard|first2=Claude|last3=Cavadias|first3=George S.|last4=Bobée|first4=Bernard|date=2001-12-10|title=Regional flood frequency estimation with canonical correlation analysis|journal=Journal of Hydrology|volume=254|issue=1–4|pages=157–173|doi=10.1016/S0022-1694(01)00488-7|bibcode=2001JHyd..254..157O}}{{Cite journal|last1=Ribeiro-Corréa|first1=J.|last2=Cavadias|first2=G.S.|last3=Clément|first3=B.|last4=Rousselle|first4=J.|title=Identification of hydrological neighborhoods using canonical correlation analysis|journal=Journal of Hydrology|language=en|volume=173|issue=1–4|pages=71–89|doi=10.1016/0022-1694(95)02719-6|bibcode=1995JHyd..173...71R|year=1995}} storm characterization,{{Cite journal|last=Marshall|first=R.J.|title=The estimation and distribution of storm movement and storm structure, using a correlation analysis technique and rain-gauge data|journal=Journal of Hydrology|language=en|volume=48|issue=1–2|pages=19–39|doi=10.1016/0022-1694(80)90063-3|bibcode=1980JHyd...48...19M|year=1980}}{{Cite journal|last1=Nathan|first1=R. J.|last2=McMahon|first2=T. A.|date=1990-07-01|title=Evaluation of automated techniques for base flow and recession analyses|journal=Water Resources Research|language=en|volume=26|issue=7|pages=1465–1473|doi=10.1029/WR026i007p01465|issn=1944-7973|bibcode=1990WRR....26.1465N}} and groundwater flow in karst systems.{{Cite journal|title=Contribution of correlation and spectral analyses to the regional study of a large karst aquifer (Charente, France) |language=en|doi=10.1016/S0022-1694(97)00155-8|bibcode=1998JHyd..205..217L|volume=205|issue=3–4 |journal=Journal of Hydrology|pages=217–231 | last1 = Larocque | first1 = M.|year=1998 }}
Regression analysis is used in hydrology to determine whether a relationship may exist between independent and dependent variables. Bivariate diagrams are the most commonly used statistical regression model in the physical sciences, but there are a variety of models available from simplistic to complex.{{Cite journal|last=Geological Survey (U.S.)|date=1950-01-01|title=Geological Survey water-supply paper.|journal=Geological Survey Water-supply Paper.|language=en|issn=0083-1131|oclc=1422999}} In a bivariate diagram, a linear or higher-order model may be fitted to the data.
Factor analysis and principal component analysis are multivariate statistical procedures used to identify relationships between hydrologic variables.{{Cite journal|last1=Matalas|first1=N. C.|last2=Reiher|first2=Barbara J.|date=1967-03-01|title=Some comments on the use of factor analyses|journal=Water Resources Research|language=en|volume=3|issue=1|pages=213–223|doi=10.1029/WR003i001p00213|issn=1944-7973|bibcode=1967WRR.....3..213M}}{{cite journal | last1 = Pearson | first1 = K | year = 1901 | title = LIII. On lines and planes of closest fit to systems of points in space | url = https://zenodo.org/record/1430636| journal = The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science | volume = 2 | issue = 11| pages = 559–572 | doi=10.1080/14786440109462720}}
Convolution is a mathematical operation on two different functions to produce a third function. With respect to hydrologic modeling, convolution can be used to analyze stream discharge's relationship to precipitation. Convolution is used to predict discharge downstream after a precipitation event. This type of model would be considered a “lag convolution”, because of the predicting of the “lag time” as water moves through the watershed using this method of modeling.
Time-series analysis is used to characterize temporal correlation within a data series as well as between different time series. Many hydrologic phenomena are studied within the context of historical probability. Within a temporal dataset, event frequencies, trends, and comparisons may be made by using the statistical techniques of time series analysis.Salas, Jose D. Applied modeling of hydrologic time series. Water Resources Publication, 1980. The questions that are answered through these techniques are often important for municipal planning, civil engineering, and risk assessments.
Markov chains are a mathematical technique for determine the probability of a state or event based on a previous state or event.{{Cite web|url=http://setosa.io/ev/markov-chains/|title=Markov Chains explained visually|website=Explained Visually|access-date=2017-04-21}} The event must be dependent, such as rainy weather. Markov Chains were first used to model rainfall event length in days in 1976,{{Cite journal|last1=Haan|first1=C. T.|last2=Allen|first2=D. M.|last3=Street|first3=J. O.|date=1976-06-01|title=A Markov Chain Model of daily rainfall|journal=Water Resources Research|language=en|volume=12|issue=3|pages=443–449|doi=10.1029/WR012i003p00443|issn=1944-7973|bibcode=1976WRR....12..443H}} and continues to be used for flood risk assessment and dam management.
Data-driven models
Data-driven models in hydrology emerged as an alternative approach to traditional statistical models, offering a more flexible and adaptable methodology for analysing and predicting various aspects of hydrological processes. While statistical models rely on rigorous assumptions about probability distributions, data-driven models leverage techniques from artificial intelligence, machine learning, and statistical analysis, including correlation analysis, time series analysis, and statistical moments, to learn complex patterns and dependencies from historical data. This allows them to make more accurate predictions and provide insights into the underlying processes.Abrahart, R. P., See, L. M., & Solomatine, D. P. (2008). Practical Hydroinformatics: Computational Intelligence and Technological Developments in Water Applications. Springer.
Since their inception in the latter half of the 20th century, data-driven models have gained popularity in the water domain, as they help improve forecasting, decision-making, and management of water resources. A couple of notable publications that use data-driven models in hydrology include "Application of machine learning techniques for rainfall-runoff modelling" by Solomatine and Siek (2004),Solomatine, D. P., & Siek, M. B. (2004). Application of machine learning techniques for rainfall-runoff modeling. Hydroinformatics: A wide range of technologies, 333-342. and "Data-driven modelling approaches for hydrological forecasting and prediction" by Valipour et al. (2021).Valipour, M., Shahsavani, D., & Choubin, B. (2021). Data-driven modeling approaches for hydrological forecasting and prediction. Journal of Hydrology, 597, 125996. These models are commonly used for predicting rainfall, runoff, groundwater levels, and water quality, and have proven to be valuable tools for optimizing water resource management strategies.
Conceptual models
Conceptual models represent hydrologic systems using physical concepts. The conceptual model is used as the starting point for defining the important model components. The relationships between model components are then specified using algebraic equations, ordinary or partial differential equations, or integral equations. The model is then solved using analytical or numerical procedures.
Conceptual models are commonly used to represent the important components (e.g., features, events, and processes) that relate hydrologic inputs to outputs.{{Cite journal|last1=Enemark|first1=Trine|last2=Peeters|first2=Luk J.M.|last3=Mallants|first3=Dirk|last4=Batelaan|first4=Okke|date=February 2019|title=Hydrogeological conceptual model building and testing: A review|url=https://linkinghub.elsevier.com/retrieve/pii/S0022169418309387|journal=Journal of Hydrology|language=en|volume=569|pages=310–329|doi=10.1016/j.jhydrol.2018.12.007|bibcode=2019JHyd..569..310E|hdl=2328/38835|s2cid=135032550 |hdl-access=free}} These components describe the important functions of the system of interest, and are often constructed using entities (stores of water) and relationships between these entitites (flows or fluxes between stores). The conceptual model is coupled with scenarios to describe specific events (either input or outcome scenarios).
For example, a watershed model could be represented using tributaries as boxes with arrows pointing toward a box that represents the main river. The conceptual model would then specify the important watershed features (e.g., land use, land cover, soils, subsoils, geology, wetlands, lakes), atmospheric exchanges (e.g., precipitation, evapotranspiration), human uses (e.g., agricultural, municipal, industrial, navigation, thermo- and hydro-electric power generation), flow processes (e.g., overland, interflow, baseflow, channel flow), transport processes (e.g.,
sediments, nutrients, pathogens), and events (e.g., low-, flood-, and mean-flow conditions).
Model scope and complexity is dependent on modeling objectives, with greater detail required if human or environmental systems are subject to greater risk. Systems modeling can be used for building conceptual models that are then populated using mathematical relationships.
Example 1
The linear-reservoir model (or Nash model) is widely used for rainfall-runoff analysis. The model uses a cascade of linear reservoirs along with a constant first-order storage coefficient, K, to predict the outflow from each reservoir (which is then used as the input to the next in the series).
The model combines continuity and storage-discharge equations, which yields an ordinary differential equation that describes outflow from each reservoir. The continuity equation for tank models is:
:
which indicates that the change in storage over time is the difference between inflows and outflows. The storage-discharge relationship is:
:
where K is a constant that indicates how quickly the reservoir drains; a smaller value indicates more rapid outflow. Combining these two equation yields
:
and has the solution:
:
Example 2
Instead of using a series of linear reservoirs, also the model of a non-linear reservoir can be used.[https://www.waterlog.info/pdf/reservoir.pdf Rainfall-runoff modelling using a non-linear reservoir]
In such a model the constant K in the above equation, that may also be called reaction factor, needs to be replaced by another symbol, say α (Alpha), to indicate the dependence of this factor on storage (S) and discharge (q).
In the left figure the relation is quadratic:
:α = 0.0123 q2 + 0.138 q - 0.112
= Governing equations =
Governing equations are used to mathematically define the behavior of the system. Algebraic equations are likely often used for simple systems, while ordinary and partial differential equations are often used for problems that change in space in time. Examples of governing equations include:
Manning's equation is an algebraic equation that predicts stream velocity as a function of channel roughness, the hydraulic radius, and the channel slope:
:
Darcy's law describes steady, one-dimensional groundwater flow using the hydraulic conductivity and the hydraulic gradient:
:
Groundwater flow equation describes time-varying, multidimensional groundwater flow using the aquifer transmissivity and storativity:
:
Advection-Dispersion equation describes solute movement in steady, one-dimensional flow using the solute dispersion coefficient and the groundwater velocity:
:
Poiseuille's law describes laminar, steady, one-dimensional fluid flow using the shear stress:
:
Cauchy's integral is an integral method for solving boundary value problems:
:
= Solution algorithms =
== Analytic methods ==
== Numeric methods ==
Many real-world mathematical models are too complex to meet the simplifying assumptions required for an analytic solution. In these cases, the modeler develops a numerical solution that approximates the exact solution. Solution techniques include the finite-difference and finite-element methods, among many others.
Specialized software may also be used to solve sets of equations using a graphical user interface and complex code, such that the solutions are obtained relatively rapidly and the program may be operated by a layperson or an end user without a deep knowledge of the system. There are model software packages for hundreds of hydrologic purposes, such as surface water flow, nutrient transport and fate, and groundwater flow.
Commonly used numerical models include SWAT, MODFLOW, FEFLOW, [https://www.acricfd.com/software/porflow/default.htm PORFLOW], MIKE SHE, and WEAP.
= Model calibration and evaluation =
File:Rainfall_runoff_relation_simulated_by_non-linear_resrvoir_model.png
Physical models use parameters to characterize the unique aspects of the system being studied. These parameters can be obtained using laboratory and field studies, or estimated by finding the best correspondence between observed and modelled behavior.{{cite journal |last1=Musyoka |first1=F.K |last2=Strauss |first2=P |last3=Zhao |first3=G |last4=Srinivasan |first4=R |last5=Klik |first5=A |title=Multi-Step Calibration Approach for SWAT Model Using Soil Moisture and Crop Yields in a Small Agricultural Catchment |journal=Water |date=2021 |doi=10.3390/w13162238 |language=en|volume=16|issue=13|page=2238 |issn=2073-4441 |doi-access=free|hdl=10261/253007 |hdl-access=free }}{{Cite journal|date=2014-07-01|title=A simple, regionally parameterized model for predicting nonpoint source areas in the northeastern US|journal=Journal of Hydrology: Regional Studies|language=en|volume=1|pages=74–91|doi=10.1016/j.ejrh.2014.06.003|issn=2214-5818|doi-access=free|last1=Archibald|first1=J.A.|last2=Buchanan|first2=B.P.|last3=Fuka|first3=D.R.|last4=Georgakakos|first4=C.B.|last5=Lyon|first5=S.W.|last6=Walter|first6=M.T.|bibcode=2014JHyRS...1...74A }}{{Cite journal|last1=Archibald|first1=Josephine A.|last2=Walter|first2=M. Todd|date=2014|title=Do Energy-Based PET Models Require More Input Data than Temperature-Based Models? — An Evaluation at Four Humid FluxNet Sites|url=https://onlinelibrary.wiley.com/doi/abs/10.1111/jawr.12137|journal=JAWRA Journal of the American Water Resources Association|language=en|volume=50|issue=2|pages=497–508|doi=10.1111/jawr.12137|bibcode=2014JAWRA..50..497A|s2cid=129781883 |issn=1752-1688}}{{Cite journal|last1=Knighton|first1=James|last2=Singh|first2=Kanishka|last3=Evaristo|first3=Jaivime|date=2020|title=Understanding Catchment-Scale Forest Root Water Uptake Strategies Across the Continental United States Through Inverse Ecohydrological Modeling|journal=Geophysical Research Letters|language=en|volume=47|issue=1|pages=e2019GL085937|doi=10.1029/2019GL085937|bibcode=2020GeoRL..4785937K|s2cid=213914582|issn=1944-8007|doi-access=free}} Between neighbouring catchments which have physical and hydrological similarities, the model parameters varies smoothly suggesting the spatial transferability of parameters.{{Cite journal|last1=Nepal|first1=Santosh|last2=Flügel|first2=Wolfgang-Albert|last3=Krause|first3=Peter|last4=Fink|first4=Manfred|last5=Fischer|first5=Christian|date=2017-07-30|title=Assessment of spatial transferability of process-based hydrological model parameters in two neighbouring catchments in the Himalayan Region|journal=Hydrological Processes|language=en|volume=31|issue=16|pages=2812–2826|doi=10.1002/hyp.11199|issn=1099-1085|bibcode=2017HyPr...31.2812N|doi-access=free}}
Model evaluation is used to determine the ability of the calibrated model to meet the needs of the modeler. A commonly used measure of hydrologic model fit is the Nash-Sutcliffe efficiency coefficient.
See also
References
External links
- http://drought.unl.edu/MonitoringTools/DownloadableSPIProgram.aspx
{{Agricultural water management}}
{{Computer modeling}}