Proximal Gamma-Ray Spectroscopy: An Effective Tool to Discern Rain from Irrigation

: Proximal gamma-ray spectroscopy is a consolidated technology for a continuous and real-time tracing of soil water content at ﬁeld scale. New developments have shown that this method can also act as an unbiased tool for remotely distinguishing rainwater from irrigation without any meteorological support information. Given a single detector, the simultaneous observation in a gamma spectrum of a transient increase in the 214 Pb signal, coupled with a decrease in the 40 K signal, acts as an effective proxy for rainfall. A decrease in both 214 Pb and 40 K signals is, instead, a reliable ﬁngerprint for irrigation. We successfully proved this rationale in two data-taking campaigns performed on an agricultural test ﬁeld with different crop types (tomato and maize). The soil moisture levels were assessed via the 40 K gamma signal on the basis of a one-time setup calibration. The validation against a set of gravimetric measurements showed excellent results on both bare and vegetated soil conditions. Simultaneously, the observed rain-induced increase in the 214 Pb signal permitted to identify accurately the rain and irrigation events occurred in the 8852 h of data taking.


Introduction
Due to continuing human population growth and the consequent increasing demand for available water, the management of water resources has become a pressing issue for governments and international organizations. Irrigated agriculture accounts for more than 70% of global water withdrawals [1]. Improvement of farming practices is therefore imperative to ensure availability and sustainable management of water in the future, one of the goals of the 2030 Agenda for Sustainable Development [2].
A concrete solution for an adequate irrigation schedule is the use of a decision support system based on novel technologies for a non-invasive and smart monitoring of soil water content (SWC) [3,4]. The knowledge of the actual SWC represents a strategic element allowing users to manage and use the right amount of water at the right time, leading to an optimal irrigation performance. In a global scenario, a crucial step towards the reduction of water wastage and the increase of water security is the awareness of irrigated lands and of the amount of irrigation water that is really distributed in the soil [5,6].
In this context, the SWC assumes a great importance since it offers double information: how much water is necessary for a growing crop and how much water enters the soil after an irrigation or rainfall event. To exploit these potentialities, it appears evident that knowing which of these two events caused the water precipitation in the soil is mandatory. A unique equipment that measures the field scale SWC and simultaneously distinguishes irrigation and rain events could be the keystone to face these challenges. Proximal gammaray spectroscopy (PGRS) offers the most effective tool to reach these desiderata, recording the gamma signal coming from 40 K and 214 Pb decays to study the evolution of SWC in time and discriminating between rain-and irrigation-induced variations, respectively.
In this paper, we present the results of a two-stage PGRS experiment performed on an agricultural test field in Emilia-Romagna (Italy). In addition to the results obtained for a tomato crop in [7,8], the new dataset of SWC estimated on the basis of gamma data acquired on a maize crop is presented and validated. The PGRS is a novel method that is successfully applied on two different crops adopting a single calibration. For the first time, applying the model developed in [9], the 214 Pb data acquired for both stages of the experiment are used for discriminating rain and irrigation, providing original and reliable results.

