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.

Knowledge on groundwater–surface water interaction and especially on exchange fluxes between streams and aquifers is an important prerequisite for the study of transport and fate of contaminants and nutrients in the hyporheic zone. One possibility to quantify groundwater–surface water exchange fluxes is by using heat as an environmlental tracer. Modern field equipment including multilevel temperature sticks and the novel open‐source analysis tool LPML make this technique ever more attractive. The recently developed LPML method solves the one‐dimensional fluid flow and heat transport equation by combining a local polynomial method with a maximum likelihood estimator. In this study, we apply the LPML method on field data to quantify the spatial and temporal variability of vertical fluxes and their uncertainties from temperature–time series measured in a Belgian lowland stream. Over several months, temperature data were collected with multilevel temperature sticks at the streambed top and at six depths for a small stream section. Long‐term estimates show a range from gaining fluxes of −291 mm day−1 to loosing fluxes of 12 mm day−1; average seasonal fluxes ranged from −138 mm day−1 in winter to −16 mm day−1 in summer. With our analyses, we could determine a high spatial and temporal variability of vertical exchange fluxes for the investigated stream section. Such spatial and temporal variability should be taken into account in biogeochemical cycling of carbon, nutrients and metals and in fate analysis of contaminant plumes. In general, the stream section was gaining during most of the observation period. Two short‐term high stream stage events, seemingly caused by blockage of the stream outlet, led to a change in flow direction from gaining to losing conditions. We also found more discharge occurring at the outer stream bank than at the inner one indicating a local flow‐through system. With the conducted analyses, we were able to advance our understanding of the regional groundwater flow system. Copyright © 2015 John Wiley & Sons, Ltd.

Complex interactions occur in riparian wetlands between groundwater, surface water and climatic conditions. Knowledge of the hydrology of these systems is necessary to understand their functioning and their value and models are a useful and probably essential tool to capture their hydrological complexity. In this study, a 2D-model describing saturated-unsaturated water flow is applied to a transect through a groundwaterfed riparian wetland located along the middle reach of the river Dijle. The transect has high levees close to the river and a depression further into the floodplain. Scaling factors are introduced to describe the variability of soil hydraulic properties along the transect. Preliminary model calculations for one year show a good agreement between model calculations and measurements and demonstrate the capability of the model to capture the internal groundwater dynamics. Seasonal variations in soil moisture are reproduced well by the model thus translating external hydrological boundary conditions to root zone conditions. The model proves to be a promising tool for assessing effects of changes in hydrological boundary conditions on vegetation type distribution and to gain more insight in the highly variable internal flow processes of riparian wetlands.

In places where tidal marshes were formerly embanked for agricultural land use, these marshes are nowadays increasingly restored with the aim to regain important ecosystem services. However, there is growing evidence that restored tidal marshes and their services develop slowly and differ from natural tidal marshes in many aspects. Here we focus on groundwater dynamics, because these affect several key ecosystem functions and services, such as nutrient cycling and vegetation development. We hypothesize that groundwater dynamics in restored tidal marshes are reduced as compared to natural marshes because of the difference in soil structure. In the macro-tidal Schelde estuary (Belgium), in both a natural and a restored (since 2006) freshwater tidal marsh, we measured depth profiles of soil properties (grain size distribution, LOI (loss on ignition), moisture content and bulk density) and temporal dynamics of groundwater levels along a transect with increasing distance from a tidal creek. X-ray micro CT-scanning was used to quantify soil macroporosity. The restored marsh has a twolayered soil stratigraphy with a topsoil of freshly accreted sediment (ranging in depth between 10 and 60 cm, deposited since 2006) and a subsoil of compact relict agricultural soil. We found that both the soil in the natural marsh and the topsoil of the restored marsh consist of loosely packed sediment rich in macropores and organic matter, whereas the relict agricultural soil in the restored marsh is densely packed and has few macropores. Our results show that groundwater level fluctuations in the restored marsh are restricted to the top layer of newly deposited sediment (i.e. on average 0.08 m depth). Groundwater level fluctuations in the natural marsh occur over a larger depth of the soil profile (i.e. on average 0.28 m depth). As a consequence, the reduced groundwater dynamics in restored tidal marshes are expected to alter the subsurface fluxes of water and nutrients, the sourcesink function and the development of marsh vegetation.