Premium
LPMLE3: A novel 1‐D approach to study water flow in streambeds using heat as a tracer
Author(s) -
Schneidewind U.,
van Berkel M.,
Anibas C.,
Vandersteen G.,
Schmidt C.,
Joris I.,
Seuntjens P.,
Batelaan O.,
Zwart H. J.
Publication year - 2016
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1002/2015wr017453
Subject(s) - estimator , flow (mathematics) , thermal diffusivity , boundary (topology) , series (stratigraphy) , boundary value problem , mathematics , environmental science , computer science , hydrology (agriculture) , mathematical optimization , algorithm , geology , mathematical analysis , statistics , geometry , thermodynamics , geotechnical engineering , physics , paleontology
Abstract We introduce LPMLE3, a new 1‐D approach to quantify vertical water flow components at streambeds using temperature data collected in different depths. LPMLE3 solves the partial differential equation for coupled water flow and heat transport in the frequency domain. Unlike other 1‐D approaches it does not assume a semi‐infinite halfspace with the location of the lower boundary condition approaching infinity. Instead, it uses local upper and lower boundary conditions. As such, the streambed can be divided into finite subdomains bound at the top and bottom by a temperature‐time series. Information from a third temperature sensor within each subdomain is then used for parameter estimation. LPMLE3 applies a low order local polynomial to separate periodic and transient parts (including the noise contributions) of a temperature‐time series and calculates the frequency response of each subdomain to a known temperature input at the streambed top. A maximum‐likelihood estimator is used to estimate the vertical component of water flow, thermal diffusivity, and their uncertainties for each streambed subdomain and provides information regarding model quality. We tested the method on synthetic temperature data generated with the numerical model STRIVE and demonstrate how the vertical flow component can be quantified for field data collected in a Belgian stream. We show that by using the results in additional analyses, nonvertical flow components could be identified and by making certain assumptions they could be quantified for each subdomain. LPMLE3 performed well on both simulated and field data and can be considered a valuable addition to the existing 1‐D methods.