Background
The only direct method available for measuring the SWC is currently the gravimetric technique [10], a disruptive and time-consuming procedure requiring operator intervention. Conversely, indirect techniques consist of in situ non-disruptive measurements at a punctual scale (~dm 2 ) or non-invasive remote sensing (RS) methods covering larger areas (up to~km 2 ). PGRS [7,11] is, together with cosmic-ray neutron sensing (CRNS) [12][13][14], a field scale non-invasive technique that fills the aforementioned scale gap. CRNS is based on the measurement of cosmic-ray neutrons that scatter on hydrogen atoms in the ground. These atoms mainly come from the water in the soil, so the CRNS technique can infer the SWC level by studying the scattered neutron flux. Based on the observation of cosmic-ray products, this technique suffers from uncertainties tied to the estimation of the unscattered flux. PGRS is sensitive to the gamma radiation produced by the nuclear decays of 40 K and of radionuclides belonging to 232 Th and 238 U chains, which are naturally present in the soil. The SWC is measured in real time on the basis of the temporal variations of the gamma signal produced by the decays of the 40 K that is naturally and homogeneously (in space and time) distributed in cultivated soils. Since water distributed in the terrain shields the terrestrial gamma flux, an inverse relation between SWC and gamma signal measured by a spectrometric station can be clearly observed and used for a quantitative estimation of SWC [7].
The knowledge of irrigated lands cannot rely on the available datasets since most of the existing comprehensive maps of irrigation water are based only on statistical surveys or can identify just the areas equipped for irrigation, rather than the truly irrigated areas [15]. These maps are not global and are not up-to-date, since illegal pumping is not included; there are self-reporting relevant biases; and temporal/spatial coverage is not homogenous [16,17]. Some vegetation indexes, based on optical and visible RS techniques, are tentatively used as proxies for irrigation monitoring, since irrigated and non-irrigated lands show different spectral responses [18,19]. Sometimes these methods fail because climate conditions, natural vegetation regimes, and agricultural practices are often unpredictable variables. Recent attempts to integrate these approaches seem to give promising results [20], although temporal and spatial heterogeneity of the information remains a persistent Achilles's heel.
In recent years, an alternative approach for global irrigation monitoring has been developed. Diverse algorithms, one of which is SM2RAIN [21], are based on the inversion of the soil-water balance equation to derive the amount of irrigation water, once the rainfall water fraction is identified and subtracted from the total amount of water. PGRS comes into play here: a single detector can retrieve the SWC from the 40 K gamma signal and simultaneously act as an effective proxy for rainfall occurrences. Indeed, PGRS is sensitive to the 214 Pb gamma emitter, a radon progeny concentrated in rain droplets [22,23]. A transient increase in the 214 Pb gamma signal, coupled with a decrease in the 40 K signal, can hence be used as an effective proxy for rain events [9]. Irrigation is, instead, characterized by a decrease in both 214 Pb and 40 K gamma signals.

