Data assimilation of in situ soil moisture measurements in hydrological models: third annual doctoral progress report, work plan and achievements

Efficient water utilization and optimal water supply/distribution to increase food and fodder productivity are of utmost importance in confronting worldwide water scarcity, climate change, growing populations and increasing water demands. In this respect, irrigation efficiency, which is influenced by the type of irrigation and irrigation scheduling, is an essential issue for achieving higher productivity. To improve irrigation strategies in precision agriculture, soil water status can be more accurately described using a combination of advanced monitoring and modeling. Our study focuses on the combination of high resolution hydrological data with hydrological models that predict water flow and solute (pollutants and salts) transport and water redistribution in agricultural soils under irrigation. Field plots of a potato farmer in a sandy region in Belgium were instrumented to continuously monitor soil moisture and water potential before, during and after irrigation in dry summer periods. The aim is to optimize the irrigation process by assimilating online sensor field data into process based models. This research is part of Activity 305 ‘Precision agriculture and remote sensing’ of the VITO GWO and is also part of the strategic cooperation with UGent within the platform ‘Managing Natural Resources’. Over the past 2 years, we applied a combination of in-situ monitoring and numerical modeling -Hydrus 1D- to estimate water content fluctuations in a heterogeneous sandy grassland soil under irrigation with water table fluctuating between 80 and 155 cm. Over the last year, more sampling and analyses were carried out to further characterize the hydraulic properties over the entire field. Modeling results for the field demonstrated clearly the profound effect of the position of the GWL, and to a lesser extent, the effect of spatially variable soil hydraulic properties (Ks, n and α) on the estimated water content in the sandy two-layered soil under grass. Our results show that currently applied uniform water distribution using sprinkler irrigation seems not to be efficient since at locations with shallow groundwater, the amount of water applied will be excessive as compared to the plant requirements while in locations with a deeper GWL, requirements will not be met. To derive the optimal parameter set best describing the measured soil moisture content, 37 optimization scenarios were conducted with two to six parameters using various parameter combinations for the two soil layers. The best performing parameter optimization scenario was a 2-parameter scenario with Ks optimized for each layer. The results showed a better identifiability of the parameters (less correlations among parameters) with equal performance as compared to three, four or six parameter optimization. Model predictions using the calibrated model (with data from 2012) for an independent data set of soil moisture data in the validation period (2013) showed satisfactory performance of the model in view of irrigation management purposes. Comparing the degree of water stress for different optimization scenarios of groundwater depth, showed that grass was exposed to water stress in summer in 2013 but not for such a long period as compared to the 2012 growing season. The degree of water stress simulated with Hydrus 1D suggested to increase the irrigation amount in 2012 and 2013 and at least one or two times in the summer (June and July) and further distributing the amount of irrigation during the growing season, instead of using a huge amount of irrigation later in the season, as is common practice by the farmer. A second part of the study focused on finding a relation between measured soil hydraulic properties and apparent electrical conductivity ECa. Our measurements of hydraulic properties of the field clearly confirm that there is considerable spatial variability in the field and that this has an impact on the simulation of soil moisture content. Therefore this should be taken into account when upscaling soil hydraulic properties to the field scale in order to in understand and model flow, solute and energy fluxes in the field and develop strategies for efficient irrigation. Upscaling soil hydraulic properties to the field scale can be done by linking them to apparent electrical conductivity (ECa), which can be measured efficiently and inexpensively so a spatially dense dataset for describing within-field spatial soil variability can be generated. In this study relations between the spatial variation of soil hydraulic properties and apparent soil electrical conductivity ECa measured with EM38 and DUALEM-21S sensors at two depths of explorations (DOE) 0-50 and 0-100 cm were investigated. Two predictive modelling approaches, i.e. i) a simple regression and ii) applying Archie’s laws for saturated and unsaturated conditions in combination with MVG equations, were developed and it was compared how they were able to explain the observed values of hydraulic parameters. Results demonstrated the spatial variability and heterogeneity of ECa and soil hydraulic properties Ks, α and n. We derived a regression relationship between log Ks and ECa measured with DUALEM (r2≥0.70) and with EM38 (r2>0.46) sensors. The predicted results were tested vs measured data and confirmed that the performance of DUALEMp,100-Ks model is relatively better than that of the same sensor with lower DOE and of the EM38 sensor (RMSE = 1.31 cmh-1, R2 = 0.55). The relationships between MVG shape parameters and ECa datasets were generally poor (0.05<R2<0.26). In the second approach, we showed that the water retention curve can be translated to ECa-(h) and ECa-Se relations by combining the MVG equations and Archie’s law. Results also show that reformulating the MVG equations based on ECa-Se relationships can help to estimate unsaturated hydraulic conductivity at the field scale. In the third year, a second study site has been set up in a nearby field where potatoes are grown and has been instrumented with soil moisture sensors, tensiometers, groundwater level loggers and a weather station. Field hydraulic properties for the field will be derived using the equations developed for the first study site and the modeling approach developed for the first field will be tested here. Also quasi 3D-modelling of water flow at the field scale will be conducted. In this modeling set-up, the field will be modeled as a collection of 1D-columns representing the different field conditions (combination of soil properties, GWL, root zone depth). Combining this model with crop based models such as LINGRA-N or Aquacrop gives a direct simulation of the impact of irrigation strategies on crop yield at the field scale

