Next Article in Journal
Effect of Simulated Tillage in Combination with Post-Shattering Temperature Conditions on Senna obtusifolia and Xanthium strumarium Seed Survival, Seedling Emergence and Seedbank Potential
Next Article in Special Issue
Factors Affecting Collective Actions in Farmer-Managed Irrigation Systems of Nepal
Previous Article in Journal
Characterisation of Castor (Ricinus communis L.) Seed Quality Using Fourier Transform Near-Infrared Spectroscopy in Combination with Multivariate Data Analysis
Previous Article in Special Issue
Harvesting Method Affects Water Dynamics and Yield of Sweet Orange with Huanglongbing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling Soil Water Content in a Tomato Field: Proximal Gamma Ray Spectroscopy and Soil–Crop System Models

1
Department of Physics and Earth Sciences, University of Ferrara, Via Saragat 1, 44121 Ferrara, Italy
2
INFN, Ferrara Section, Via Saragat 1, 44121 Ferrara, Italy
3
INFN, Legnaro National Laboratories, Viale dell’Università 2, 35020 Padua, Italy
4
Consorzio Bonifica CER, Via Masi 8, 40137, Bologna, Italy
5
Department of Agricultural and Food Sciences (DISTAL), University of Bologna, Viale Fanin 44, 40127 Bologna, Italy
6
Servizio Idro-Meteo-Clima di Bologna Agenzia Regionale Prevenzione, Ambiente ed Energia, Via Po 5, 40139 Bologna, Italy
*
Author to whom correspondence should be addressed.
Agriculture 2018, 8(4), 60; https://doi.org/10.3390/agriculture8040060
Submission received: 12 March 2018 / Revised: 6 April 2018 / Accepted: 17 April 2018 / Published: 19 April 2018
(This article belongs to the Special Issue Modelling for Water Management in Agriculture Systems)

Abstract

:
Proximal soil sensors are taking hold in the understanding of soil hydrogeological processes involved in precision agriculture. In this context, permanently installed gamma ray spectroscopy stations represent one of the best space–time trade off methods at field scale. This study proved the feasibility and reliability of soil water content monitoring through a seven-month continuous acquisition of terrestrial gamma radiation in a tomato test field. By employing a 1 L sodium iodide detector placed at a height of 2.25 m, we investigated the gamma signal coming from an area having a ~25 m radius and from a depth of approximately 30 cm. Experimental values, inferred after a calibration measurement and corrected for the presence of biomass, were corroborated with gravimetric data acquired under different soil moisture conditions, giving an average absolute discrepancy of about 2%. A quantitative comparison was carried out with data simulated by AquaCrop, CRITeRIA, and IRRINET soil–crop system models. The different goodness of fit obtained in bare soil condition and during the vegetated period highlighted that CRITeRIA showed the best agreement with the experimental data over the entire data-taking period while, in presence of the tomato crop, IRRINET provided the best results.

1. Introduction

In the context of current global warming and uncertainty about future climate conditions, variation in rainfall amounts and dry spells frequency are expected to have negative impacts on the vegetation water availability and consequently on crop yields and water productivity [1,2]. Irrigation water, the primary input for agriculture development, is about two-thirds of the total fresh water assigned to human uses [3,4]. In areas of water scarcity, such as those of the Mediterranean basin, farmers are encouraged to adopt management strategies aimed at reducing water wastes. Significant improvements typically come from advanced technological systems aimed at optimizing water uptake both in terms of amounts of water required for the specific crop and of irrigation frequency [5]. Soil–crop systems models, designed for practical applications and decision support, are useful tools for providing a real-time estimate of the amount of water available to the crop [6]. Nevertheless, the high temporal variation of soil properties, weather, and environmental conditions makes the modelling of soil water dynamics complex and problematic if not assisted with the running of expensive and time-consuming field monitoring of inter-temporal and site-specific parameters.
The continuous determination of soil moisture dynamics at field scale with non-invasive and contactless measurement techniques is at the time the most promising challenge for optimizing agricultural management and in particular for a sustainable use of water [7]. Among the proximal sensing techniques, gamma ray spectroscopy is recognized as one of the best space–time trade off methods [8,9]. Thanks to the permanent installation of measurement stations, soil water content can be estimated in real time on the basis of temporal changes in gamma ray intensity to which soil moisture is negatively correlated. Although the possibility of performing non-destructive and fast measurements has been investigated in the past decades [10,11,12,13], recent developments in cyber-physical systems boosted by investments in Industry 4.0 sectors and supported by the cloud technology are opening new perspectives in the field. In this scenario, the gamma ray spectroscopy method excellently matches the current needs for high accuracy nonstop soil water content monitoring representing a joining link between punctual and satellite fields of view. Indeed, continuous time series of water content at intermediate spatial scales are usually lacking with traditional techniques. From one side, time domain reflectometry [14] and remote sensing such as microwave scatterometers [15] provide the temporal evolution of soil water budget at punctual or catchment scale. On the other side, ground penetrating radar [16] and electrical resistivity tomography [17] act at the field spatial scale with the limitation of providing data referred only to the time of the survey.
The aim of this study is to investigate the potentialities of proximal gamma ray spectroscopy for a real-time and continuous monitoring of soil water content in the framework of an ad hoc experiment. A permanent measurement station was installed in a tomato field in the Emilia-Romagna region (Italy) and remotely controlled during the entire seven months data-taking period. Tomatoes are among the most water-intensive vegetable crops and this region represents the largest hub for its production in Italy [18]. Limited irrigation capacity, as well as water volumes excesses, can negatively influence the development of plants and fruits and provokes different kinds of diseases. Experimental daily values of soil water content, estimated on the basis of gamma-ray measurements, were compared with the results obtained with three different soil–crop system models: (i) CRITeRIA, a physically-based numerical model for soil water balance and crop systems, developed by the Regional Environmental Protection Agency (ARPA) of the Emilia-Romagna region [19]; (ii) AquaCrop, the Food and Agriculture Organization (FAO) crop model simulating the yield response to water [20]; and (iii) IRRINET, a system for irrigation scheduling developed by the Emiliano-Romagnolo Canal (CER) irrigation consortium [21].

2. Materials and Methods

2.1. Experimental Site