Rationale
PGRS is a method particularly sensitive to the SWC, since, at typical gamma photon energies (~MeV), the mass attenuation coefficient of water is slightly higher than those of typical minerals commonly present in the soil [24]. As a result, small increments in the SWC can be indirectly estimated by measuring the attenuation in the gamma signal coming from the soil. The detection of the 40 K gamma signal photopeak (1.46 MeV) is convenient in agricultural lands, which are generally rich in this radioisotope. The relation between the 40 K photopeak net count rate C K (t) (in cps) and the soil water content w γ (t) in kg/kg, measured at time t, can be expressed by an inverse proportionality: where C Cal K = 11.7 ± 0.2 (cps) and w Cal g = 0.163 ± 0.008 (kg/kg) are the 40 K net count rate and the gravimetric soil moisture level at calibration time (18 September 2017), respectively [25]. These parameters are site-and detector-dependent and need to be properly measured for an accurate estimation of w γ . The adimensional coefficient Ω = 0.899 represents the ratio between the mass attenuation coefficient of the solid portion of the soil and the one of water, calculated according to the elemental composition of the soil of the experimental site [25]. In the absence of a detailed mineralogical analysis, a Ω = (0.903 ± 0.011) mean value can be employed [25]. The w γ (t) was converted in volumetric soil water content (SWC γ ) in m 3 /m 3 by multiplying it by 1.345, namely the ratio between soil bulk density (Table 1) and water density (∼ 1000 kg/m 3 ).
The adimensional parameter Λ represents the time-dependent count rate attenuation function that accounts for the biomass water content (BWC), which is an attenuation factor of the gamma signal due to the presence of vegetation. Differently from other vegetation indexes, the BWC refers to the amount of water contained in the entire plant (stems, leaves, and fruits). The effect of the BWC on the C K is nontrivial and, in principle, it is indistinguishable from that generated by an increase in the SWC γ . The BWC can be modeled as a layer of some mm of water covering the ground. The presence of a 1.0 mm water layer leads to an overestimation of~10% for the SWC γ value (see Figure 7b of [26]): this attenuation needs to be studied and accounted for in order to avoid systematic errors. The quantification of the gamma signal attenuation as a function of the modeled water layer thickness was studied through Monte Carlo simulations [26]. Since the temporal evolution of BWC is tightly related to crop's growth stage, we estimated the overall BWC on the base of gravimetric measurements on stems, leaves, and fruits samples collected at four and six different stages of tomato and maize maturity, respectively. Modeling the evolution of the entire organism by means of a Gompertz sigmoid function [27], we calculated the time-dependent correction to be applied to the measured gamma signal.
The calibration of the PGRS station was performed only once in 2017 and successfully used for both 2017 and 2020 data-taking campaigns. The value of w Cal g in Equation (1) was obtained from the measurements of 16 sampling points distributed within 15 m from the PGRS station ( Figure 4 of [26]), which is the area where~85% of the signal received by the station originates. For each sampling position, three samples were collected in the depth intervals of 0-10 cm, 10-20 cm, and 20-30 cm for a total number of 48 samples. After drying the samples at 105 • C for about 24 h and inferring their gravimetric water content, the vertical weighted mean for each of the 16 sampling points was obtained by assigning to each depth interval a weight based on its contribution to the gamma signal (calculated according to Equation (3) of [25]). Finally, the value w Cal g was calculated as the mean of the 16 sampling points, with the uncertainty given by their standard deviation, which is of the order of 5%.
As estimated in [26], the uncertainty on w γ (Equation (1)) is dominated by that of w Cal g , rather than by the statistical error of C K (t) (<1%). Hence, the absolute uncertainty on w γ is a constant value equal to 0.017 (kg/kg), corresponding to a 0.023 (m 3 /m 3 ) uncertainty on the estimated volumetric SWC γ values. Since the abundance of 40 K in the ground remains constant in time, due to the inverse relation between the C K signal and w γ expressed by Equation (1), a sudden decrease observed by the PGRS station in C K (t) is related to a SWC increase.
To distinguish between irrigation and rainfall, we focused on the photopeak of 214 Pb (352 keV), a radon progeny belonging to 238 U decay chain. Because of the high volatility of 222 Rn, this gaseous radioisotope easily exhales from the ground and enters the atmosphere. When 222 Rn decays inside a cloud, raindrops capture the aerosol clusters formed by its daughters ( 218 Po and then 214 Pb). During a rainfall, the additional 214 Pb nuclei fallen to the ground cause a sudden increase in the 214 Pb count rate (C Pb ) over the environmental background ( Figure 1). A weak rainfall rate of 5 mm/h causes a C Pb jump from a typical environmental background of ∼1.0 cps to a value of ∼2.2 cps [9]. This increase in C Pb is proportional to the square root of the rain rate R (mm/h), as expressed by the empirically observed relationship: where ∆T = 0.25 h is the temporal bin and the parameter A = 2.15 ± 0.15 (cps mm −0.50 h −0.50 ) is detector dependent [9]. The responses of gamma spectra in the 100-2800 keV energy range changes according to different weather conditions ( Figure 1). In the case of a rain event (Figure 1c), a decrease in the 40 K photopeak is coupled with an increase in the 214 Pb gamma signal ( Figure 1d). During irrigation (Figure 1e), the gamma spectrum ( Figure 1f) does not depict this contrasting behavior, showing a common decrease between the 40 K and 214 Pb photopeaks. Due to the different responses of the 40 K and 214 Pb photopeaks to rainfall and irrigation, it is therefore possible to distinguish between these two events.

Experimental Site
Data were acquired during two different growing seasons: the first data-taking campaign was conducted in the period 4 April-2 November 2017, while the second was carried out in the period 5 March-31 August 2020. During the first campaign (hereafter referenced as T2017), the site was dedicated to a tomato crop, while in the second (M2020) maize was grown on the field. Both data takings refer to the agricultural test field of Acqua Campus, a research center of the Emiliano-Romagnolo Canal (CER) irrigation district in the Emilia-Romagna region of Italy (44.57 • N, 11.53 • E, 16 m above sea level) (Figure 2). The study area is classified, according to the Köppen-Geiger classification [28], in the "Cfa" group, which is characterized by 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. Table 1 lists the main physical and hydraulic parameters of the soil, characterized by a loamy texture and a 1.26% organic matter content [7].
Irrigation water was applied with a fixed sprinkler system in T2017 and with a hose reel irrigator in M2020, following the irrigation schedule supplied by the irrigation advisory service IRRINET [29][30][31]. The employed surface water was drawn from CER, a derivation of the Po river.
The reliability of the PGRS station is proven by its high duty cycle (i.e., the ratio of effective hours of data taking over the total number of hours), which is above 90% for both T2017 and M2020 ( Table 2). The major sources of station downtime were not related to the reliability of the detector, but were, rather, due to the extraordinary maintenance required after extreme weather conditions. Table 1. Physical and hydraulic parameters of the soil of the experimental site for the depth interval of 0-30 cm. Sand, silt, and clay percentage, as well as bulk density and organic matter, were determined from direct measurements. The soil classification is based on the United States Department of Agriculture method. The soil hydraulic properties, i.e., wilting point (θ WP ), field capacity (θ FC ), and saturated hydraulic conductivity (K S ), are derived from [7].

Parameter Value
Sand

Experimental Setup
The experimental setup included a PGRS station and an agro-meteorological station (MeteoSense 2.0, Netsens) (Figure 2), both equipped with internet connection and powered by solar panels.
The PGRS station consisted of an external steel box welded on top of a 2.25 m pole and enclosing a 1 L sodium iodide (NaI(Tl)) gamma-ray spectrometer: the entire station was designed and built specifically for the experiment [25]. The height of the pole was studied to adapt the field of view (FOV) of the PGRS station to the field size. As a result, 95% of the detected gamma signal was produced by a disk that had a radius of approximately 25 m (Figure 2 of [25]), covering, therefore, the entire field width.
The NaI(Tl) crystal was coupled with a photo-multiplier tube base and the output was then processed by a digital multi-channel analyzer (MCA, CAEN γstream) that had 2048 acquisition channels. The MCA was paired with a small integrated computer that ran the software required for managing the acquisition parameters, such as the start time, the spectral gain, and the operating voltage. Additional software was developed and implemented to make the data taking continuous and to remotely inspect the data acquired in real time.
The gamma-ray spectrometer was able to detect the photon radiation produced by the decay of naturally occurring radionuclides ( 40 K, 238 U, and 232 Th decay chains) in the soil. The statistical uncertainty on C K (t), measured in cps, for a gamma spectrum with a temporal length of one hour is typically lower than 1%. The continuously acquired radiometric data was, then, sent every hour to a dedicated server through a 3G internet connection, where an energy calibration was performed. As a result, 15 min and 1 h acquisition time spectra were obtained. From the hourly produced spectra, the server inferred the net count rate in the main 40 K, 214 Pb, 214 Bi ( 238 U), and 208 Tl ( 232 Th) photopeak energy windows [32] and, from their processing, it was able to obtain the volumetric SWC (expressed in m 3 /m 3 ; see Section 3.1). The 15 min spectra were, instead, used for the 214 Pb analysis, with the aim of identifying rain events. All processed data were made available to be queried through a Web API together with the site coordinates and the sensor validity area. The PGRS could thus be easily integrated with any decision support system.