Data assimilation of in situ soil moisture measurements in hydrological models: first annual doctoral progress report, work plan and achievements

Water scarcity and the presence of water of good quality is a serious public concern since it determines the availability of water to society. Water scarcity especially in arid climates and due to extreme droughts related to climate change drive water use technologies such as irrigation to become more efficient and sustainable. Plant root water and nutrient uptake is one of the most important processes in subsurface unsaturated flow and transport modeling, as root uptake controls actual plant evapotranspiration, water recharge and nutrient leaching to the groundwater, and exerts a major influence on predictions of global climate models. To improve irrigation strategies, water flow needs to be accurately described using advanced monitoring and modeling. Our study focuses on the assimilation of hydrological data in hydrological models that predict water flow and solute (pollutants and salts) transport and water redistribution in agricultural soils under irrigation. Field plots of a potato farmer in a sandy region in Belgium were instrumented to continuously monitor soil moisture and water potential before, during and after irrigation in dry summer periods. The aim is to optimize the irrigation process by assimilating online sensor field data into process based models. Over the past year, we demonstrated the calibration and optimization of the Hydrus 1D model for an irrigated grassland on sandy soil. Direct and inverse calibration and optimization for both heterogeneous and homogeneous conceptualizations was applied. Results show that Hydrus 1D closely simulated soil water content at five depths as compared to water content measurements from soil moisture probes, by stepwise calibration and local sensivity analysis and optimization the Ks, n and α value in the calibration and optimization analysis. The errors of the model, expressed by deviations between observed and modeled soil water content were, however, different for each individual depth. The smallest differences between the observed value and soil-water content were attained when using an automated inverse optimization method. The choice of the initial parameter value can be optimized using a stepwise approach. Our results show that statistical evaluation coefficients (R2, Ce and RMSE) are suitable benchmarks to evaluate the performance of the model in reproducing the data. The degree of water stress simulated with Hydrus 1D suggested to increase irrigation at least one time, i.e. at the beginning of the simulation period and further distribute the amount of irrigation during the growing season, instead of using a huge amount of irrigation later in the season. In the next year, we will further look for to the best method (using soft data and methods for instance PTFs, EMI, Penetrometer) to derive and predict the spatial variability of soil hydraulic properties (saturated hydraulic conductivity) of the soil and link to crop yield at the field scale. Linear and non-linear pedotransfer functions (PTFs) have been assessed to predict penetrometer resistance of soils from their water status (matric potential, ψ and degree of saturation, S) and bulk density, ρb, and some other soil properties such as sand content, Ks etc. The geophysical EMI (electromagnetic induction) technique provides a versatile and robust field instrument for determining apparent soil electrical conductivity (ECa). ECa, a quick and reliable measurement, is one of ancillary properties (secondary information) of soil, can improve the spatial and temporal estimation of soil characteristics e.g., salinity, water content, texture, prosity and bulk density at different scales and depths. According to previous literature on penetrometer measurements, we determined the effective stress and used some models to find the relationships between soil properties, especially Ks, and penetrometer resistance as one of the prediction methods for Ks. The initial results obtained in the first yearshowed that a new data set would be necessary to validate the results of this part. In the third year, quasi 3D-modelling of water flow at the field scale will be conducted. In this modeling set -up, the field will be modeled as a collection of 1D-columns representing the different field conditions (combination of soil properties, groundwater depth, root zone depth). The measured soil properties are extrapolated over the entire field by linking them to the available spatially distributed data (such as the EMI-images). The data set of predicted Ks and other soil properties for the whole field constructed in the previous steps will be used for parameterising the model. Sensitivity analysis ‘SA’ is essential to the model optimization or parametrization process. To avoid overparameterization, the use of global sensitivity analysis (SA) will be investigated. In order to include multiple objectives (irrigation management parameters, costs, …) in the parameter optimization strategy, multi-objective techniques such as AMALGAM have been introduced. We will investigate multi-objective strategies in the irrigation optimization