The proximal gamma ray sensing experiment was conducted in the period 4 April–2 November 2017 in a 40 × 108 m tomato test field of the Acqua Campus (44.57° N, 11.53° E; 16 m above sea level), a research Centre close to Bologna of the CER irrigation district in Emilia-Romagna [22] (Figure 1). The climate of the experimental area is classified according to the Koppen–Geiger classification as Cfa [23], i.e., a temperate climate, without dry seasons and with hot summers, with a mean temperature of 14.0 °C and mean annual precipitation of about 700 mm.
The main physical and hydraulic parameters of the soil, characterized by a loamy texture, are reported in Table 1. The analysis by the sieving and the hydrometer method followed the standard procedure described in [24] and was performed with the aim of measuring the particle size distribution (PSD) function and, therefore, the percentage of silt, clay, and sand (Table 1). The soil classification was based on the United States Department of Agriculture (USDA) method. The PSD was used to derive the soil water retention and hydraulic conductivity functions.
After having experimentally obtained the PSD, we employed a mono-modal log-normal distribution, where the particle diameter is replaced by its natural logarithm. Details about the methodology are provided in [25]. The µ = 22.07 µm mean and a σ = 2.81 µm standard deviation show that the distribution is within the silt range, when using the USDA textural classifications limits [25]. The obtained PSD parameters were adopted for the retrieving the parameters of the Campbell’s equation [26]:
ψ e = 0.61 ln μ 3.9
b = 8.25 1.26 ln μ
where ψ e is the air entry potential (J/kg), b is the slope parameter in the Campbell’s water retention curve equation [27]. The saturated water content was obtained from the relationship θ s = 1 ρ b ρ s , where ρ b is the measured bulk density (Table 1), and ρ s is the particle density (a value of 2600 kg/m3 was assumed). The Campbell’s water retention curve (Figure 2) for the soil water content θ as function of the ψ potential was then obtained as:
θ ( ψ ) = { θ s ( ψ m ψ e ) 1 b ,   if   ψ m ψ e θ s ,   if   ψ m < ψ e
The soil hydraulic properties in terms of wilting point (θWP), field capacity (θFC), and saturated hydraulic conductivity (Ks) (Table 1) were then inferred from the retention curve of Figure 2.
Tomato plants were transplanted on the 23 May in coupled rows: the plants geometry was such that the distance between two couples was 1.5 m, the distance between two rows of the same couple 0.45 m and the distance between two tomato plants of the same row 0.38 m. On the basis of this configuration, and assuming a homogeneous plant distribution on the experimental site area, a 3.5 plants/m2 planting density was estimated. The anthesis of the plants took place approximately on the 9 June, the berries maximum maturity occurred on the 30 August and the harvesting on the 14 September (Figure 3a). The crop irrigation, performed via sprinklers (Figure 1), was scheduled following the criteria provided by the decision support tool of IRRINET (see Section 2.3.3).

2.2. Experimental Setup

The experimental setup comprised a commercial agro-meteorological station (MeteoSense 2.0, Netsens) and a gamma ray spectroscopy station specifically designed and built for the experiment (Figure 4).
Measured weather data included air temperature, relative air humidity (RH), wind direction and speed, precipitation, and short wave incoming radiation (SWIR), the latter measured with a silicon-photodiode pyranometer. We report for the entire data-taking period the temporal profiles of daily values recorded for the minimum (Tmin) and maximum temperatures (Tmax) (respectively ranging in the Tmin = (1.3–22.7) °C and Tmax = (13.5–39.3) °C intervals) (Figure 3a), the SWIR (ranging from 34.7 to 257.3 W/m2) (Figure 3b), the RH (ranging from 44.3% to 95.2%), the rainfall amount (which reached a maximum of 56.2 mm), the irrigation water (which reached a maximum of 35 mm) (Figure 3c). The reference evapotranspiration (ET0) (Figure 3b) was calculated through the Hargreaves equation [28] by adopting the weather parameters recorded during the data-taking by the agro-meteorological station.
A permanent gamma station, specifically designed and built for the experiment, with a 1 L sodium iodide (NaI) crystal was placed inside a steel box mounted on top of a 2.25 m high steel pole [29]. The gamma spectrometer was coupled to a photomultiplier tube base which output is processed by a digital multi-channel analyzer (MCA, CAEN γstream); the whole system was powered by a solar panel (Figure 4). The detector measured the photon radiation produced in the decays of natural occurring radionuclides (40K, 238U, and 232Th) and recorded a gamma spectrum, i.e., a histogram representing the energy distribution of photons emitted by the source. Since each gamma decay has a specific emission energy, a gamma spectrum is characterized by the presence of distinctive structures (photopeaks) which allow for the identification and quantification of 40K, 238U, and 232Th abundances in the soil. The integrated number of events inside the energy ranges associated to the main photopeaks [30] were used to determine the counts per second (cps), associated to 40K, 238U, and 232Th activities in the soil. The statistical uncertainty on the measured cps for a gamma spectrum with a temporal length of one hour is typically lower than 1%.
The lateral and vertical horizons of proximal gamma-ray spectroscopy can be defined on the basis of the probability law that governs photon survival in traversing a given material, as reported in [31]. As photon propagation is ruled by the density of traversed materials, a gamma spectrometer is sensitive to about 25 m far in the horizontal direction (Figure 4) and approximately 30 cm deep in the soil (Figure 5).
Weather data were recorded by the agro-meteorological station with an average frequency of about five minutes, while the gamma station provided a list mode output, i.e., a continuous logging of individual photons arrival time and energy. Therefore, a dedicated software was developed to group and synchronize all the acquired data in a single time-referenced dataset. As both stations were equipped with a GPRS connection, it was possible to remotely preprocess the data in real time. The resulting dataset had a temporal resolution of 15 min and merged 10 data fields from the agro-meteorological station temporally aligned with 34 data fields from the processed energy spectra. Both stations were operative during the entire data taking period, for a 94.8% overlapping duty cycle and a 260 GB global amount of uncompressed data, the major part being raw gamma spectra.

2.3. Soil Water Content with Soil Crop System Models

The daily values of soil water content simulated with CRITeRIA, AquaCrop, and IRRINET soil–crop system models were compared to the experimental daily values inferred from proximal gamma ray spectroscopy. For the simulation we used the soil properties reported in Table 1, daily temperatures and precipitation amounts measured by the agro-meteorological station, scheduled irrigations, growth stages of the tomato crop and reference evapotranspiration (ET0) calculated on the basis of the Hargreaves method [28] (Figure 3). The CRITeRIA and AquaCrop values, referred to the entire data-taking period (4 April–2 November 2017), were obtained by running the models with a six month spin-up which was performed by adopting meteorological data published at [32] in order to adjust for initial conditions [25]. As the IRRINET simulation is strictly bounded to the crop development, soil water content data are referred to the sole tomato crop season (23 May–14 September 2017).

2.3.1. CRITeRIA

CRITeRIA is a suite of soil water balance and crop modelling systems developed by the Regional Environmental Protection Agency (ARPA) of the Emilia-Romagna region (Italy). Here we used CRITeRIA-1D, which is a one-dimensional model simulating the soil water balance, the nitrogen balance and crop development. The CRITeRIA-1D model simulates soil water movement by using a numerical solution of Richards’ equation as described in [19]. The model implements many herbaceous crops and fruit trees and requires as input daily weather data (temperature and precipitation), soil texture and hydraulic properties, bulk density and crop management information. For this study, water retention and hydraulic conductivity curves were determined on the basis of PSD experimental data (see Section 2.1).
CRITeRIA implements a crop growth model based on the day degree sum. The relevant variables are the development of the leaf system (expressed by the leaf area index parameter) for the epigeal part and a root growth model for increase and spatial distribution of the root system. As the tomato season was characterized by anomalously high temperatures, parameters regulating the leaf area index increase and decrease were properly calibrated to follow the effective tomato plant growth at the experimental site.

2.3.2. AquaCrop

AquaCrop is a crop water productivity model implemented by the Land and Water Division of FAO [20]. In the perspective of being applied to diverse worldwide agricultural systems characterized by large crop and soil variability and by the unavailability of extended number of site specific input variables, AquaCrop is structured as an almost ready to use simulation model, requiring a low number of explicit parameters.
The model simulates the water balance and yield response to water of herbaceous crops and it is typically used in situations where water is a key limiting factor in crop production [20]. AquaCrop develops a structure (sub-model components) that includes: (i) the soil, with its water balance; (ii) the crop, with its development, growth, yield and management; (iii) the atmosphere, with its thermal regime, rainfall, evaporative demand, and carbon dioxide concentration. The calculation procedure is grounded on basic and complex biophysical processes to guarantee an accurate simulation of the response of the crop in the plant–soil system. The computation of the soil water transport is based on a tipping-bucket conceptual model employing soil hydraulic properties in terms of θFC and θWP.

2.3.3. IRRINET

IRRINET is one of the tools provided to farmers in the framework of the Emilia-Romagna Action Plan for Rural Development 2007–2013. It is a model for irrigation management backed by the results of more than 50 years of research on plant/water relation and sustainable irrigation management [21]. IRRINET was developed with the aim of progressively reducing water use for irrigation without harming farmers’ income, therefore saving water and optimizing water productivity. The freely available service is web and GIS based and provides decisional criteria for irrigations for a large number of water demanding crops. The model can be employed by using several data sources as meteorological and soil data from local services and crop parameters as defined by CER, including application of the most effective crop tailored irrigation strategy.
Since 2009, IRRINET implements economic calculation of the irrigation profitability assessing the economic benefit related to the next irrigation. Users are provided with optimal irrigation volume and interval, via web or mobile phone text message.
The crop water balance is calculated at field scale with a daily frequency and referred to the crop characteristics, simulated, or input by the farmer. The model structure deals with the soil–plant–atmosphere continuum and includes (i) the soil water balance with capillary rise [33], infiltration rate and run off [34]; (ii) the development of plant with crop coefficients [35]; (iii) the atmosphere (thermal regime, rainfall, and evaporative demand); and (iv) the irrigation system adopted. For sprinkler systems, IRRINET guides the irrigation in order to maintain soil moisture between 30% and 70% of the available water (θFC–θWP). Mean soil water content is estimated in a progressively deeper soil range which linearly follows the growth of the root system from an initial depth of 15 cm (at the time of transplanting) up to 65 cm.

3. Results

The soil water content was determined on the basis of proximal gamma-ray spectroscopy for the entire seven month data-taking period. The basic principle for the determination of soil water content using natural terrestrial gamma radiation recorded by a permanent measurement station consists in the quantification of relative differences between gamma signals measured under different moisture level conditions. While 208Tl (232Th decay chain) and 40K are distributed solely in the soil, gamma radiation originated by the decay of 214Bi (238U decay chain) comes both from 214Bi in the soil and from 214Bi in the atmosphere produced by the decay of 222Rn gas exhaled from rocks and soils. For this reason, we chose 40K as natural gamma emitter for soil water content assessment purposes, given also the highest net counting statistics in its photopeak.
Since the water mass attenuation coefficient is significantly higher than those of typical minerals commonly present in the soil, a gamma spectroscopy measurement is extremely sensitive to different soil moisture levels. The gamma photon flux at the soil–air interface and, consequently, the measured gamma signal are inversely proportional to soil water content: by simultaneously performing an independent gravimetric calibration measurement and a radiometric acquisition, it is possible to evaluate the gravimetric soil water content wγ (kg/kg) at a specific time t as [29,36]:
w γ K ( t ) = S K Cal S K ( t ) × [ Ω + w Cal ] Ω
where wCal is the gravimetric water content at calibration time, SKCal and SK(t) are the gamma signals in cps recorded in the 40K energy window respectively at calibration time and at time t.
The calibration measurement wCal = (0.163 ± 0.008) kg/kg is the mean gravimetric water content estimated after a gravimetric campaign performed in the experimental field in bare soil condition in the 0–30 cm depth range at 16 planar sampling points homogeneously distributed within 15 m from the gamma station. The gravimetric measurements of the calibration process are described in detail in [36]. The gamma calibration measurement SKCal refers to a 2 h late morning data acquisition (10.00 a.m.–12.00 a.m.) concomitant with the gravimetric sampling, distant from rainfall and irrigation events and with stable atmospheric conditions. The adimensional coefficient Ω is a constant determined on the basis of the ratio between the mass attenuation coefficient in solids and water. In this study, we adopted a value of Ω = 0.899, estimated on the basis of the specific composition of the soil of the experimental site [30]. If a soil mineralogical analysis lacks, a mean value together with its standard deviation Ω = (0.903 ± 0.011) estimated in [29] can be adopted. The Ω variability is expected to have second order effects on the estimation of gravimetric water content as the variation introduced by different Ω values would be undistinguishable within the ~0.017 kg/kg absolute uncertainty on the estimated w values, which is essentially dominated by the systematic component associated to the Scal and wCal calibration reference values [36].
While in bare soil condition gamma rays propagate only through air after having reached the soil surface, the presence of growing vegetation introduces a sizable extra attenuation that can be modelled in terms of biomass water content (BWC), i.e., an equivalent water layer which thickness varies in time as the crop evolves during its life-cycle. In order to avoid a systematic overestimation of the gravimetric water content, a correction for the BWC presence needs to be applied. The shielding effect produced by the crop system can be estimated by modeling stems, leaves and fruits as an equivalent water layer characterized by a given thickness which we express as a BWC in units of mm. The gravimetric water content corrected for the attenuation to the vegetation w γ K Λ at time t [36] is given by:
w γ K Λ ( t ) = S K Cal × Λ K ( BWC ( t ) ) S K ( t ) × [ Ω + w G Cal ] Ω
where the BWC correction factor Λk (cps/cps) changes in time and it is given by the ratio between the SK recorded at the time t when there was a given BWC at the experimental site and the signal recorded for BWC = 0. In this study, we adopted the Λk curve determined in [36] by adopting the Monte Carlo simulation method described in [29]. Λk is identically equal to 1 for null vegetation cover, as expected, and is equal to 0.88 for the 9.7 mm maximum BWC estimated during the field experiment. If neglected, this missing correction would have biased the gravimetric water content by approximately a factor of 2.
The hourly average results in terms of 40K gamma count rates (in cps) and of volumetric water content θγ (in m3/m3) are reported in Figure 6. The latter was obtained by multiplying the gravimetric water content wγ for the experimental site soil bulk density (Table 1). As expected, the gamma signal and the volumetric water content are respectively negatively and positively correlated with the amount of precipitations that include both rainfall and irrigation water.
The correction for the BWC and the overall reliability of the method were tested with an ad hoc sampling campaign performed in presence of biomass and in three different soil moisture conditions (Table 2), i.e., one day before an irrigation event (24 July), one (26 July) and three days (28 July) after the same event. One set of additional validation measurements was carried out also in bare soil condition two days after a rainfall event (21 September) [36]. The gamma water contents θγ, inferred from gamma signals measured in the same time interval of samples collection, were compared with the weighted average volumetric water contents θG inferred from the sampling campaign performed with the same scheme of the calibration process. The 0–10 cm, 10–20 cm, and 20–30 cm θG values were combined with weights respectively equal to 0.79, 0.16, and 0.05, determined on the basis of the gamma signal depth profile.
An average discrepancy between the experimental measurements (θG) and estimated values (θγ) of −2.1% is observed. By taking into account the statistical dispersions of the sets of gravimetric measurements, the θγ values are compatible with θG values at 1σ level for all the validation days. Systematic errors leading to underestimations or overestimation of the soil water content are to be excluded also in the presence of the tomato crop at the experimental site.

4. Discussions

As soil water content values provided by CRITeRIA (θC), AquaCrop (θA), and IRRINET (θI) are referred to a daily frequency, we averaged the 24 hourly experimental volumetric water contents (θγ) in order to obtain soil moisture dynamics curves with the same temporal resolution (Figure 7). θC and θA were determined by combining the 0–10 cm, 10–20 cm, and 20–30 cm values according to the 0.79, 0.16, and 0.05 weights described in Section 3. θI refers to the mean value in the variable depth horizon modelled by IRRINET (see Section 2.3.3). The four datasets show consistent temporal trends, positively correlated with the amount of precipitated water (Figure 3c).
Although the three simulated datasets never reach the saturation value θs (Table 1), we observed different variability ranges that reflect the different computational methods of water depth redistribution. The widest dynamic observed for the θA values is attributable to the tipping-bucket conceptual model implemented in AquaCrop, making the water content of the most superficial layer exceed field capacity θFC (Table 1) and reach a maximum value (0.44 m3/m3) close to θs before draining water to the deeper layers. In the summer period, when the depletion of water content due to plant water uptake becomes relevant, AquaCrop provides the lowest values (0.07 m3/m3), which are below the wilting point θWP (Table 1). Conversely, CRITeRIA exhibits a more homogenous depth water distribution as physically based models usually provide a more accurate computation of water fluxes across soil strata. As in CRITeRIA each soil layer is characterized by a limited water content variability (~25%), the excursion of θC values is confined in the (0.11–0.31) m3/m3 interval, almost corresponding to the (θWP–θFC) range. A narrower variability (0.14–0.23 m3/m3) is observed in the temporal trend of IRRINET values (θI), which are intrinsically averaged over a variable maximum depth defined as the tomato root system vertical horizon. The θγ experimental values lie in the (0.12–0.33) m3/m3 range, compatible with soil hydraulic properties in terms of θWP and θFC (Table 1), and are consistent with models variability almost over the entire data taking. Limited temporal misalignments between experimental and simulated daily values can be due to the averaging of gamma inferred volumetric water contents, which are sensitive to the occurring of impulsive rainfall and irrigation events (Figure 6). Conversely, simulation models cluster the output water contents with a daily resolution, making the variations unresolvable due to distinct precipitation events within a day.
The general agreement between experimental and simulated values was assessed in terms of linear regression between the datasets over the entire data-taking period and separately over the bare soil and vegetated periods (Figure 8). Together with the slope (m ± δm), intercept (q ± δq) and coefficient of determination (r2), we estimated the Nash–Sutcliffe coefficient (NS) [37] (Table 3). The NS quantifies the relative magnitude of the residual variance between experimental (θγ and simulated values (θA, θC, and θI) compared to the measured data variance. A good model performance is indicated by a NS value close to 1. When NS = 0 the model predictions are as accurate as the mean of the observed data, whereas for NS < 0 the observed mean is a better predictor than the model, since the residual variance is greater than the data variance. It is a useful indicator for time series as it also captures the model performance over time.
Considering the m, r2, q, and NF values over the entire data taking period, CRITeRIA shows a better agreement with experimental measurements with respect to AquaCrop (Table 3). Indeed AquaCrop performed poorly, especially in presence of tomato crop when a negative NS value (−1.04) is observed and θA values increase on average by a factor of 2 with respect to θγ values. During the vegetated period, when also IRRINET simulated data are available, the best agreement is observed between θγ and θI. Indeed, the goodness of fit between experimental and IRRINET simulated data is proved by a q value compatible with 0 at a 1.5σ level and an average 20% increase of soil water content with respect to experimental values. For CRITeRIA, AquaCrop, and IRRINET we obtain r2 values moderately different from 1 (Table 3) which, however, do not indicate fitting models different from linear trends but highlight a relatively high dispersion of the data around the regression line (Figure 8a–e). The mean and standard deviation of percentage differences between modelled and experimental values are (1 ± 14)% for CRITeRIA, (−9 ± 35)% for AquaCrop, and (−3 ± 16)% for IRRINET. Although for all the models relatively high normalized differences respect to θγ are observed in the vegetated period (Figure 8b–f), the time series of percentage discrepancies show a scattered behavior which leads to exclude evident biases, e.g., due to the correction for BWC attenuation applied to experimental signals.
The volumetric water content values θG determined on the basis of gravimetric campaigns carried out in July and in September (Section 3) represent a corroboration reference for evaluating the performances of simulation models and of the gamma ray spectroscopy method. It is worth underlining that Figure 9 reports θG values representative of the water content of the soil at the time of the gravimetric campaign (~2 h) while θγ, θA, θC, and θI refer to daily moisture levels. In the presence of tomato crops (Figure 9a), the gamma ray spectroscopy method provides the best agreement with gravimetric measurements, characterized by a 9.8% absolute mean relative discrepancy, while all the simulation models show a deviation larger than 15%. Although in bare soil condition (Figure 9b) gravimetric measurements are lower than simulated and experimental values, CRITeRIA and AquaCrop show coherent temporal trends, enclosing the time series inferred from gamma ray measurements.

5. Conclusions

In the perspective of optimizing agricultural management with a sustainable use of water, one of the current challenges is performing continuous and reliable soil moisture monitoring at field scale with non-destructive and real time techniques. In this context, we performed a proximal gamma ray spectroscopy experiment by installing a permanent station in a tomato field which allowed for assessing soil water content dynamics over a seven-month period. The main results of this work are as follows:
(i)
Proximal gamma ray spectroscopy is an excellent method for a non-stop tracing of soil water content at an intermediate spatial scale between punctual and satellite fields of view;
(ii)
Once a reliable calibration is provided through direct measurements, soil water contents inferred from gamma ray spectroscopy do not require detailed soil and crop parameterization and are characterized by relatively low uncertainties;
(iii)
While soil–crop system models simulate soil dynamics with a daily resolution, the proposed method is able to provide reliable higher frequency estimations sensitive to transient soil moisture levels, as proved by the excellent agreement with direct gravimetric measurements;
(iv)
Proximal gamma ray spectroscopy gives a satisfactory description of soil water content over time also when compared to simulation data, showing that the combination of accurate soil water content measurements and water budget computation with crop models can be effective tools for water resources and irrigation planning.

Acknowledgments

This work was partially founded by the National Institute of Nuclear Physics (INFN) through the ITALian RADioactivity project (ITALRAD) and by the Theoretical Astroparticle Physics (TAsP) research network. The authors would like to acknowledge the support of the Project Agroalimentare Idrointelligente CUP D92I16000030009, of the Geological and Seismic Survey of the Umbria Region (UMBRIARAD), of the University of Ferrara (Fondo di Ateneo per la Ricerca scientifica FAR 2016), and of the MIUR (Ministero dell’Istruzione, dell’Università e della Ricerca) under MIUR-PRIN-2012 project. The authors thank the staff of GeoExplorer Impresa Sociale s.r.l. for their support and Renzo Valloni, Enrico Calore, Fabio Sebastiano Schifano, Claudio Pagotto, Ivan Callegari, Gerti Xhixha, and Merita Kaçeli Xhixha for their collaboration which made possible the realization of this study. The authors show their gratitude to Giovanni Fiorentini, Ilaria Marzola, Michele Montuschi, and Barbara Ricci for useful comments and discussions. We would like to thank Carmine Gerardo Gragnano for particle size distribution analysis. The authors thank the Università degli Studi di Ferrara and INFN-Ferrara for the access to the COKA GPU cluster.

Author Contributions

V.S., M.A., M.Ba., and F.M. conceived and designed the work as it is. V.S., M.A., S.A., M.Ba., M.Bi., C.B., E.C., K.G.C.R., D.S., and F.M. realized the experimental setup and performed the measurements with the support of B.F. and V.G. The software for processing gamma data and the spectroscopy analysis was carried out and tuned by M.Ba. and E.C.; the simulations with soil–crop system models were performed by M.Bi., F.T., and G.V. (for CRITeRIA); by M.A. and K.G.C.R. (for AquaCrop); and by S.A. and D.S. (for IRRINET). V.S., M.Ba., and F.M. took the lead in designing and composing the manuscript with the input from all the authors. All the authors critically revised and provided final approval of the version submitted.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Saadi, S.; Todorovic, M.; Tanasijevic, L.; Pereira, L.S.; Pizzigalli, C.; Lionello, P. Climate change and mediterranean agriculture: Impacts on winter wheat and tomato crop evapotranspiration, irrigation requirements and yield. Agric. Water Manag. 2015, 147, 103–115. [Google Scholar] [CrossRef]
  2. Pereira, L.S. Challenges on water resources management when searching for sustainable adaptation to climate change focusing agriculture. Eur. Water 2011, 34, 41–54. [Google Scholar]
  3. Sultana, S.N.; Mani Hazarika, U.; Misra, U.K. Improvement of water use efficiency and remote sensing applications for surface soil moisture monitoring. ADBU-J. Eng. Technol. Sult. 2017, 6, 2348–7305. [Google Scholar]
  4. Patanè, C.; Tringali, S.; Sortino, O. Effects of deficit irrigation on biomass, yield, water productivity and fruit quality of processing tomato under semi-arid mediterranean climate conditions. Sci. Hortic. 2011, 129, 590–596. [Google Scholar] [CrossRef]
  5. Ozbahce, A.; Tari, A.F. Effects of different emitter space and water stress on yield and quality of processing tomato under semi-arid climate conditions. Agric. Water Manag. 2010, 97, 1405–1410. [Google Scholar] [CrossRef]
  6. Liang, H.; Hu, K.; Batchelor, W.D.; Qi, Z.; Li, B. An integrated soil–crop system model for water and nitrogen management in north China. Sci. Rep. 2016, 6, 1–20. [Google Scholar] [CrossRef] [PubMed]
  7. Ochsner, T.E.; Cosh, M.H.; Cuenca, R.H.; Dorigo, W.A.; Draper, C.S.; Hagimoto, Y.; Kerr, Y.H.; Njoku, E.G.; Small, E.E.; Zreda, M. State of the art in large-scale soil moisture monitoring. Soil Sci. Soc. Am. J. 2013, 77, 1888. [Google Scholar] [CrossRef] [Green Version]
  8. Rouze, G.S.; Morgan, C.L.S.; McBratney, A.B. Understanding the utility of aerial gamma radiometrics for mapping soil properties through proximal gamma surveys. Geoderma 2017, 289, 185–195. [Google Scholar] [CrossRef]
  9. Bogena, H.R.; Huisman, J.A.; Güntner, A.; Hübner, C.; Kusche, J.; Jonard, F.; Vey, S.; Vereecken, H. Emerging methods for noninvasive sensing of soil moisture dynamics from field to catchment scale: A review. Wiley Interdiscip. Rev. Water 2015, 2, 635–647. [Google Scholar] [CrossRef]
  10. Carroll, T. Airborne soil moisture measurement using natural terrestrial gamma radiation. Soil Sci. 1981, 132, 358–366. [Google Scholar] [CrossRef]
  11. Grasty, R.L. Radon emanation and soil moisture effects on airborne gamma-ray measurements. Geophysics 1997, 62, 1379–1385. [Google Scholar] [CrossRef]
  12. McHenry, J.R.; Gill, A.C. Measurement of soil moisture with a portable gamma ray scintillation spectrometer. Water Resour. Res. 1970, 6, 989–992. [Google Scholar] [CrossRef]
  13. Loijens, H.S. Determination of soil water content from terrestrial gamma radiation measurements. Water Resour. Res. 1980, 16, 565–573. [Google Scholar] [CrossRef]
  14. Robinson, D.A.; Gardner, C.M.K.; Cooper, J.D. Measurement of relative permittivity in sandy soils using TDR, capacitance and theta probes: Comparison, including the effects of bulk soil electrical conductivity. J. Hydrol. 1999, 223, 198–211. [Google Scholar] [CrossRef]
  15. Brocca, L.; Hasenauer, S.; Lacava, T.; Melone, F.; Moramarco, T.; Wagner, W.; Dorigo, W.; Matgen, P.; Martínez-Fernández, J.; Llorens, P.; et al. Soil moisture estimation through ASCAT and AMSR-E sensors: An intercomparison and validation study across Europe. Remote Sens. Environ. 2011, 115, 3390–3408. [Google Scholar] [CrossRef]
  16. Gerhards, H.; Wollschläger, U.; Yu, Q.; Schiwek, P.; Pan, X.; Roth, K. Continuous and simultaneous measurement of reflector depth and average soil-water content with multichannel ground-penetrating radar. Geophysics 2008, 73, J15–J23. [Google Scholar] [CrossRef]
  17. Alamry, A.S.; van der Meijde, M.; Noomen, M.; Addink, E.A.; van Benthem, R.; de Jong, S.M. Spatial and temporal monitoring of soil moisture using surface electrical resistivity tomography in mediterranean soils. Catena 2017, 157, 388–396. [Google Scholar] [CrossRef]
  18. Tavola Cpom: Superficie (Ettari) e Produzione (Quintali): Pomodoro, Pomodoro da Industria. Dettaglio per Provincia-Anno 2016. Available online: http://agri.istat.it/jsp/dawinci.jsp?q=plCPO0000020000213200&an=2016&ig=1&ct=418&id=15A|18A|28A (accessed on 22 February 2018).
  19. Bittelli, M.; Tomei, F.; Pistocchi, A.; Flury, M.; Boll, J.; Brooks, E.S.; Antolini, G. Development and testing of a physically based, three-dimensional model of surface and subsurface hydrology. Adv. Water Resour. 2010, 33, 106–122. [Google Scholar] [CrossRef]
  20. Steduto, P.; Hsiao, T.C.; Raes, D.; Fereres, E. Aquacrop-the fao crop model to simulate yield response to water: I. Concepts and underlying principles. Agron. J. 2009, 101, 426–437. [Google Scholar] [CrossRef]
  21. Mannini, P.; Genovesi, R.; Letterio, T. Irrinet: Large scale DSS application for on-farm irrigation scheduling. Proced. Environ. Sci. 2013, 19, 823–829. [Google Scholar] [CrossRef]
  22. Munaretto, S.; Battilani, A. Irrigation water governance in practice: The case of the canale emiliano romagnolo district, Italy. Water Policy 2014, 16, 578–594. [Google Scholar] [CrossRef]
  23. Peel, M.C.; Finlayson, B.L.; McMahon, T.A. Updated world map of the Koppen–Geiger climate classification. Hydrol. Earth Syst. Sci. 2007, 11, 1633–1644. [Google Scholar] [CrossRef]
  24. Head, K.; Epps, R.H. Manual of Soil Laboratory Testing, 2nd ed.; CRC Press: Boca Raton, FL, USA, 1992; Volume 1. [Google Scholar]
  25. Bittelli, M.; Campbell, G.S.; Tomei, F. Soil Physics with Python: Transport in the Soil-Plant-Atmosphere System; Oxford University Press: Oxford, UK, 2015. [Google Scholar]
  26. Campbell, G.; Shiozawa, S. Prediction of hydraulic properties of soils using particle-size distribution and bulk density data. In Proceedings of the International Workshop on Indirect Methods for Estimating the Hydraulic Properties of Unsaturated Soils, Riverside, CA, USA, 11–13 October 1989; University of California: Riverside, CA, USA, 1992; pp. 317–328. [Google Scholar]
  27. Campbell, G.S. A simple method for determining unsaturated conductivity from moisture retention data. Soil Sci. 1974, 117, 311–314. [Google Scholar] [CrossRef]
  28. Hargreaves, G.H.; Samani, Z.A. Reference crop evapotranspiration from temperature. Appl. Eng. Agric. 1985, 1, 96–99. [Google Scholar] [CrossRef]
  29. Baldoncini, M.; Albéri, M.; Bottardi, C.; Chiarelli, E.; Raptis, K.G.C.; Strati, V.; Mantovani, F. Investigating the potentialities of monte carlo simulation for assessing soil water content via proximal gamma ray spectroscopy. J. Environ. Radioact. 2018, in press. [Google Scholar]
  30. Nuclear Energy Department, IAEA. Guidelines for Radioelement Mapping Using Gamma Ray Spectrometry Data; Nuclear Fuel Materials Section, IAEA: Vienna, Austria, 2003; p. 179. [Google Scholar]
  31. Feng, T.; Cheng, J.; Jia, M.; Wu, R.; Feng, Y.; Su, C.; Chen, W. Relationship between soil bulk density and pvr of in situ γ spectra. Nuclear Instrum. Methods Phys. Res. Sect. A 2009, 608, 92–98. [Google Scholar] [CrossRef]
  32. Dext3r. Available online: http://www.Smr.Arpa.Emr.It/dext3r/ (accessed on 27 February 2018).
  33. Battilani, A.; Mannini, P. Influence of water table depth on the yield and quality of processing tomatoes. In Proceedings of the International Society for Horticultural Science (ISHS), Leuven, Belgium, 14–18 July 1994; pp. 295–298. [Google Scholar]
  34. Driessen, P.M. The Water Balance of the Soil in Modelling of Agricultural Production: Weather, Soil and Crops; Van keulen, H., Wolf, E.J., Eds.; Center Agricultural Pub: Wageningen, The Netherlands, 1986. [Google Scholar]
  35. Doorenbos, J. Land and Water Division Guidelines for Predicting Crop Water Requirements; FAO Irrigation and Drainage Paper; FAO: Rome, Italy, 1977. [Google Scholar]
  36. Baldoncini, M.; Albéri, M.; Bottardi, C.; Chiarelli, E.; Raptis, K.G.C.; Strati, V.; Mantovani, F. Biomass water content effect in soil water content assessment via proximal gamma ray spectroscopy. Geoderma 2018, in press. [Google Scholar]
  37. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
Figure 1. We report the geographic location of the experimental site of the Acqua Campus of CER (Emilia-Romagna, Italy), the dimensions of the tomato test field and the positions of the gamma (γ) and agro-meteorological (w) stations (panel a). The tomato crop with focuses on a tomato plant and on a sprinkler are shown (panel b).
Figure 1. We report the geographic location of the experimental site of the Acqua Campus of CER (Emilia-Romagna, Italy), the dimensions of the tomato test field and the positions of the gamma (γ) and agro-meteorological (w) stations (panel a). The tomato crop with focuses on a tomato plant and on a sprinkler are shown (panel b).
Agriculture 08 00060 g001
Figure 2. Campbell’s water retention curve obtained on the basis of the mean particle size determined after the parameterization of the PSD.
Figure 2. Campbell’s water retention curve obtained on the basis of the mean particle size determined after the parameterization of the PSD.
Agriculture 08 00060 g002
Figure 3. Weather parameters recorded in the 4 April–2 November 2017 data-taking period by the agro-meteorological station installed at the experimental site: maximum (Tmax) and minimum (Tmin) temperature (panel a), short wave incoming radiation (SWIR) (panel b) and relative air humidity (RH), rainfall and irrigation water amount (panel c). The reference evapotranspiration (ET0) (panel b) is calculated on the basis of Hargreaves equation [28]. The arrows in panel a) indicate the four major crop maturity phases, i.e., planting (P, 23 May), anthesis (A, 9 June), maturity (M, 30 August), and harvesting (H, 14 September).
Figure 3. Weather parameters recorded in the 4 April–2 November 2017 data-taking period by the agro-meteorological station installed at the experimental site: maximum (Tmax) and minimum (Tmin) temperature (panel a), short wave incoming radiation (SWIR) (panel b) and relative air humidity (RH), rainfall and irrigation water amount (panel c). The reference evapotranspiration (ET0) (panel b) is calculated on the basis of Hargreaves equation [28]. The arrows in panel a) indicate the four major crop maturity phases, i.e., planting (P, 23 May), anthesis (A, 9 June), maturity (M, 30 August), and harvesting (H, 14 September).
Agriculture 08 00060 g003
Figure 4. Gamma (γ) and agro-meteorological (w) stations installed at the experimental site. About 95% of the signal received by the gamma spectrometer located at 2.25 m above the ground is produced within a ~25 m radial distance [31].
Figure 4. Gamma (γ) and agro-meteorological (w) stations installed at the experimental site. About 95% of the signal received by the gamma spectrometer located at 2.25 m above the ground is produced within a ~25 m radial distance [31].
Agriculture 08 00060 g004
Figure 5. Percentage contribution to the 40K ground level unscattered photon flux as function of soil depth. The contribution is computed, adopting the theoretical background reported in [31], for four soil density values: as the density increases from 1.2 g/cm3 to 1.8 g/cm3, the depth corresponding to 95% of the signal decreases from 28 cm to 19 cm.
Figure 5. Percentage contribution to the 40K ground level unscattered photon flux as function of soil depth. The contribution is computed, adopting the theoretical background reported in [31], for four soil density values: as the density increases from 1.2 g/cm3 to 1.8 g/cm3, the depth corresponding to 95% of the signal decreases from 28 cm to 19 cm.
Agriculture 08 00060 g005
Figure 6. In panel (a) is reported the 40K gamma signal in cps (SK). In panel (b) the volumetric water content (θγ) estimated on the basis of gamma spectroscopy measurements and corrected for the attenuation due to the biomass water content. The data are hourly averaged. In both panels the daily amount of rainfall and irrigation water are reported in mm.
Figure 6. In panel (a) is reported the 40K gamma signal in cps (SK). In panel (b) the volumetric water content (θγ) estimated on the basis of gamma spectroscopy measurements and corrected for the attenuation due to the biomass water content. The data are hourly averaged. In both panels the daily amount of rainfall and irrigation water are reported in mm.
Agriculture 08 00060 g006
Figure 7. Daily volumetric water content as function of time over the entire data-taking period inferred from gamma ray spectroscopy measurements (θγ), CRITeRIA (θC) and AquaCrop (θA), and over the tomato crop season for IRRINET (θI). The arrows indicate the tomato crop planting (P) and harvesting (H).
Figure 7. Daily volumetric water content as function of time over the entire data-taking period inferred from gamma ray spectroscopy measurements (θγ), CRITeRIA (θC) and AquaCrop (θA), and over the tomato crop season for IRRINET (θI). The arrows indicate the tomato crop planting (P) and harvesting (H).
Agriculture 08 00060 g007
Figure 8. Scatter plot of the daily volumetric water contents inferred from gamma ray spectroscopy (θγ) versus the daily volumetric water content estimated with CRITeRIA (θC) (panel a), AquaCrop (θA) (panel c) and IRRINET (θI) (panel e). For each scatter plot the regression line is superimposed and its corresponding equation is reported. The percentage differences with respect to θγ values are shown for θC (panel b), θA (panel d), and θI (panel f). Different colors are assigned to data points in bare soil condition (brown) and in presence of the tomato crop (orange).
Figure 8. Scatter plot of the daily volumetric water contents inferred from gamma ray spectroscopy (θγ) versus the daily volumetric water content estimated with CRITeRIA (θC) (panel a), AquaCrop (θA) (panel c) and IRRINET (θI) (panel e). For each scatter plot the regression line is superimposed and its corresponding equation is reported. The percentage differences with respect to θγ values are shown for θC (panel b), θA (panel d), and θI (panel f). Different colors are assigned to data points in bare soil condition (brown) and in presence of the tomato crop (orange).
Agriculture 08 00060 g008
Figure 9. Weighted average volumetric water contents derived from the gravimetric campaigns (θG) carried out in the late morning (10 a.m.–12 a.m.) during the vegetated period (panel a) and in bare soil condition (panel b). The θG error bars correspond to the weighted average standard deviation. Daily volumetric water contents (θγ) are reported together with their 1σ uncertainty. Daily simulated values with the CRITeRIA (θC), AquaCrop (θA), and IRRINET (θI) models are also plotted.
Figure 9. Weighted average volumetric water contents derived from the gravimetric campaigns (θG) carried out in the late morning (10 a.m.–12 a.m.) during the vegetated period (panel a) and in bare soil condition (panel b). The θG error bars correspond to the weighted average standard deviation. Daily volumetric water contents (θγ) are reported together with their 1σ uncertainty. Daily simulated values with the CRITeRIA (θC), AquaCrop (θA), and IRRINET (θI) models are also plotted.
Agriculture 08 00060 g009
Table 1. Physical and hydraulic parameters of the experimental site soil for the depth horizon [0–30] cm. Sand, silt, and clay percentage as well as bulk density and organic matter were determined from direct measurements. The wilting point (θWP), field capacity (θFC), saturation value (θs), and saturated hydraulic conductivity (Ks) are inferred on the basis of the methods discussed in Section 2.1.
Table 1. Physical and hydraulic parameters of the experimental site soil for the depth horizon [0–30] cm. Sand, silt, and clay percentage as well as bulk density and organic matter were determined from direct measurements. The wilting point (θWP), field capacity (θFC), saturation value (θs), and saturated hydraulic conductivity (Ks) are inferred on the basis of the methods discussed in Section 2.1.
ParameterValue
Sand (%)45
Silt (%)40
Clay (%)15
Soil textural classLoam
Soil bulk density (kg/m3)1345
Organic matter (%)1.26
Wilting Point (θWP) (m3/m3)0.09
Field Capacity (θFC) (m3/m3)0.32
Saturation (θs) (m3/m3)0.48
Ks (cm/day)23
Table 2. Results of the four corroboration measurements of volumetric water content. Measurements in presence of biomass were performed one day before an irrigation event (24 July), one (26 July) and three days (28 July) after the same event. The fourth measurement (21 September) was performed in bare soil condition two days after a rainfall event. For each measurement we report the volumetric water content inferred from proximal gamma ray spectroscopy measurements (θγ) together with its 1σ uncertainty, the weighted average volumetric content (θG) determined from 16 planar sampling points homogeneously distributed within 15 m from the gamma station, the daily values simulated with the CRITeRIA (θC) and AquaCrop (θA) models; for IRRINET the mean daily volumetric water content (θI) is reported only during the tomato crop season (see Section 2.3).
Table 2. Results of the four corroboration measurements of volumetric water content. Measurements in presence of biomass were performed one day before an irrigation event (24 July), one (26 July) and three days (28 July) after the same event. The fourth measurement (21 September) was performed in bare soil condition two days after a rainfall event. For each measurement we report the volumetric water content inferred from proximal gamma ray spectroscopy measurements (θγ) together with its 1σ uncertainty, the weighted average volumetric content (θG) determined from 16 planar sampling points homogeneously distributed within 15 m from the gamma station, the daily values simulated with the CRITeRIA (θC) and AquaCrop (θA) models; for IRRINET the mean daily volumetric water content (θI) is reported only during the tomato crop season (see Section 2.3).
DayθG [m3/m3]θγ [m3/m3]θC [m3/m3]θA [m3/m3]θI [m3/m3]
24 July0.167 ± 0.0280.170 ± 0.0230.1340.1070.154
26 July0.265 ± 0.0280.243 ± 0.0230.2240.2470.189
28 July0.189 ± 0.0290.179 ± 0.0230.1680.1520.168
21 July0.237 ± 0.0150.245 ± 0.0230.2550.310/
Table 3. Parameters of the linear regressions between the daily volumetric water content inferred from gamma ray spectroscopy measurements and from simulation models. Together with the slope (m ± δm), intercept (q ± δq) and coefficient of determination (r2), the Nash–Sutcliffe efficiency (NS) is reported for three different datasets: in presence of the tomato crop, in bare soil condition, and for the whole data-taking period.
Table 3. Parameters of the linear regressions between the daily volumetric water content inferred from gamma ray spectroscopy measurements and from simulation models. Together with the slope (m ± δm), intercept (q ± δq) and coefficient of determination (r2), the Nash–Sutcliffe efficiency (NS) is reported for three different datasets: in presence of the tomato crop, in bare soil condition, and for the whole data-taking period.
CRITeRIAAquaCropIRRINET
In presence of the tomato cropm ± δm0.68 ± 0.050.37 ± 0.031.21 ± 0.11
q ± δq [m3/m3]0.06 ± 0.010.12 ± 0.01−0.03 ± 0.02
r20.650.530.55
NS0.51−1.040.51
Bare soil conditionm ± δm1.15 ± 0.040.57 ± 0.03/
q ± δq [m3/m3]−0.04 ± 0.010.08 ± 0.01/
r20.870.80/
NS0.810.27/
Whole periodm ± δm0.81 ± 0.040.44 ± 0.02/
q ± δq [m3/m3]0.04 ± 0.010.11 ± 0.01/
r20.720.67/
NS0.69−0.27/