SWC Estimation
The PGRS station was able to continuously measure the SWC γ with hourly resolution during the T2017 and M2020 measurement campaigns. The time evolution of the SWC γ values recorded during the two data-taking periods (Figure 3) highlighted the ability of the PGRS technique to reliably follow the phases of dampening and subsequent drying following rain and irrigation events, both in bare and vegetated soil. The sensor response to water precipitation was immediate: the SWC γ increased as soon as water fell on the ground and decreased with a week-long timescale immediately after water stopped precipitating. The highest value attained by the measured SWC γ during both periods (0.46 m 3 /m 3 ) never exceeded the 0.48 m 3 /m 3 saturation value estimated for the soil, which is consistent with the soil hydraulic and physical properties ( Table 1). The minimum recorded value of 0.09 m 3 /m 3 matched the θ WP (Table 1) exactly and was immediately followed by an irrigation.
A validation campaign performed via gravimetric sampling under different soil moisture and crop conditions was conducted to test the reliability of the adopted PGRS technique and of the BWC correction (Table 3). In T2017, for each of the 4 sets of gravimetric measurements, 48 samples were collected using a soil-auger and the average SWC g was calculated adopting the same procedure used for the calibration of the PGRS station (Section 3.1). During M2020, the validation campaign was extended to 7 sets of gravimetric measurements in which the SWC g was measured on a single bulk sample of the first 40 cm of soil. The obtained 11 SWC g validation measurements were then compared to the SWC γ estimated by the PGRS station in the same dates for both bare soil and vegetated periods ( Table 3 and Figure 4). A single calibration carried out on bare soil (Section 3.1) permitted the PGRS station to perform unbiased measurements of the SWC γ during both data-taking campaigns. The estimated SWC γ and SWC g exhibited an excellent agreement for both T2017 and M2020: all measurements proved to be compatible at 1σ level ( Table 3). The accuracy of the PGRS station appears even more appreciable from the correlation plot between SWC γ and SWC g (Figure 4). The best fit line linking SWC g measurements to SWC γ exhibits a slope and an intercept compatible at less than 1σ level to 1 and 0, respectively. The accuracy of the method extends to both bare soil and tomato-maize covered soil, confirming that, once the BWC is properly modeled, the crop type does not affect SWC γ estimations through PGRS. Table 3. Results of the T2017 and M2020 validation campaigns. For both bare and vegetated soil, the table reports the measured SWC g values together with the corresponding SWC γ values measured by the PGRS station in the same dates. The reported SWC γ correspond to the average SWC γ recorded by the PGRS station during the entire period covering the gravimetric sampling. The uncertainty on SWC γ was set to the constant value of 0.023 m 3 /m 3 , as explained in Section 3.1. The SWC g uncertainties for T2017 samples were calculated according to the procedure used for the calibration of PGRS station (Section 3.1). For M2020 measurements, the relative uncertainty assumed for SWC g data was calculated by averaging the relative standard deviation exhibited by the T2017 validation campaign. The SWC g measurement of 18/09/2017 (in italics) has been used to calibrate SWC γ values according to Equation (1)

Irrigation and Rain Discrimination
The PGRS station offers the ability to independently understand the trigger of water precipitation. In particular, the observation of the time series of C Pb enables an unambiguous identification of irrigation and rain events. This ability appears evident from two significant episodes during T2017 ( Figure 5) and M2020 (Figure 6), when irrigation and rain happened shortly after each other. Whilst the estimated SWC γ increased, correspondingly, during both irrigation and rainfall, the measured C Pb decreased during the former and drastically increased in the latter case. These increases in C Pb , amounting to more than 6 cps, appear clearly distinguishable from both the typical statistical fluctuations in the 214 Pb window (~0.2 cps) and the variability ascribed to well-known day-night oscillations (~0.3 cps) related to atmospheric radon [33][34][35].
Both the estimations of SWC γ and of the rain-induced enhancement of C Pb come from two complementary uses of the same technique, although the time variation of these quantities is related to completely distinct processes characterized by different timescales. After each precipitation, the increase in SWC γ was followed by a week-long decrease caused by the drying of the soil. These very slow variations do not require high temporal resolutions: a 1 h sampling represents an optimal compromise between low statistical errors and a good temporal sampling. The increase of C Pb due to the rainwater precipitation of additional 214 Pb atoms to the ground exponentially decreases following a 26.8 min half-life. A 15 min resolution thus allows for a proper characterization of both the C Pb variation due to the rain event and of the exponential decrease due to 214 Pb decay.  The above illustrated behavior was regularly found during all precipitations recorded during the two data-taking periods (Figures 7 and 8, and Table 4). Regardless of the origin of the precipitated water, the observed SWC γ increased for both rain and irrigation events, with a growth (typically of~0.10 m 3 /m 3 in response to~10 mm of precipitated water) clearly distinguishable from the detected systematic uncertainties (0.023 m 3 /m 3 ).
As soon as rainfall started, C Pb immediately increased (Figures 7a-c and 8a-c), reaching peak values even five times larger than values obtained before rainfall. This increase was proportional to the square root of the rain rate R (Equation (2)) and could not be mistaken for statistical fluctuations, typically of the order of~0.2 cps [9]. This appeared particularly evident from the positive variation of the average C Pb value recorded before and after rains (Table 4), which proved an effective smoking gun for rainfalls. Furthermore, during irrigation events (Figures 7d-f and 8d-f) water precipitation was not followed by any transient increase in C Pb ; the additional water layer precipitated to the soil surface shielded part of the gamma radiation coming from the ground, leading to a negative variation in the average C Pb values recorded before and after irrigation. Figure 7. Temporal evolution of SWC γ (in green) and C Pb (in orange) during three rainfalls (a-c) and three irrigation events (d-f) occurred in T2017, modified after [36]. The SWC γ (in green) is hourly estimated, while the C Pb (in orange) and the water amount (in blue for rain and in gray for irrigated water) have a 15 min time resolution. Figure 8. Temporal evolution of SWC γ (in green) and C Pb (in orange) during three rainfalls (a-c) and three irrigation events (d-f) occurred in M2020. The SWC γ (in green) is hourly estimated, while the C Pb (in orange) and the water amount (in blue for rain and in gray for irrigated water) have a 15 min time resolution. Table 4. Relative variations (in percentage) of C Pb and SWC γ before and after the rains and irrigation episodes shown in Figures 7 and 8, together with their start date and time, duration, and total water precipitated. SWC γ variations are calculated by comparing the average values measured in the 4 h before rain/irrigation and in the 4 h after their end. C Pb variations are calculated considering the 4 h before and the 4 h after rain/irrigation started.  In the absence of precipitation or irrigation water, C Pb statistically fluctuated around the typical background of~1 cps (Figures 7 and 8). However, the background in the 214 Pb energy region is not constant in time: it follows a diurnal and seasonal variation caused by the oscillation of its progenitor 222 Rn content in the atmosphere. This phenomenon is tightly bound to factors affecting 222 Rn mobility, mainly pressure and temperature, and causes a 1-day period oscillation cycle in the observed C Pb (Figure 9). The non-stop measurement of C Pb during the 268 sunny days occurred in T2017, and M2020 allowed a high-statistics characterization of these oscillations. On average, the C Pb variability exhibited in each hour of the day during the entire data taking (i.e., the width of the 1σ band of Figure 9) was 0.33 cps.
Diurnal C Pb oscillations did not have any effect on our ability to identify rain events. Firstly, the discrimination of rains and irrigations is based on the relative increase of C Pb over the background, so the increase remains visible regardless of the before-rain C Pb value. Secondly, these oscillations occur with a period of 24 h, which is much longer than the timescales of rain-induced C Pb variations. Finally, the observed variability of 0.33 cps (obtained in the worst-case scenario by aggregating all data) is lower than rain-induced variations: the C Pb increase due to a rain rate of 1 mm/h is 0.5 cps, while for a rain rate of 5 mm/h is 1.2 cps (Equation (2)). This means that it is still possible to identify a 1 mm/h rain from the background at 1σ, and a 5 mm/h rain at 3σ level.