Sensitivity of water stress in a two-layered sandy grassland soil to variations in groundwater depth and soil hydraulic parameters

Monitoring and modelling tools may improve irrigation strategies in precision agriculture. We used non-invasive soil moisture monitoring, a crop growth and a soil hydrological model to predict soil water content fluctuations and crop yield in a heterogeneous sandy grassland soil under supplementary irrigation. The sensitivity of the soil hydrological model to hydraulic parameters, water stress, crop yield and lower boundary conditions was assessed after integrating models. Free drainage and incremental constant head conditions were implemented in a lower boundary sensitivity analysis. A time-dependent sensitivity analysis of the hydraulic parameters showed that changes in soil water content are mainly affected by the soil saturated hydraulic conductivity K-s and the Mualem-van Genuchten retention curve shape parameters n and alpha. Results further showed that different parameter optimization strategies (two-, three-, four- or six-parameter optimizations) did not affect the calculated water stress and water content as significantly as does the bottom boundary. In this case, a two-parameter scenario, where K-s was optimized for each layer under the condition of a constant groundwater depth at 135-140 cm, performed best. A larger yield reduction, and a larger number and longer duration of stress conditions occurred in the free drainage condition as compared to constant boundary conditions. Numerical results showed that optimal irrigation scheduling using the aforementioned water stress calculations can save up to 12-22 % irrigation water as compared to the current irrigation regime. This resulted in a yield increase of 4.5-6.5 %, simulated by the crop growth model

Development of soft computing and applications in agricultural and biological engineering

Soft computing is a set of “inexact” computing techniques, which are able to model and analyze very complex problems. For these complex problems, more conventional methods have not been able to produce cost-effective, analytical, or complete solutions. Soft computing has been extensively studied and applied in the last three decades for scientific research and engineering computing. In agricultural and biological engineering, researchers and engineers have developed methods of fuzzy logic, artificial neural networks, genetic algorithms, decision trees, and support vector machines to study soil and water regimes related to crop growth, analyze the operation of food processing, and support decision-making in precision farming. This paper reviews the development of soft computing techniques. With the concepts and methods, applications of soft computing in the field of agricultural and biological engineering are presented, especially in the soil and water context for crop management and decision support in precision agriculture. The future of development and application of soft computing in agricultural and biological engineering is discussed

A site-specific and dynamic modeling system for zoning and optimizing variable rate irrigation in cotton

Cotton irrigation has been rapidly expanding in west Tennessee during the past decade. Variable rate irrigation is expected to enhance water use efficiency and crop yield in this region due to the significant field-scale soil spatial heterogeneity. A detailed understanding of the soil available water content within the effective root zone is needed to optimally schedule irrigation. In addition, site-specific crop-yield mathematical relationships should be established to identify optimum irrigation management. This study aimed to design and evaluate a site-specific modeling system for zoning and optimizing variable rate irrigation in cotton. The specific objectives of this study were to investigate (i) the spatial variability of soil attributes at the field-scale, (ii) site-specific cotton lint yieldwater relationships across all soil types, and (iii) multiple zoning strategies for variable rate irrigation scenarios. The field (73 ha) was sampled and apparent soil electrical conductivity (ECa) was measured. Landsat 8 satellite data was acquired, processed, and transformed to compare indicators of vegetation and soil response to cotton lint yields, variable irrigation rates, and the spatial variability of soil attributes. Multiple modeling scenarios were developed and examined. Although experiments were performed during two wet years, supplemental irrigation enhanced cotton yield across all soil types in comparison with rain-fed conditions. However, length of cropping season and rainfall distribution remarkably affected cotton response to supplemental irrigation. Geostatistical analysis showed spatial variability in soil textural components and water content was significant and correlated to yield patterns. There was as high as four-fold difference between available water content between coarse-textured and fine-textured soils on the study site. A good agreement was observed (RMSE = 0.052 cm3 cm-3 [cubic centimeter per cubic centimeter] and r = 0.88) between predicted and observed water contents. ECa and space images were useful proximal data to investigate soil spatial variability. The site-specific water production functions performed well at predicting cotton lint yield with RMSE equal to 0.131 Mg ha-1 [megagram per hectare] and 0.194 Mg ha-1 in 2013 and 2014, respectively. The findings revealed that variable rate irrigation with pie shape zones could enhance cotton lint yield under supplemental irrigation in west Tennessee