Share and Cite

MDPI and ACS Style

Strati, V.; Albéri, M.; Anconelli, S.; Baldoncini, M.; Bittelli, M.; Bottardi, C.; Chiarelli, E.; Fabbri, B.; Guidi, V.; Raptis, K.G.C.; et al. Modelling Soil Water Content in a Tomato Field: Proximal Gamma Ray Spectroscopy and Soil–Crop System Models. Agriculture 2018, 8, 60. https://doi.org/10.3390/agriculture8040060

AMA Style

Strati V, Albéri M, Anconelli S, Baldoncini M, Bittelli M, Bottardi C, Chiarelli E, Fabbri B, Guidi V, Raptis KGC, et al. Modelling Soil Water Content in a Tomato Field: Proximal Gamma Ray Spectroscopy and Soil–Crop System Models. Agriculture. 2018; 8(4):60. https://doi.org/10.3390/agriculture8040060

Chicago/Turabian Style

Strati, Virginia, Matteo Albéri, Stefano Anconelli, Marica Baldoncini, Marco Bittelli, Carlo Bottardi, Enrico Chiarelli, Barbara Fabbri, Vincenzo Guidi, Kassandra Giulia Cristina Raptis, and et al. 2018. "Modelling Soil Water Content in a Tomato Field: Proximal Gamma Ray Spectroscopy and Soil–Crop System Models" Agriculture 8, no. 4: 60. https://doi.org/10.3390/agriculture8040060

APA Style

Strati, V., Albéri, M., Anconelli, S., Baldoncini, M., Bittelli, M., Bottardi, C., Chiarelli, E., Fabbri, B., Guidi, V., Raptis, K. G. C., Solimando, D., Tomei, F., Villani, G., & Mantovani, F. (2018). Modelling Soil Water Content in a Tomato Field: Proximal Gamma Ray Spectroscopy and Soil–Crop System Models. Agriculture, 8(4), 60. https://doi.org/10.3390/agriculture8040060

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop