Here, we present the protocol for in situ measurement of soil carbon using the neutron-gamma technique for single point measurements (static mode) or field averages (scanning mode). We also describe system construction and elaborate data treatment procedures.
The herein described application of the inelastic neutron scattering (INS) method for soil carbon analysis is based on the registration and analysis of gamma rays created when neutrons interact with soil elements. The main parts of the INS system are a pulsed neutron generator, NaI(Tl) gamma detectors, split electronics to separate gamma spectra due to INS and thermo-neutron capture (TNC) processes, and software for gamma spectra acquisition and data processing. This method has several advantages over other methods in that it is a non-destructive in situ method that measures the average carbon content in large soil volumes, is negligibly impacted by local sharp changes in soil carbon, and can be used in stationary or scanning modes. The result of the INS method is the carbon content from a site with a footprint of ~2.5 – 3 m2 in the stationary regime, or the average carbon content of the traversed area in the scanning regime. The measurement range of the current INS system is >1.5 carbon weight % (standard deviation ± 0.3 w%) in the upper 10 cm soil layer for a 1 hmeasurement.
Knowledge of soil carbon content is required for optimization of soil productivity and profitability, understanding the impact of agricultural land use practices on soil resources, and evaluating strategies for carbon sequestration1,2,3,4. Soil carbon is a universal indicator of soil quality5. Several methods have been developed for soil carbon measurements. Dry combustion (DC) has been the most widely used method for years6; this method is based on field sample collection and laboratory processing and measurement that is destructive, labor intensive, and time consuming. Two newer methods are laser-induced breakdown spectroscopy, and near and mid infrared spectroscopy7. These methods are also destructive and only analyze the very near-surface soil layer (0.1 – 1 cm soil depth). In addition, these methods only yield point measurements of carbon content for small sample volumes (~60 cm3 for DC method, and 0.01-10 cm3 for infrared spectroscopy methods). Such point measurements make it difficult to extrapolate results to field or landscape scales. Since these methods are destructive, recurring measurements are also impossible.
Previous researchers at Brookhaven National Laboratory suggested applying neutron technology for soil carbon analysis (INS method)7,8,9. This initial effort developed the theory and practice of using neutron gamma analysis for soil carbon measurement. Starting in 2013, this effort was continued at the USDA-ARS National Soil Dynamics Laboratory (NSDL). The expansion of this technological application over the last 10 years is due to two main factors: the availability of relatively inexpensive commercial neutron generators, gamma detectors, and corresponding electronics with software; and state of the art neutron-nuclei interaction reference databases. This method has several advantages over others. An INS system, placed on a platform, could be maneuvered over any type of field that requires measurement. This non-destructive in-situ method can analyze large soils volumes (~300 kg) that can be interpolated to a whole agricultural field using just a few measurements. This INS system is also capable of operating in a scanning mode that determines the average carbon content of an area based on scanning over a predetermine grid of the field or landscape.
1. Construction of the INS system
Figure 1. INS System Geometry. Please click here to view a larger version of this figure.
Figure 2. Overview of the INS System.
A) first block contains neutron generator, neutron detector, and power system; B) second block contains three NaI (Tl) detectors; C) third block contains equipment for system operation; D) general view of the first block showing individual components; and E) close up view of the gamma detectors.10 Please click here to view a larger version of this figure.
2. Caution and Personal Requirements
3. Preparation of the INS system for Measurement
4. Calibration of the INS System
Figure 3. View of Pit with Sand and Pit with 10 Cw% Sand-carbon Mixture. Please click here to view a larger version of this figure.
5. Conducting Field Soil Measurements in Static Mode
6. Conducting Field Soil Measurements in the Scanning Mode
Soil INS & TNC and TNC gamma spectra
A general view of the measured soil gamma spectra is shown in Figure 4. The spectra consist of a set of peaks on a continuous background. The main peaks of interest have centroids at 4.44 MeV and 1.78 MeV in the INS & TNC spectra. The second peak can be attributed to silicon nuclei contained in soil, and the first peak is an overlapping peak from carbon and silicon nuclei. The procedure for net carbon peak area extraction from these spectra is described above. This procedure should be used in all cases for determining the net carbon peak area due solely to carbon nuclei.11
Figure 4. A Typical Gamma Spectra for Soil Measured by the INS System. Please click here to view a larger version of this figure.
INS system background measurements
Net INS spectra measured at various system elevation heights above the ground surface are shown in Figure 5.11 Dependencies of the peak areas with centroids at 1.78 MeV, 4.44 MeV, and 6.13 MeV (oxygen peak) with height are illustrated in Figure 6. As shown in this Figure, the spectra no longer change at heights greater than 4 m above the ground surface. Accordingly, the spectra at heights greater than 4 m can be attributed to gamma spectra that appear due to the interaction of neutrons with system construction materials. We used one of these spectra (at H=6 m) as the system background spectrum in our data processing.
Figure 5. a) Net-INS spectra at different INS system heights above the ground; b) Fragment of the net-INS spectra around 1.78 MeV; and c) Fragment of the net-INS spectra around 4.44 MeV. Arrow designates increasing height.11 Please click here to view a larger version of this figure.
Figure 6. Dependencies of Peaks Areas with Centroids at 1.78, and 4.44 MeV in the Net-INS Spectra for INS System with Changing Heights Above the Ground. Please click here to view a larger version of this figure.
Calibration
The net INS spectra generated during INS system calibration are shown in Figure 7a.11 Fragments of the net INS spectra near the 1.78 MeV and 4.44 MeV peaks are shown on a larger scale in Figures 7b and 7c, respectively. As can be seen, the peak with a centroid of 4.44 MeV increases with increasing carbon content in the pit. At the same time, the peak with a centroid at 1.78 MeV slightly decreases as carbon in the pit increases. The dependency of the net carbon peak area (calculated from these spectra) with carbon content in pits (expressed in weight %) is shown in Figure 8.11 As can be seen, this can be represented by a direct proportional dependency passing through the origin (0, 0 point) within experimental error limits. This dependency was used to calibrate further measurements.
Figure 7. a) Net INS spectra for pits with sand-carbon mixtures at 0, 2.5, 5, and 10 carbon w% (uniform mixture); b) Fragment of the net INS spectra around 1.78 MeV; c) Fragment of the net INS around 4.44 MeV.11 Please click here to view a larger version of this figure.
Figure 8. Dependence of the Net Carbon Peak Area with Carbon Concentration in Pits (Points with Error Bars), and INS System Calibration Line (Solid Line).11 Please click here to view a larger version of this figure.
Field measurements of carbon content in static mode
Carbon content measurements in static mode were conducted at several field sites. Results from the Alabama Agricultural Experiment Station Piedmont Research Unit, Camp Hill, AL (110 m x 30 m) are presented in Table 1. Field measurements were conducted at the intersections of a 3 by 5 grid with equal distances between gridlines (total 15 sites). As can be seen from the table, the carbon content for individual intersection points varied between 1.4 to 3.1 w% with the standard deviation of all measurements being ~0.3 w%. For comparison, destructive soil samples were also taken at each location for determining soil carbon content using the standard DC method. These data are also presented in Table 1. Comparison of the two data sets showed good agreement between both methods for each location and for the average value over the whole field.
Location | INS measurements | Dry Combustion measurements | |||||
Site # | |||||||
Carbon, w% | STD, | Plot Average | Carbon, w% | STD, | Plot Average | ||
w% | ±STD, w% | w% | ±STD, w% | ||||
Camp Hill | OF1 | 2.2 | 0.29 | 2.23±0.45 | 2.85 | 0.25 | 2.25±0.51 |
OF2 | 2.51 | 0.29 | 2.54 | 0.31 | |||
OF3 | 1.76 | 0.22 | 1.91 | 0.13 | |||
OF4 | 1.88 | 0.23 | 2.99 | 0.94 | |||
OF5 | 2.82 | 0.25 | 3.03 | 0.37 | |||
OF6 | 2.15 | 0.21 | 1.99 | 0.26 | |||
OF7 | 2.77 | 0.32 | 1.92 | 0.41 | |||
OF8 | 2.52 | 0.25 | 2.44 | 0.15 | |||
OF9 | 2.06 | 0.26 | 1.79 | 0.27 | |||
OF10 | 2.17 | 0.27 | 2.25 | 0.45 | |||
OF11 | 2.39 | 0.22 | 2.23 | 0.3 | |||
OF12 | 3.11 | 0.31 | 2.91 | 0.47 | |||
OF13 | 1.44 | 0.25 | 1.49 | 0.42 | |||
OF14 | 1.93 | 0.29 | 1.8 | 0.19 | |||
OF15 | 1.86 | 0.27 | 1.67 | 0.25 |
Table 1. Average Weight Percent in the Upper Soil Layer by Dry Combustion and INS Methods.
It is interesting to compare the carbon distribution maps of the field based on the INS and DC methods (Figure 9 and 10). Both maps look very similar, but it should be noted that 2 days were spent on INS mapping, while ~2 months were required to process samples to create the DC map.
Figure 9. Carbon Distribution Map of the Camp Hill Field Based on the INS Method. Please click here to view a larger version of this figure.
Figure 10. Carbon Distribution Map of the Camp Hill Field Based on the DC Method. Please click here to view a larger version of this figure.
Field measurements of carbon content in scanning mode
Soil scientists are often interested in determining carbon content for large areas (e.g., 100 m x 100 m). Instead of determining carbon in locations 10 m apart (requiring 1 h per measurement using INS), it is possible to determine the average carbon content for a 100 m x 100 m field using the INS scanning mode. In scanning mode, it is possible to take INS measurements while passing over the entire field. This scanning measurement can be conducted in the same amount of time needed to measure a single location in static mode (1 h). The proof and principle of the INS scanning mode are demonstrated in this article.
It should be noted that the first attempt for measuring carbon in the scanning mode was less than satisfactory. The acquired scanning spectra were visibly different from the INS & TNC and TNC static mode spectra; the peaks of interest were broader and shorter with peak areas being much less than observed in the static mode. Investigations determined that this distortion was due to the influence of the Earth's magnetic field on the gamma detector's photomultiplier12. To resolve this problem, a magnetic screen (mu-metal) was used to shield the gamma-detector. Testing showed that the gamma spectra of a Co-60 control source was nearly identical regardless of the orientation of the screened gamma detector (vertical, horizontal, inclined), while peak centroids and peaks widths changed depending on the orientation of the unscreened detector. These result demonstrated that the effect of Earth's magnetic field on the photomultiplier can be suppressed by using a magnetic screen. Magnetic screening eliminated peak broadening and produced a scanning gamma spectra that looked very similar to the static mode spectra.
To compare static and scanning modes, static measurements of carbon content were performed (1 h each) at 5 random locations within a 15 m x 45 m field and measurement in scanning mode (1 h total) were performed on the same field which had a fairly uniform carbon content. A map of the field showing individual measurement locations and the scanning path is illustrated in Figure 11. The net INS spectra of the 5 static mode locations and that of the scanning mode are shown in Figure 12. As shown in the Figure 12, the scanning mode spectrum looks similar to the static mode spectra and falls in the mid-range of all static spectra.
Figure 11. Map of the Field Showing Static Measurement Locations (stars) and the Scanning Path (Lines). Please click here to view a larger version of this figure.
Figure 12. Net INS spectra for static and scanning modes; inset is a fragment of the net INS spectra around 4.44 MeV. Please click here to view a larger version of this figure.
Results of the net carbon peak area calculations are shown in Table 2. As can be seen from the presented data, the value of the net carbon peak area measured in scanning mode agrees with the average static mode value within the limits of experimental error. These results prove that INS scanning mode measurements can be used to define the average carbon content in a field. It is important to note that 5 h were spent determining average carbon content in static mode, whereas only 1 h was required in scanning mode.
Mode | Site # | Net Carbon | STD, | Field Average |
Peak Area, cps | cps | ±STD, cps | ||
Static | 1 | 64.8 | 3.9 | 63.3±3.8 |
2 | 58.1 | 3.5 | ||
3 | 65.4 | 3.4 | ||
4 | 68.9 | 4.1 | ||
5 | 59.4 | 4.1 | ||
Scanning | over field | 64.4 | 3.3 | 64.4±3.3 |
Table 2. Net Carbon Peak Area for Static and Scanning Modes.
Building on the foundation established by previous researchers, the NSDL staff addressed questions critical to the practical and successful use of this technology in real world field settings. Initially, NSDL researchers demonstrated the necessity to account for the INS system background signal when determining net carbon peak areas.11 Another effort showed that the net carbon peak area characterizes the average carbon weight percent in the upper 10 cm soil layer (regardless of carbon depth distribution shape) by direct proportional dependency. In addition, equipment required for INS system calibration (i.e., 1.5 m x 1.5 m x 0.6 m pits with different sand-carbon mixtures) was constructed and calibration procedures necessary for real world applications were developed and performed. The resultant calibration line renders it possible to determine soil carbon content from measured net carbon peak area. While NSDL researchers have incorporated many INS system design improvements, the recent addition of magnetic field shielding of gamma detectors allows for the practical use of the INS system scanning mode for large scale investigations of soil carbon.
Experience in applying the INS method for soil carbon analysis revealed several critical protocol steps. To obtain correct measurement results, it is critical to carefully check and adjust detector parameters using reference sources; this is very important for system stability and reproducing measurement results. The system background and calibration measurements are also critical steps for accurate determination of soil carbon content. Note that the detector parameters should be the same for both system background and calibration measurements. It is expedient to conduct calibration measurements (pits and system background) for several hours to increase accuracy of calibration coefficients. Installing magnetic screens on the detectors is critical for accurate measurement in the scanning mode since unscreened detectors produce very large errors due to the influence of the earth's magnetic field. In addition, magnetic screening improves results in the static mode.
The significance of using the INS method versus the "gold standard" DC method was demonstrated during field mapping. The speed of defining the carbon content by the INS method was ~30 times greater than the DC method. Other advantages of the INS method were discussed in the Introduction section.
Despite the demonstrated agreement between the INS and DC ("gold standard") methods, the current modification of the INS technique has one main limitation which is the minimal detectible level (1.5 w%). Since soil carbon content can be less than this, future efforts will concentrate on improving the sensitivity of the INS system by increasing the number of gamma detectors and optimizing the overall system design or by applying target neutron methods.13
Despite this limitation, the current modification of the INS system can be recommended for soil carbon determination of individual locations and for carbon distribution mapping of field terrains. Possible future work using the INS method may explore measuring other soil elements such as nitrogen, iron, and hydrogen.
The authors have nothing to disclose.
The authors are indebted to Barry G. Dorman, Robert A. Icenogle, Juan Rodriguez, Morris G. Welch, and Marlin Siegford for technical assistance in experimental measurements, and to Jim Clark and Dexter LaGrand for assistance with computer simulations. We thank XIA LLC for allowing the use of their electronics and detectors in this project. This work was supported by NIFA ALA Research Contract No ALA061-4-15014 “Precision geospatial mapping of soil carbon content for agricultural productivity and lifecycle management”.
Neutron Generator | Thermo Fisher Scientific, Colorado Springs, CO DNC software |
MP320 | |
Gamma-detector: | na | ||
– NaI(Tl) crystal | Scionix USA, Orlando, FL | ||
– Electronics | XIA LLC, Hayward, CA | ||
– Software | ProSpect | ||
Battery | Fullriver Battery USA, Camarillo, CA | DC105-12 | |
Invertor | Nova Electric, Bergenfield, NJ | CGL 600W-series | |
Charger | PRO Charging Systems, LLC, LaVergne, TN | PS4 | |
Block of Iron | Any | na | |
Boric Acid | Any | na | |
Laptop | Any | na | |
mu-metal | Magnetic Shield Corp., Bensenville, IL | MU010-12 | |
Construction sand | Any | na | |
Coconut shell | General Carbon Corp., Patterson, NJ | GC 8 X 30S | |
Reference Cs-137 source | Any | na |