Using Cosmic-Ray Neutron Probes to Monitor Landscape Scale Soil Water Content in Mixed Land Use Agricultural Systems
With an ever-increasing demand for natural resources and the societal need to understand and predict natural disasters, soil water content (SWC) observations remain a critical variable to monitor in order to optimally allocate resources, establish early warning systems, and improve weather forecasts. However, routine agricultural production practices of soil cultivation, planting, and harvest make the operation and maintenance of direct contact point sensors for long-term monitoring challenging. In this work, we explore the use of the newly established Cosmic-Ray Neutron Probe (CRNP) and method to monitor landscape average SWC in a mixed agricultural land use system in northeast Austria. The calibrated CRNP landscape SWC values compare well against an independent in situ SWC probe network (MAE = 0.0286 m3/m3) given the challenge of continuous in situ monitoring from probes across a heterogeneous agricultural landscape. The ability of the CRNP to provide real-time and accurate landscape SWC measurements makes it an ideal method for establishing long-term monitoring sites in agricultural ecosystems to aid in agricultural water and nutrient management decisions at the small tract of land scale as well as aiding in management decisions at larger scales.
An accurate understanding of the spatiotemporal dynamics of near surface soil water content (SWC) is essential for a greater understanding of the surface energy balance , the degree of land surface atmospheric coupling [2, 3] and improvements in short-term weather forecasting . With respect to water balance, SWC is a key state variable in determining the partitioning of infiltration and surface runoff  and thus predicting stream hydrograph response and making optimal water management decisions. Moreover, SWC is a key parameter in helping understand and predicting the timing and severity of natural disasters such as drought [6–8] and landslides [9, 10].
Given the importance of SWC observations for understanding energy balance, water balance, and natural disasters, a wide variety of national and state level monitoring networks using point sensors (i.e., Soil Climate Analysis Network, Climate Reference Network, Oklahoma Mesonet, Automated Weather Data Network, etc.) have been established along with remote sensing from airborne (i.e., Airborne Microwave Observatory of Subcanopy and Subsurface) and satellite platforms (i.e., Soil Moisture Ocean Salinity, Soil Moisture Active Passive) [cf. ]. However, a significant gap still exists in fully reconciling differences between point observations with remote sensing  given the natural heterogeneity  and spatial organization SWC fields exhibit .
With the unresolved issues between point sensors and remote sensing, and critical spatiotemporal gaps from remote sensing observations, production agriculture has typically relied on direct insertion point based SWC monitoring technologies to aid in decision-making [15, 16], albeit with recent acceptance of indirect sensors . A key weakness of direct insertion point sensors in production agriculture is often logistical, as routine management practices of planting, soil cultivation, and harvest make installation and maintenance of point sensor networks costly and time consuming to continuously manage.
In this work, we will explore the use of the indirect Cosmic-Ray Neutron Probe (CRNP)  for providing a landscape average SWC value in heterogeneous agricultural landscapes. The CRNP method was recently developed with its main applications in academic research through the establishment of national monitoring networks in the USA , Australia , UK, and South Africa, with probe installations mostly concentrated in natural ecosystems. Here we will investigate the use of the CRNP in a patchy agricultural landscape with mixed winter and summer crops in northeast Austria and compare the landscape average SWC against an independent in situ Time-Domain Transmissivity (TDT) monitoring network. The landscape average SWC value is critical for aiding in both making optimal management decisions of agricultural water and nutrient application at the parcel scale (i.e., a small tract of land) as well as providing data for larger scale management decisions of flood and drought prediction in human dominated landscapes via data integration with modeling and or early warning systems.
The remainder of the paper is organized as follows. First, we will present an overview of the heterogeneous agricultural land use system and soils of the study site in northeast Austria. Next, we will describe the layout of TDT network and briefly summarize the CRNP method for estimating landscape average SWC. We will then compare the landscape average SWC values between the TDT network and CRNP. Finally, we will discuss practical uses and recommendations of using CRNP in agricultural water management systems for operational use and long-term monitoring.
2. Study Area
A CRNP (Model # CRS 1000/B, HydroInnova LLC, Albuquerque, NM, USA) was installed at the study area in northeast Austria (48.1547°N, 15.1483°E, elevation 277 m, Figure 1) on December 11, 2013. The study site, Hydrological Open Air Laboratory (HOAL) , is a cooperation project between the Federal Agency for Water Management (BAW Petzenkirchen) and the Technical University Vienna (TU Vienna), is located in Petzenkirchen, about 100 km west of Vienna and receives an annual average 823 mm of rainfall mostly between April and September. The average annual temperature is 9.5°C. The research station is located in an undulating agricultural landscape, characterized by Cambisols (56%), Planosols (21%), Anthrosols (17%), Gleysols (6%), and Histosols (<1%). Infiltration capacities tend to be medium to low, water storage capacities tend to be high, and shrinking cracks may occur in summer due to high clay contents (see Table A1 and Figure 4 in ). The main crops are winter wheat, barley, maize, and rape. The land use at the study site consists of various parcel sizes making up a patchwork of different crops. Figure 2(a) illustrates the location of the CRNP and the various land use parcels within the CRNP’s measurement area (~28 ha) of a ~300 m radius circle (see Section 3.1). Table 1 summarizes the 2014 planting, soil cultivation, and harvest dates of 11 of the 12 land use parcels within the study area. For full details of the study site, available datasets, overarching research questions, and specific hypotheses, see .
In addition to the weather station and CRNP, a network of Time-Domain Transmissivity (TDT) sensors (SPADE, Julich, Germany) were installed in the second half of 2013. The TDT sensors record hourly SWC at a point and were installed at 31 sites distributed around the study area (Figure 2(b) illustrates the 16 sites within the CRNP measurement area; see Figure 6 in  for full details). At each site 4 TDT sensors were installed horizontally at 4 depths (representing soil layers of ~0–5 cm, 5–10 cm, 15–20 cm, and 45–50 cm). Depending on routine agricultural operations and location of the stations, the TDT sensors are removed at various times throughout the year. The TDT sensors are later reinserted for continued monitoring but the full 2014 data were not available. Instead the network of available TDT sensors (between December 12, 2013, and May 1, 2014) were used to independently compare against the CRNP observations of landscape SWC. We note that given the limited distribution of sensors and spatially varying SWC  that establishing a “true” landscape average SWC is challenging and a comparison against the CRNP should be framed within the expected uncertainty of the mean given the inherent limitations of “spatial representativeness” of averaging a few point sensors in an area.
3. Materials and Methods
3.1. Summary of Cosmic-Ray Neutron Method for Measuring Soil Water Content
The principles and practice of measuring SWC with active source neutrons is well established in soil science and agricultural research [22, 23]. The passive CRNP measures change in the naturally occurring amount of low-energy environmental neutrons to quantify changes in landscape SWC [18, 19, 24]. Because the neutron scattering process at these energies (~1 MeV to 0.5 eV) occurs at high velocities (>10 km/s) , and over tens of meters per collision , the effective radius of measurement is around ~240 m at sea level in dry air [27, 28]. However, the effective CRNP radius is larger at higher elevations (i.e., less air mass means neutrons can travel further during scattering) and smaller with increases in absolute humidity (i.e., more hydrogen in the air reduces travel distance during scattering; see  for full details and equations for calculating footprint). Given the elevation of the study site (277 m.a.s.l.) and absolute humidity measurements, the effective CRNP radius will be ~300 m  and will be adopted for this work herein. Similarly, the effective penetration depth of the CRNP varies from ~15 cm in fully saturated soils (0.40 m3/m3) to ~75 cm in pure silica (SiO2). Full details of the neutron scattering theory, neutron modeling, and coupling to a physically based unsaturated zone model can be found elsewhere [28, 29]. Given the variable methods used for CRNP calibration and validation (i.e., gravimetric, TDR, TDT), variable individual support volumes of calibration method samples (~0–20 cm for TDT and gravimetric and ~0–15 cm for TDR), variable number of TDT sensors used for landscape average, and primary focus of this work for practical applications, we will assume the CRNP has an effective penetration depth of ~20 cm for all observed neutron counts. However, we note that the effective depth is dependent on the depth of the calibration dataset used to parameterize the calibration function (see Section 3.2). Full details and sample calculations on horizontal and depth weighting of CRNP are provided elsewhere .
3.2. Conversion of Observed Neutron Counts into Soil Water Content
The CRNP installed at the site recorded hourly values of moderated neutron counts (counts per hour, cph), atmospheric pressure (hPa), air temperature (°C), and relative humidity (%) (raw data available in real-time at http://cosmos.hwr.arizona.edu/Probes/StationDat/087/index.php). The moderated neutron counts were first corrected for location (i.e., neutron scaling factor), incoming high-energy particles, atmospheric pressure, and absolute humidity following established protocols [19, 30]. The corrected moderated neutron counts were then converted to volumetric pore water content using the calibration function originally proposed by  and further modified by . The calibration function is given bywhere is gravimetric water content (g/g), is lattice water content (g/g), is soil organic carbon water content equivalent (g/g), is the corrected moderated neutron counts per time interval (cph), and is a specific calibrated parameter that represents the count rate over dry silica soils (cph). We note that Soil Water Content (SWC) = volumetric water content, and that (m3/m3), where is the dry soil bulk density (g/cm3) and is the density of water (=1 g/cm3). Following , soil organic carbon water content equivalent can be estimated from on-site soil chemistry sampling aswhere TC is the soil total carbon (g/g), TIC is the inorganic carbon determined by measuring CO2 after the sample is acidified (g/g), 1.724 is a constant to convert total organic carbon into total organic matter, and is the stoichiometric ratio of H2O to organic carbon (assuming organic carbon is cellulose C6H10O5) .
From grid sampling at 50 m resolutions,  found the study site had an average clay weight fraction of 0.20 (g/g), thus yielding g/g by using regression analyses from a catalog of lattice water samples and clay weight fractions (0.03 to 0.38 g/g) from across the globe (Franz unpublished data). We estimated the study site had g/g using a 1 km global SOC product , noting that local samples are also available  but vary widely with land use practice. Lastly, we found from direct sampling at the site on December 12, 2013, the dry soil bulk density was 1.43 g/cm3, with a standard deviation of 0.091 g/cm3, from 61 samples. We note that in agricultural landscapes bulk density may vary significantly across land use and time depending on routine practices. Future work with the CRNP in agricultural settings should address how bulk density changes through time and with land use.
4. Results and Discussion
4.1. Spatiotemporal Variation of Time-Domain Transmissivity Soil Water Content
The time series of daily average SWC between December 2013 and September 2014 for all 16 TDT sites organized by depth are illustrated in Figure 3. Because of varying planting, soil cultivation, and harvest date, TDT profiles from different land use parcels were removed at different dates (Table 1). The TDT probes were later reinserted but the data was not available for the entire time period. Therefore, we restricted our main analyses between December 12, 2013, and May 1, 2014, when 12 of the 16 sites were all available. Figure 3 illustrates the wide range of expected SWC that occur at all 3 measured depths across the various land uses. Figure 4 illustrates the landscape average mean SWC of all sensors by depth, standard error of the mean, and standard deviation. The key points drawn from the TDT network are as follows: () relative changes of TDT response to rainfall across sites are consistent, () estimates of the landscape SWC are uncertain (~0.02 m3/m3 standard error of the mean and 0.07 m3/m3 standard deviation for range of SWC and all soil depths), () absolute values of SWC for a single site are not representative of the landscape SWC for all depths, and () comparison of TDT values for individual sites within a single land use parcel (i.e., S11, 24, 25, 26, 28, 29, and 31 within L4 and S20, 21, and 22 within L15) show a similarly wide range of variability as all TDT sites within the study area. The wide spatial variability of SWC at this site is reported elsewhere using higher density TDR surveys (see Figure in ).
(a) Depth 0–5 cm
(b) Depth 5–10 cm
(c) Depth 15–20 cm
(d) Depth 45–50 cm
(a) Depth 0–5 cm
(b) Depth 5–10 cm
(c) Depth 15–20 cm
(d) Depth 45–50 cm
The observation that a single sensor behaves consistently against the landscape average is consistent with the concept of temporal stability [34, 35] and the ability of the site to be representative of the landscape average. However, this landscape average is not known a priori and requires either a network of sensors to establish how good the mean is (here we found an average of 16 sensors gave a standard error of the mean as ~0.02 m3/m3) or repeated but labor-intensive field campaigns with a portable system such as a TDR. The observation that the intra-land use SWC variability was similar to the inter-land use SWC variability was consistent with other direct point sensor work in agriculture  and has relevance to using characteristics like vegetation type, structure, or leaf area index to select a small number of representative point sensor locations . Due to the fractal nature of SWC fields  and the expected variance at all length scales , point sensors with a small support volume (~0.01 m3) that are located relatively close together (~2 m) may still exhibit very different absolute SWC values due to the natural variation that exists. This implies that the practice of removal and replacement of point sensors following normal production activities can be challenging to recreate a similar SWC response to rainfall, thus requiring a new calibration and new temporal stability analysis to justify the representativeness of the point sensor against the landscape SWC.
4.2. Temporal Variation of Cosmic-Ray Neutron Probe Soil Water Content
Figure 5 illustrates the time series of daily rainfall and moderated neutron counts at the study site. We notice the corrected moderated counts decrease sharply with precipitation and increase slowly following an exponential shape. Using 1 gravimetric calibration and 2 TDR calibrations we estimated using (1). The gravimetric campaign consisted of averaging 61 individual samples, by collecting samples at 3 depths (0–5, 5–10, and 15–20 cm) for 16 locations (where the TDT are located) on December 12, 2013, yielding cph. The TDR campaign consisted of sampling the 16 TDT site locations using a portable device with 15 cm rod lengths inserted vertically from the surface. The TDR calibrations yielded cph on April 5, 2014, and cph on April 30, 2014. The average cph was used to convert moderated neutron counts into SWC for the entire time period. Figure 6 illustrates the time series of daily average SWC from the CRNP, landscape average TDT from different depths, 0–20 cm landscape average TDT, and SWC from the 3 Grav/TDR calibration dates (0–20 cm and 0–15 cm, resp.).
In this work we did not consider any influence of vegetation on the neutron counts given the wide variety of crop types, different timing of planting and harvesting (Table 1), and expected relatively low landscape average fresh standing biomass (<2 kg/m2). For a discussion of the effects of vegetation (i.e., mass changes in above-ground hydrogen) on the expected neutron counts, we refer the reader elsewhere [39–41]. Moreover, given the small changes of values for the 3 different calibration dates (26 cph or 1.8% difference), we found the effects of vegetation on neutron count were likely small for this environment given the accuracy of the gravimetric calibration datasets themselves (standard error of mean 0.02 m3/m3). Finally, the recommendation of using a minimum of 3 calibration sampling dates at different water contents to estimate was reported elsewhere  and is based on neutron particle transport modeling and an error propagation analysis. For best practices, we recommend a minimum of 3 gravimetric calibration periods to estimate in agricultural environments.
4.3. Comparison of Landscape Soil Water Content
In Figure 6, we find that the CRNP compares well against the independent TDT network observations given the standard error of the mean at 0.02 m3/m3 for the TDT landscape average. Most importantly, the CRNP and shallow TDT sensors all respond to precipitation (Figure 5(a)) and decrease at similar rates. Table 2 summarizes a comparison between the 3 TDT and Grav/TDR calibration dates, and estimation of landscape SWC using the CRNP and various calibration datasets. The analysis is based on the spatial averaging and error analysis presented elsewhere . Reference  describes 9 different sampling strategies and error estimates for 3 different study sites around the globe. Here we find that the mean absolute error (MAE) of the landscape average SWC between the 3 TDT and Grav/TDR sampling campaigns is 0.0229 m3/m3 (Table 2). Table 2 also summarizes a cross calibration analysis of the CRNP if individual calibration dates using either TDT or Grav/TDR are used. For a single calibration of , we find that the estimates of SWC from the CRNP can vary widely for different sampling dates (range: −0.074 to 0.0582 m3/m3). However, when 3 calibration dates are used to estimate , we find the range is much narrower (−0.0392 to 0.0354 m3/m3), with MAE = 0.0255 m3/m3 for the TDT data, and 0.0209 m3/m3 for the Grav/TDR samples. Finally, the MAE of the CRNP data using the 3 Grav/TDR calibration datasets versus the independent TDT data is 0.0286 m3/m3. The reduction of error depending on single versus multiple calibration and absolute versus dynamic calibration, and the number of points used in the spatial average is consistent with the findings in . Given the continued investment into the HOAL experiment , we expect future work to more fully describe the various sources of instrument and averaging error.
In terms of the absolute SWC comparison between the Grav/TDR calibrations and the TDT data, we find that the MAE = 0.0286 m3/m3 is comparable to other studies in various natural ecosystems (mixed montane forest , semiarid shrubland , and deciduous forests in the eastern USA  and Germany ) and is on the same order of magnitude as the TDT sensors averaged by depth (Figure 4). The difference between the methods are likely due in part to differences between the TDT factory calibration and local field conditions, uncertainty in the CRNP observations, and natural variability due to spatial averaging heterogeneous environments. Overall the comparison between the TDT network average and CRNP was within acceptable error of <0.04 m3/m3 used in validating remote sensing products against ground observations [45, 46].
In this work, we present a comparison of landscape SWC values between a CRNP and independent TDT sensor network in a mixed agricultural land use system in northeast Austria. Routine agricultural production practices of soil cultivation, planting, and harvest make the operation and maintenance of direct contact point sensor networks challenging. Here, the CRNP footprint consisted of 12 different land use parcels with a mixture of winter and summer crops. From the 3 gravimetric/TDR calibration dates, we found relatively small changes (1.8%) of due to the different land uses and vegetation conditions for this mixed agricultural system. This is consistent with other studies that found minimal changes of vegetation (<2 kg/m2) on expected neutron counts. After using the average from the 3 gravimetric/TDR calibration datasets, we find the MAE of 0.0286 m3/m3 between the CRNP and the independent TDT network is on the same order of magnitude as other studies in natural ecosystems, within the uncertainty for well-established SWC point sensors (0.02 m3/m3 standard error of mean), and within the 0.04 m3/m3 recommendation used in validating remote sensing SWC products against ground observations. For best practice using CNRP, we recommend a minimum of 3 calibration datasets at different wetness levels to estimate instead of relying on a single calibration dataset. Moreover, additional calibration datasets could be used to confirm or locally establish the other three calibration function coefficients given in (1). Because the CRNP can be placed out of the way of routine production practices, we see it as a more suitable method for establishing long-term SWC monitoring in agricultural settings or used in irrigation applications. The ability of the CRNP to provide real-time and accurate landscape SWC measurements in agricultural ecosystems make it an ideal method for aiding in agricultural water and nutrient management decisions at the parcel level as well as aiding in management decisions at larger scales.
Conflict of Interests
Darin Desilets is coowner of HydroInnova LLC. HydroInnova is the only company that commercially produces the cosmic-ray neutron probe. The remaining authors declare that there is no conflict of interests regarding the publication of this paper.
Trenton E. Franz and Ammar Wahbi equally share authorship.
Trenton E. Franz would like to acknowledge the support of the Daugherty Water for Food Institute at the University of Nebraska.
L. S. Dingman, Physical Hydrology, Prentice-Hall, Upper Saddle River, NJ, USA, 2nd edition, 2002.
E. D. Hunt, M. Svoboda, B. Wardlow, K. Hubbard, M. Hayes, and T. Arkebauer, “Monitoring the effects of rapid onset of drought on non-irrigated maize with agronomic data and climate-based drought indices,” Agricultural and Forest Meteorology, vol. 191, pp. 1–11, 2014.View at: Publisher Site | Google Scholar
D. A. Robinson, A. Binley, N. Crook et al., “Advancing process-based watershed hydrological research using near-surface geophysics: a vision for, and review of, electrical and magnetic geophysical methods,” Hydrological Processes, vol. 22, no. 18, pp. 3604–3635, 2008.View at: Publisher Site | Google Scholar
S. Glasstone and M. C. Edlund, Elements of Nuclear Reactor Theory, Van Nostrand, New York, NY, USA, 1952.
D. Desilets, Sandia Report: SAND2011-1101, Radius of Influence for a Cosmic-Ray Soil Moisture Probe: Theory and Monte Carlo Simulations, Sandia National Laboratories, Albuquerque, NM, USA, 2011.
R. Rosolem, W. J. Shuttleworth, M. Zreda et al., “The effect of atmospheric water vapor on the cosmic-ray soil moisture signal,” Journal of Hydrometeorology, vol. 14, no. 5, pp. 1659–1671, 2013.View at: Google Scholar
D. W. Nelson and L. E. Sommers, “Total carbon, organic carbon, and organic matter,” in Methods of Soil Analysis, Part 2, A. L. Page, Ed., pp. 961–1010, American Society of Agronomy, 2nd edition, 1996.View at: Google Scholar
B. B. Mandelbrot, “Stochastic models for the Earth's relief, the shape and the fractal dimension of the coastlines, and the number-area rule for islands,” Proceedings of the National Academy of Sciences of the United States of America, vol. 72, no. 10, pp. 3825–3828, 1975.View at: Publisher Site | Google Scholar | MathSciNet
J. Iwema, R. Rosolem, R. Baatz, T. Wagener, and H. R. Bogena, “Investigating temporal field sampling strategies for site-specific calibration of three soil moisture–neutron intensity parameterisation methods,” Hydrology and Earth System Sciences, vol. 19, no. 7, pp. 3203–3216, 2015.View at: Publisher Site | Google Scholar