Engineering applications of soil information system (sis): precision irrigation and drainage systems design

The objectives of this internship were to demonstrate and apply the knowledge and technical training obtained during Doctor of Engineering coursework and to become familiar with the organizational approach to problems. These objectives were fulfilled in three commercial research and development projects in the field of precision agriculture. The first project involved optimization of a center pivot irrigation system in coordination with the system’s manufacturing company in order to apply irrigation water to maintain uniform soil water content across the field. An optimization-simulation model was developed for this purpose using a dynamic programming approach. The simulations run by the optimization model showed that the existing pivot speed prescription resulted in a more uniform soil water content across the field reducing the crop yield losses. The objective of the second internship project was to analyze the agricultural drainage industry for identification of the potential applications of the spatial soil information into agricultural drainage design and construction. In order to accomplish this task, a comprehensive literature survey was conducted with an emphasis on the drainage approximate equation and numerical methods. Software tools that are currently employed in drainage design and construction were evaluated. A detailed market analysis was conducted with a focus on the industry stakeholders. A strategic strengths, weaknesses, opportunities and threats (SWOT) analysis was conducted for the agricultural drainage industry using Porter’s five forces method. The last internship project involved investigation of the potential for using soil information obtained by SIS in the assessment of soil salinity. A correlation analysis was conducted between the soil paste extract electrical conductivity values measured in the laboratory on collected soil samples and those estimated using the soil resistivity values collected by the SIS, which is a measure of soil salinity. The results showed no clear correlations. While the internship projects provided the intern the opportunity to apply some of the analytical methods learned as part of the Doctor of Engineering coursework, they also provided invaluable experience for the intern to understand research and development projects in a business environment, which was one of the major objectives of the internship

Opportunities for improving irrigation efficiency with quantitative models, soil water sensors and wireless technology

Increasingly serious shortages of water make it imperative to improve the efficiency of irrigation in agriculture, horticulture and in the maintenance of urban landscapes. The main aim of the current review is to identify ways of meeting this objective. After reviewing current irrigation practices, discussion is centred on the sensitivity of crops to water deficit, the finding that growth of many crops is unaffected by considerable lowering of soil water content and, on this basis, the creation of improved means of irrigation scheduling. Subsequently, attention is focused on irrigation problems associated with spatial variability in soil water and the often slow infiltration of water into soil, especially the subsoil. As monitoring of soil water is important for estimating irrigation requirements, the attributes of the two main types of soil water sensors and their most appropriate uses are described. Attention is also drawn to the contribution of wireless technology to the transmission of sensor outputs. Rapid progress is being made in transmitting sensor data, obtained from different depths down the soil profile across irrigated areas, to a PC that processes the data and on this basis automatically commands irrigation equipment to deliver amounts of water, according to need, across the field. To help interpret sensor outputs, and for many other reasons, principles of water processes in the soil–plant system are incorporated into simulation models that are calibrated and tested in field experiments. Finally, it is emphasized that the relative importance of the factors discussed in this review to any particular situation varies enormously

Internet of Things in Agricultural Innovation and Security

The agricultural Internet of Things (Ag-IoT) paradigm has tremendous potential in transparent integration of underground soil sensing, farm machinery, and sensor-guided irrigation systems with the complex social network of growers, agronomists, crop consultants, and advisors. The aim of the IoT in agricultural innovation and security chapter is to present agricultural IoT research and paradigm to promote sustainable production of safe, healthy, and profitable crop and animal agricultural products. This chapter covers the IoT platform to test optimized management strategies, engage farmer and industry groups, and investigate new and traditional technology drivers that will enhance resilience of the farmers to the socio-environmental changes. A review of state-of-the-art communication architectures and underlying sensing technologies and communication mechanisms is presented with coverage of recent advances in the theory and applications of wireless underground communications. Major challenges in Ag-IoT design and implementation are also discussed