Developing a reliable strategy to infer the effective soil hydraulic properties from field evaporation experiments for agro-hydrological models

The Richards equation has been widely used for simulating soil water movement. However, the take-up of agro-hydrological models using the basic theory of soil water flow for optimizing irrigation, fertilizer and pesticide practices is still low. This is partly due to the difficulties in obtaining accurate values for soil hydraulic properties at a field scale. Here, we use an inverse technique to deduce the effective soil hydraulic properties, based on measuring the changes in the distribution of soil water with depth in a fallow field over a long period, subject to natural rainfall and evaporation using a robust micro Genetic Algorithm. A new optimized function was constructed from the soil water contents at different depths, and the soil water at field capacity. The deduced soil water retention curve was approximately parallel but higher than that derived from published pedo-tranfer functions for a given soil pressure head. The water contents calculated from the deduced soil hydraulic properties were in good agreement with the measured values. The reliability of the deduced soil hydraulic properties was tested in reproducing data measured from an independent experiment on the same soil cropped with leek. The calculation of root water uptake took account for both soil water potential and root density distribution. Results show that the predictions of soil water contents at various depths agree fairly well with the measurements, indicating that the inverse analysis is an effective and reliable approach to estimate soil hydraulic properties, and thus permits the simulation of soil water dynamics in both cropped and fallow soils in the field accurately

LPMLE3 : a novel 1-D approach to study water flow in streambeds using heat as a tracer

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

A chromosome-level Amaranthus cruentus genome assembly highlights gene family evolution and biosynthetic gene clusters that may underpin the nutritional value of this traditional crop

Traditional crops historically provided accessible and affordable nutrition to millions of rural dwellers but have been neglected, with most modern agricultural systems over reliant on a small number of internationally-traded crops. Traditional crops are typically well-adapted to local agro-ecological conditions and many are nutrient-dense. They can play a vital role in local food systems through enhanced nutrition (especially where diets are dominated by starch crops), food security and livelihoods for smallholder farmers, and a climate-resilient and biodiverse agriculture. Using short-read, long-read and phased sequencing technologies we generated a high-quality chromosome-level genome assembly for Amaranthus cruentus, an under-researched crop with micronutrient- and protein-rich leaves and gluten-free seed, but lacking improved varieties, with respect to productivity and quality traits. The 370.9 MB genome demonstrates a shared whole genome duplication with a related species, Amaranthus hypochondriacus. Comparative genome analysis indicates chromosomal loss and fusion events following genome duplication that are common to both species, as well as fission of chromosome 2 in A. cruentus alone, giving rise to a haploid chromosome number of 17 (versus 16 in A. hypochondriacus). Genomic features potentially underlying the nutritional value of this crop include two A. cruentus-specific genes with a likely role in phytic acid synthesis (an anti-nutrient), expansion of ion transporter gene families, and identification of biosynthetic gene clusters conserved within the amaranth lineage. The A. cruentus genome assembly will underpin much-needed research and global breeding efforts to develop improved varieties for economically viable cultivation and realisation of the benefits to global nutrition security and agrobiodiversity

Estimation of lacustrine groundwater discharge using heat as a tracer and vertical hydraulic gradients – a comparison

Lacustrine groundwater discharge (LGD) can play a major role in water and nutrient balances of lakes. Unfortunately, studies often neglect this input path due to methodological difficulties in the determination. In a previous study we described a method which allows the estimation of LGD and groundwater recharge using hydraulic head data and groundwater net balances based on meteorological data. The aim of this study is to compare these results with discharge rates estimated by inverse modelling of heat transport using temperature profiles measured in lake bed sediments. We were able to show a correlation between the fluxes obtained with the different methods, although the time scales of the methods differ substantially. As a consequence, we conclude that the use of hydraulic head data and meteorologically-based groundwater net balances to estimate LGD is limited to time scales similar to the calibration period

Empirical quantification of lacustrine groundwater discharge – different methods and their limitations

Groundwater discharge into lakes (lacustrine groundwater discharge, LGD) can be an important driver of lake eutrophication. Its quantification is difficult for several reasons, and thus often neglected in water and nutrient budgets of lakes. In the present case several methods were applied to determine the expansion of the subsurface catchment, to reveal areas of main LGD and to identify the variability of LGD intensity. Size and shape of the subsurface catchment served as a prerequisite in order to calculate long-term groundwater recharge and thus the overall amount of LGD. Isotopic composition of near-shore groundwater was investigated to validate the quality of catchment delineation in near-shore areas. Heat as a natural tracer for groundwater–surface water interactions was used to find spatial variations of LGD intensity. Via an analytical solution of the heat transport equation, LGD rates were calculated from temperature profiles of the lake bed. The method has some uncertainties, as can be found from the results of two measurement campaigns in different years. The present study reveals that a combination of several different methods is required for a reliable identification and quantification of LGD and groundwater-borne nutrient loads

A comparative modeling study of a dual tracer experiment in a large lysimeter under atmospheric conditions.

In this paper, five model approaches with different physical and mathematical concepts varying in their model complexity and requirements were applied to identify the transport processes in the unsaturated zone. The applicability of these model approaches were compared and evaluated investigating two tracer breakthrough curves (bromide, deuterium) in a cropped, free-draining lysimeter experiment under natural atmospheric boundary conditions. The data set consisted of time series of water balance, depth resolved water contents, pressure heads and resident concentrations measured during 800 days. The tracer transport parameters were determined using a simple stochastic (stream tube model), three lumped parameter (constant water content model, multi-flow dispersion model, variable flow dispersion model) and a transient model approach. All of them were able to fit the tracer breakthrough curves. The identified transport parameters of each model approach were compared. Despite the differing physical and mathematical concepts the resulting parameters (mean water contents, mean water flux, dispersivities) of the five model approaches were all in the same range. The results indicate that the flow processes are also describable assuming steady state conditions. Homogeneous matrix flow is dominant and a small pore volume with enhanced flow velocities near saturation was identified with variable saturation flow and transport approach. The multi-flow dispersion model also identified preferential flow and additionally suggested a third less mobile flow component. Due to high fitting accuracy and parameter similarity all model approaches indicated reliable results

From submarine to lacustrine groundwater discharge

Submarine groundwater discharge (SGD) and its role in marine nutrient cycling are well known since the last decade. The freshwater equivalent, lacustrine groundwater discharge (LGD), is often still disregarded, although first reports of LGD are more than 50 years old. We identify nine different reasons why groundwater has long been disregarded in both freshwater and marine environments such as invisibility of groundwater discharge, the size of the interface and its difficult accessibility. Although there are some fundamental differences in the hydrology of SGD and LGD, caused primarily by seawater recirculation that occurs only in cases of SGD, there are also a lot of similarities such as a focusing of discharge to near-shore areas. Nutrient concentrations in groundwater near the groundwater–surface water interface might be anthropogenically enriched. Due to spatial heterogeneity of aquifer characteristics and biogeochemical processes, the quantification of groundwater-borne nutrient loads is challenging. Both nitrogen and phosphorus might be mobile in near-shore aquifers and in a lot of case studies large groundwater-borne nutrient loads have been reported