Conclusions
This study proves the reliability of PGRS as a mature technique for proximal remote sensing applied to precision agriculture. The capability to distinguish rainwater from irrigation without any meteorological support information, paired with the effectiveness in estimating the SWC, makes the PGRS station a consolidated non-invasive technology in supplying crucial information for scheduling irrigations and managing resources with a view to reduce water wastage. In particular, we demonstrated that:

•
The simultaneous observation in a gamma spectrum of a transient increase in the 214 Pb signal, coupled with a decrease in the 40 K signal, is an effective proxy for rainfall, while a decrease in both 214 Pb and 40 K signals is a reliable fingerprint for irrigation; • During a total of 102 rainy days and 23 irrigated days, we were able to discern rain and irrigation without observing any false positive or false negative. Even low rain rates (~1 mm/h) were distinguishable from the gamma background at 1σ level. The raininduced increase in the 214 Pb signal was clearly discernible from both environmental (diurnal oscillations) and statistical fluctuations; • After a single calibration, the PGRS station successfully measured, in real time, the SWC at a field scale level for both tomato (T2017 campaign) and maize (M2020 campaign) crops for a total of~9000 h. Due to the remote-controlled data taking, the PGRS station required on-site maintenance interventions only on a few occasions due to extraordinary weather events.;

•
The accuracy of the PGRS technique was demonstrated through the validation of measurements by comparing SWC γ and SWC g estimates. The results from the two methods proved compatible within 1σ, and the regression line exhibited a slope and an intercept compatible at 1σ level with 1 and 0, respectively. The accuracy extended through bare and vegetated soil conditions, and through different crops (tomato and maize), showing the effectiveness of the correction adopted for the shielding effect of the BWC.
In future research, PGRS stations could be implemented as supporting tools in view of satellite data (e.g., SMAP, SMOS, ASCAT) calibration for quantifying the SWC over large areas, bridging the gap between proximal and remote fields of view. The footprint (field scale) and the high sensitivity of the PGRS make it an essential tool for extracting quantitative information from satellite microwave and reflectance information. Indeed, if on the one hand, new satellite measurements provide an unprecedented free source of information in terms of spatial-temporal resolutions, on the other hand, an effective ground truth of soil moisture products retrieved from satellites is mandatory.  Cosmic-Ray Neutron Sensing, a method of measuring the SWC through the detection of low-energy neutrons FOV Field of View, i.e., the effective range from which the instrument can receive a signal PGRS Proximal Gamma-Ray Spectroscopy, a technique of measuring the SWC through the detection of gamma rays emitted by the decay process of radioactive elements θ WP Wilting point, i.e., the minimum percentage of water content of the soil required by the plant to not wilt θ FC Field capacity, i.e., the maximum water-to-soil ratio that does not trigger water drainage θ S Saturation, i.e., maximum water capacity (in %) of the soil, including water interested by draining K S Hydraulic conductivity, a way of quantifying the ease for a fluid to move in the soil SWC g Volumetric Soil Water Content estimated from gravimetric measurements SWC γ Volumetric Soil Water Content estimated from PGRS measurements w g Gravimetric Soil Water Content estimated from gravimetric measurements w γ Gravimetric Soil Water Content estimated from PGRS measurements C K Gamma count rate produced from the 40 K decay with photopeak at energy 1.46 MeV C Pb Gamma count rate produced from the 214 Pb decay with photopeak at energy 295 keV Ω Parameter determined by the ratio between the mass attenuation coefficients for the solid component of the soil and its water Λ Count rate attenuation function, which expresses the correction (due to the presence of BWC) that must be accounted for during the SWC derivation R Rain rate, i.e., the amount of water accumulated by rain on the ground in a unit of time