The field study site is on a farm managed by David Drake and we thank him for his cooperation during this long-term research. We also thank Mike Murphy for his many years of field work on this project and Robin Montenieri for her help with graphics used in this paper.

1. Land Surface Elevation Data Collection
<ol>
	<li>GPS calibration for site
	<ol>
		<li>Locate or set a stable benchmark in a secure location at the survey site for use as the base station GPS for RTKGPS data collection.</li>
		<li>Set up base station for RTKGPS data collection at this local benchmark using best approximation of coordinates for the base station location (i.e., WAAS-corrected GPS position).</li>
		<li>With the rover GPS, visit at least three horizontal and vertical control point benchmarks within the radio communication limits of the RTKGPS (approximately 10 km radius) and record positions. 
		NOTE: Benchmarks described by the National Geodetic Survey can be searched online33 and were used here.</li>
		<li>Given the measured and published coordinates of the control points, use RTKGPS field software to perform a site calibration34, solving for coordinates of local benchmark to be used as base station. Check that coordinate residuals (horizontal and vertical) for the control points are within tolerable limits (&#177; 0.02 m for this calibration).</li>
	</ol>
	</li>
	<li>GPS point data collection
	<ol>
		<li>With base station GPS set on local benchmark and using local site calibration and field software, record RTKGPS position data into GPS data collector at approximately a 5-m horizontal spacing throughout the survey area.
		<ol>
			<li>Collect data efficiently by mounting rover GPS antenna at a measured fixed height above the ground surface on a vehicle and driving transects through the area (Figure 6).</li>
		</ol>
		</li>
		<li>For vehicle method, define transect endpoints to create parallel transects spaced 5 m apart. Import transect endpoints into GPS data collector for navigation of transects while driving. Collect points with the data collector automatically once per second while driving transects at approximately 5 m s-1 to obtain point data about every 5 m.</li>
		<li>Repeat point data collection at site as described above at later time (8 to 9 years later in this study) so that land surface elevation changes can be analyzed; the original GPS site calibration is used for all surveys and is not repeated.</li>
	</ol>
	</li>
</ol>
<img alt="Figure 6" class="xfigimg" src="/files/ftp_upload/56064/56064fig6.jpg" /> 
Figure 6. RTKGPS Elevation Surface Data Collection. RTKGPS elevation surface data is collected while driving a utility vehicle through the field (a), while real-time GPS corrections are provided by the on-site base station (b).
2. DEM Creation and Processing
<ol>
	<li>Creating the DEMs
	<ol>
		<li>Import position data into GIS software and interpolate to a 5-m grid DEM. Using GIS software, cross-validate measured point elevations to interpolated elevation values and choose an interpolation method that minimizes these cross-validation errors. 
		NOTE: Ordinary kriging with a Gaussian semivariogram model was the optimal interpolation method for the elevation data at this site. Cross-validation also provides a measure of elevation accuracy for the survey method35.</li>
		<li>Repeat 2.1.1 for second set of position data to create second DEM.</li>
	</ol>
	</li>
	<li>Mapping DEM change
	<ol>
		<li>Using a raster calculator tool in GIS, subtract the most recent DEM from the original DEM to create a raster map of DEM change (Figure 5), where negative values of elevation change represent erosion and positive values represent deposition.</li>
	</ol>
	</li>
	<li>Land Classification
	<ol>
		<li>Compute land surface topographic attributes (slope, aspect, contributing area) from the first grid DEM using DEM processing software.</li>
		<li>Classify land areas as summit, sideslope, or toeslope based on slope and contributing area of each DEM grid cell. 
		NOTE: Summits are represented by low slopes and low contributing areas. Sideslopes are represented by high slopes and intermediate contributing areas. Toeslopes are represented by low slopes and high contributing areas. Slope and contributing area values defining these classifications will depend on the land surface topography at the site and are qualitatively chosen to give the desired representation of each classification area for a particular site.</li>
		<li>Divide sideslope areas by the two dominant aspects, north-facing and south-facing at this site.</li>
	</ol>
	</li>
</ol>
3. Soil Sampling
<ol>
	<li>Sample Planning
	<ol>
		<li>Reference maps in GIS to plan soil sample locations. Choose a number of locations to adequately represent all landscape positions.</li>
		<li>Upload sample location coordinates to the GPS data collector so that sample sites can be located in the field.</li>
		<li>Use prior knowledge of soils at the site to guide decisions of sample depth increments in order to capture CaCO3 variability. Pre-label sealable plastic bags to indicate sample location and depth increment.</li>
	</ol>
	</li>
	<li>Field Sampling
	<ol>
		<li>Drive to the sample sites with a utility vehicle equipped with a hydraulic soil coring machine and the RTKGPS rover antenna for navigation.</li>
		<li>Using soil coring machine and sampling tube for desired soil core diameter (5.1 cm in this study), extract soil core from each sample location (Figure 7). 
		NOTE: The number of cores extracted at each location, as well as the soil core depth and increments varied in this study. In 2001, a single core to a depth of 90 cm was taken and divided into 30 cm increments. In 2012, two soil cores were taken (within 1 m of the corresponding 2001 sample) to a depth of 30 cm and divided into 15 cm increments, with the two cores being aggregated for analyses. The 2012 method is recommended.</li>
		<li>Record RTKGPS position data (x, y, z) at each sample location.</li>
		<li>Cut the soil core into desired depth increments and transfer into pre-labelled sealable plastic bags and then place in coolers for transport back to the laboratory.</li>
		<li>Repeat field sampling after significant erosion and/or deposition has occurred (11 years between samples in this study).</li>
	</ol>
	</li>
</ol>
<img alt="Figure 7" class="xfigimg" src="/files/ftp_upload/56064/56064fig7.jpg" /> 
Figure 7. Soil Sampling. Soil sample locations are navigated to using a GPS-guided utility vehicle equipped with a hydraulic soil coring machine (a) so that soil cores can be extracted (b) and divided into desired depth increments.
<ol start="3">
	<li>Position Data Processing
	<ol>
		<li>Measure differences in elevations recorded at each soil sample location between the two sample dates (198 locations sampled in 2001 and 2012 in this study). 
		NOTE: Elevations for 2001 were taken from the 2001 DEM since point elevations were not recorded at the time of soil sampling. Positive changes in elevation &#62; 0.05 m are considered depositional sites, while negative changes in elevation &#60; -0.05 m are considered erosional sites.</li>
		<li>Classify each sample location as summit, North-facing sideslope, South-facing sideslope, or toeslope based on the DEM processing (see Protocol 2.3.2); classification at a single location, as defined by the slope and contributing area criteria, may be reclassified to match the dominant classification of surrounding points.</li>
		<li>Use spatial joining tools in GIS software to assign sample locations to other spatial data layers used for analyses (management strip and soil mapping unit).</li>
	</ol>
	</li>
</ol>
4. Soils Analyses
<ol>
	<li>Soil Sample Preparation
	<ol>
		<li>Dry soil samples from the field at 60 &#176;C in a laboratory oven overnight.</li>
		<li>Grind oven-dried soils to pass through a 2 mm sieve using a motorized grinder or a mortar and pestle.</li>
	</ol>
	</li>
	<li>Modified Pressure-Calcimeter Apparatus Setup
	<ol>
		<li>Set up the modified pressure-calcimeter apparatus (Figure 8) by connecting a pressure transducer (0 - 105 kPa range, 0.03 - 5 V DC output) to a power supply with 14 gauge wire and a digital voltmeter wired in line to monitor output from the transducer.
		<ol>
			<li>Attach 9.5-mm ID tubing to the base of the pressure transducer and connect tubing to an 18 gauge Luer lock hypodermic needle with a particle filter (0.6 &#181;m) in the middle to collect any reflux from reaching the pressure transducer.</li>
		</ol>
		</li>
		<li>Use serum bottles as reaction vessels connected to the pressure transducer (Figure 9).&#160; Determine the size of the serum bottle to use by wetting a metal tablespoon with water and adding approximately 5 mL of soil that you expect to have high CaCO3 concentration.&#160; Pipet 1 mL of 0.5 N H2SO4 to this soil and observe effervescence.</li>
		<li>If effervescence is high, then assume greater than 15% CaCO3 concentration and use a 100 mL serum bottle as the reaction vessel, otherwise use a 20 mL serum bottle.</li>
	</ol>
	</li>
</ol>
<img alt="Figure 8" class="xfigimg" src="/files/ftp_upload/56064/56064fig8.jpg" /> 
Figure 8. Modified Pressure-calcimeter Apparatus. The modified pressure-calcimeter apparatus uses a serum bottle as the reaction vessel and a pressure transducer wired to a voltage meter to output the signal (from Sherrod et al.)30.
<img alt="Figure 9" class="xfigimg" src="/files/ftp_upload/56064/56064fig9.jpg" /> 
Figure 9. Reaction Vessels for the Modified Pressure-calcimeter Method. Reaction vessels for the modified pressure-calcimeter method are serum bottles containing a 0.5 dram vial with 2 mL acid reagent and a 1 g soil sample.
<ol start="3">
	<li>Carbonate Measurement
	<ol>
		<li>Place a 1 g subsample of the prepared soil (see Protocol 4.1) into a labelled reaction vessel. For soils containing greater than 50% CaCO3, use only 0.5 g of soil.</li>
		<li>Pipet 2 mL of acid reagent (6 N HCI containing 3% FeCL2O4H2O) into a 0.5 g glass vial. Place vial into reaction vessel gently so that the solution contents do not spill out by tilting the reaction vessel almost to the side position.</li>
		<li>While keeping the reaction vessel containing the soil sample and acid vial tilted, seal with gray butyl rubber stoppers and crimp with aluminum sealing ring.</li>
		<li>Shake reaction vessel with a swirling motion to insure complete mixing of soil with acid. Place the reaction vessel on the lab bench and let the reaction proceed for at least 2 hr.</li>
		<li>While waiting for reaction vessels to complete, determine a standard curve by measuring voltages of known CaCO3 concentrations using the same reaction vessel setup as soil samples (Figure 10).&#160; Mix 100% CaCO3 with glass beads or sand on a weight percentage basis to create known CaCO3 concentrations.&#160; Include a blank sample without CaCO3.</li>
		<li>After soil sample reactions are complete, pierce the rubber septum of the reaction vessel with an 18 gauge hypodermic needle and record voltage output by the pressure transducer.</li>
		<li>Solve for CaCO3 percentage given the measured voltage and the equation determined from the standard curve (Figure 10a). 
		NOTE: Increase in pressure produced by release of CO2 is linearly related to the concentration of CaCO3 present in the soil such that: % CaCO3 = (regression coefficient * change in pressure in volts) + intercept.</li>
	</ol>
	</li>
</ol>
<img alt="Figure 10" class="xfigimg" src="/files/ftp_upload/56064/56064fig10.jpg" /> 
Figure 10. CaCO3 Measurement. (a) A standard curve for CaCO3 is created using the voltage readings from the pressure transducer based on known percentages of CaCO3 (b) mixed with powder glass beads or sand.
5. Statistical Analyses
<ol>
	<li>Define two dependent variables as the change in ground surface elevations and soil surface layer CaCO3 concentrations from the first to second sample dates (2001 to 2012 in this study). Define independent or explanatory variables as management (odd- or even- numbered strip), individual strips, west or east block of strips, soil mapping unit, landscape classification, and erosional/depositional classification.</li>
	<li>Perform correlation analysis and analysis of variance to statistically quantify relationships between variables. Perform analyses in any preferred statistical package.</li>
</ol>

<ol>
	<li>Freebairn, D. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Erosion+control+-+some+observations+on+the+role+of+soil+conservation+structures+and+conservation.">Erosion control - some observations on the role of soil conservation structures and conservation.</a> Nat. Res. Mgt. 7 (1), 8-13 (2004).</li><li>Garcia-Orenes, F., Roldan, A., Mataix-Solera, J., Cerda, A., Campoy, M., Arcenegui, V., Caravaca, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Soil+structural+stability+and+erosion+rates+influenced+by+agricultural+management+practices+in+a+semi-arid+Mediterranean+agro-ecosystem.">Soil structural stability and erosion rates influenced by agricultural management practices in a semi-arid Mediterranean agro-ecosystem.</a> Soil Use and Mgt. 28, 571-579 (2012).</li><li>Hass, H. J., Willis, W. O., Bond, J. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=General+relationships+and+conclusions.">General relationships and conclusions.</a> Summer Fallow in the Western United States. USDA-ARS Conserv. Res. Rpt. No. 17. , 149-160 (1974).</li><li>Montgomery, D. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Soil+erosion+and+agricultural+sustainability.">Soil erosion and agricultural sustainability.</a> Proc. of the Nat. Acad. of Sci. of the USA. 104 (33), 13268-13272 (2007).</li><li>Skidmore, E. L., Layton, J. B., Armbrust, D. V., Hooker, M. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Soil+physical+properties+as+influenced+by+cropping+and+residue+management.">Soil physical properties as influenced by cropping and residue management.</a> Soil Sci. Soc. of Am. J. 50 (2), 415-419 (1986).</li><li>Sherrod, L. A., Erskine, R. H., Green, T. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Spatial+patterns+and+cross-correlations+of+temporal+changes+in+soil+carbonates+and+surface+elevation+in+a+winter+wheat-fallow+cropping+system.">Spatial patterns and cross-correlations of temporal changes in soil carbonates and surface elevation in a winter wheat-fallow cropping system.</a> Soil Sci. Soc. of Am. J. 79 (2), 417-427 (2015).</li><li>Stroosnijder, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Measurement+of+erosion:+Is+it+possible?.">Measurement of erosion: Is it possible?.</a> Catena. 64 (2-3), 162-173 (2005).</li><li>D&#261;bek, P., &#379;muda, R., &#262;mielewski, B., Szczepa&#324;ski, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Analysis+of+water+erosion+processes+using+terrestrial+laser+scanning.">Analysis of water erosion processes using terrestrial laser scanning.</a> Acta Geodynam. Et Geomat. 11 (1), 45-52 (2014).</li><li>Day, S. S., Gran, K. B., Belmont, P., Wawrzyniec, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Measuring+bluff+erosion+part+1:+terrestrial+laser+scanning+methods+for+change+detection.">Measuring bluff erosion part 1: terrestrial laser scanning methods for change detection.</a> Earth Surf. Proc. and Landforms. 38 (10), 1055-1067 (2013).</li><li>Eltner, A., Baumgart, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Accuracy+constraints+of+terrestrial+Lidar+data+for+soil+erosion+measurement:+Application+to+a+Mediterranean+field+plot.">Accuracy constraints of terrestrial Lidar data for soil erosion measurement: Application to a Mediterranean field plot.</a> Geomorph. 245, 243-254 (2015).</li><li>Letortu, P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Retreat+rates,+modalities+and+agents+responsible+for+erosion+along+the+coastal+chalk+cliffs+of+Upper+Normandy:+The+contribution+of+terrestrial+laser+scanning.">Retreat rates, modalities and agents responsible for erosion along the coastal chalk cliffs of Upper Normandy: The contribution of terrestrial laser scanning.</a> Geomorph. 245, 3-14 (2015).</li><li>Longoni, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Monitoring+Riverbank+Erosion+in+Mountain+Catchments+Using+Terrestrial+Laser+Scanning.">Monitoring Riverbank Erosion in Mountain Catchments Using Terrestrial Laser Scanning.</a> Rem. Sens. 8 (3), 241 (2016).</li><li>Meijer, A. D., Heitman, J. L., White, J. G., Austin, R. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Measuring+erosion+in+long-term+tillage+plots+using+ground-based+lidar.">Measuring erosion in long-term tillage plots using ground-based lidar.</a> Soil &amp; Till. Res. 126, 1-10 (2013).</li><li>Rengers, F. K., Tucker, G. E., Moody, J. A., Ebel, B. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Illuminating+wildfire+erosion+and+deposition+patterns+with+repeat+terrestrial+lidar.">Illuminating wildfire erosion and deposition patterns with repeat terrestrial lidar.</a> J. of Geophys. Res.-Earth Surf. 121 (3), 588-608 (2016).</li><li>Schubert, J. E., Gallien, T. W., Majd, M. S., Sanders, B. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Terrestrial+Laser+Scanning+of+Anthropogenic+Beach+Berm+Erosion+and+Overtopping.">Terrestrial Laser Scanning of Anthropogenic Beach Berm Erosion and Overtopping.</a> J. of Coast. Res. 31 (1), 47-60 (2015).</li><li>Stenberg, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Evaluation+of+erosion+and+surface+roughness+in+peatland+forest+ditches+using+pin+meter+measurements+and+terrestrial+laser+scanning.">Evaluation of erosion and surface roughness in peatland forest ditches using pin meter measurements and terrestrial laser scanning.</a> Earth Surf. Proc. and Landforms. 41 (10), 1299-1311 (2016).</li><li>Croke, J., Todd, P., Thompson, C., Watson, F., Denham, R., Khanal, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+use+of+multi+temporal+LiDAR+to+assess+basin-scale+erosion+and+deposition+following+the+catastrophic+January+2011+Lockyer+flood,+SE+Queensland,+Australia.">The use of multi temporal LiDAR to assess basin-scale erosion and deposition following the catastrophic January 2011 Lockyer flood, SE Queensland, Australia.</a> Geomorph. 184, 111-126 (2013).</li><li>Earlie, C., Masselink, G., Russell, P., Shail, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Sensitivity+analysis+of+the+methodology+for+quantifying+cliff+erosion+using+airborne+LiDAR+-+examples+from+Cornwall,+UK.">Sensitivity analysis of the methodology for quantifying cliff erosion using airborne LiDAR - examples from Cornwall, UK.</a> J. of Coast. Res. Spec. Iss. 65, 470-475 (2013).</li><li>Kessler, A. C., Gupta, S. C., Dolliver, H. A. S., Thoma, D. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Lidar+Quantification+of+Bank+Erosion+in+Blue+Earth+County,+Minnesota.">Lidar Quantification of Bank Erosion in Blue Earth County, Minnesota.</a> J. of Env. Quality. 41 (1), 197-207 (2012).</li><li>Pye, K., Blott, S. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Assessment+of+beach+and+dune+erosion+and+accretion+using+LiDAR:+Impact+of+the+stormy+2013-14+winter+and+longer+term+trends+on+the+Sefton+Coast,+UK.">Assessment of beach and dune erosion and accretion using LiDAR: Impact of the stormy 2013-14 winter and longer term trends on the Sefton Coast, UK.</a> Geomorph. 266, 146-167 (2016).</li><li>Thoma, D. P., Gupta, S. C., Bauer, M. E., Kirchoff, C. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Airborne+laser+scanning+for+riverbank+erosion+assessment.">Airborne laser scanning for riverbank erosion assessment.</a> Rem. Sens. of Env. 95 (4), 493-501 (2005).</li><li>Zhang, C. L., Yang, S., Pan, X. H., Zhang, J. Q. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Estimation+of+farmland+soil+wind+erosion+using+RTK+GPS+measurements+and+the+Cs-137+technique:+A+case+study+in+Kangbao+County,+Hebei+province,+northern+China.">Estimation of farmland soil wind erosion using RTK GPS measurements and the Cs-137 technique: A case study in Kangbao County, Hebei province, northern China.</a> Soil &amp; Till. Res. 112 (2), 140-148 (2011).</li><li>Neugirg, F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Erosion+processes+in+calanchi+in+the+Upper+Orcia+Valley,+Southern+Tuscany,+Italy+based+on+multitemporal+high-resolution+terrestrial+LiDAR+and+UAV+surveys.">Erosion processes in calanchi in the Upper Orcia Valley, Southern Tuscany, Italy based on multitemporal high-resolution terrestrial LiDAR and UAV surveys.</a> Geomorph. 269, 8-22 (2016).</li><li>Pineux, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Can+DEM+time+series+produced+by+UAV+be+used+to+quantify+diffuse+erosion+in+an+agricultural+watershed?.">Can DEM time series produced by UAV be used to quantify diffuse erosion in an agricultural watershed?.</a> Geomorph. 280, 122-136 (2017).</li><li>Bremer, M., Sass, O. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Combining+airborne+and+terrestrial+laser+scanning+for+quantifying+erosion+and+deposition+by+a+debris+flow+event.">Combining airborne and terrestrial laser scanning for quantifying erosion and deposition by a debris flow event.</a> Geomorph. 138 (1), 49-60 (2012).</li><li>Day, S. S., Gran, K. B., Belmont, P., Wawrzyniec, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Measuring+bluff+erosion+part+2:+pairing+aerial+photographs+and+terrestrial+laser+scanning+to+create+a+watershed+scale+sediment+budget.">Measuring bluff erosion part 2: pairing aerial photographs and terrestrial laser scanning to create a watershed scale sediment budget.</a> Earth Surf. Proc. and Landforms. 38 (10), 1068-1082 (2013).</li><li>De Rose, R. C., Basher, L. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Measurement+of+river+bank+and+cliff+erosion+from+sequential+LIDAR+and+historical+aerial+photography.">Measurement of river bank and cliff erosion from sequential LIDAR and historical aerial photography.</a> Geomorph. 126 (1-2), 132-147 (2011).</li><li>Perroy, R. L., Bookhagen, B., Asner, G. P., Chadwick, O. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Comparison+of+gully+erosion+estimates+using+airborne+and+ground-based+LiDAR+on+Santa+Cruz+Island,+California.">Comparison of gully erosion estimates using airborne and ground-based LiDAR on Santa Cruz Island, California.</a> Geomorph. 118 (3-4), 288-300 (2010).</li><li>Loeppert, R. H., Suarez, D. L., Sparks, D. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Carbonate+and+Gypsum.">Carbonate and Gypsum.</a> Methods of Soil Analysis. Part 3. Chemical Methods. , 437-474 (1996).</li><li>Sherrod, L. A., Dunn, G., Peterson, G. A., Kilberg, R. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Inorganic+carbon+analysis+by+modified+pressure-calcimeter+method.">Inorganic carbon analysis by modified pressure-calcimeter method.</a> Soil Sci. Soc. of Am. J. 66 (1), 299-305 (2002).</li><li>McCutcheon, M. C., Farahani, H. J., Stednick, J. D., Buchleiter, G. W., Green, T. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Effect+of+soil+water+on+apparent+soil+electrical+conductivity+and+texture+relationships+in+a+dryland+field.">Effect of soil water on apparent soil electrical conductivity and texture relationships in a dryland field.</a> Biosyst. Eng. 94 (1), 19-32 (2006).</li><li>Wheaton, J. M., Brasington, J., Darby, S. E., Sear, D. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Accounting+for+uncertainty+in+DEMs+from+repeat+topographic+surveys:+improved+sediment+budgets.">Accounting for uncertainty in DEMs from repeat topographic surveys: improved sediment budgets.</a> Earth Surf. Proc. and Landforms. 35 (2), 136-156 (2010).</li><li>. Survey Marks and Datasheets Available from: <a target="_blank" href="https://www.ngs.noaa.gov/datasheets/">https://www.ngs.noaa.gov/datasheets/</a> (2017)</li><li>Trimble Inc. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Trimble Access Software – General Survey. Version 1.60. Revision A. , (2011).</li><li>Erskine, R. H., Green, T. R., Ramirez, J. A., MacDonald, L. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Digital+elevation+accuracy+and+grid+cell+size:+effects+on+estimated+terrain+attributes.">Digital elevation accuracy and grid cell size: effects on estimated terrain attributes.</a> Soil Sci. Soc. of Am. J. 71, 1371-1380 (2007).</li></ol>

The authors have nothing to disclose.

Mapped changes in elevation (Figure 5) illustrate significant erosion and deposition on an agricultural field and spatial patterns indicative of multiple controlling factors over multiple scales. From field scale patterns associated with wind, down to fine scale dendritic patterns produced by water flow, processes relevant to this study are discernable. The level of elevation change detection provided by repeated RTKGPS ground surveys appears optimal. Finer detection levels, as provided by T...

Soil erosion threatens the sustainability of agricultural lands. Crop management, such as a conventionally-tilled winter wheat-fallow crop rotation, can accelerate erosion and deposition processes as bare soils during fallow periods are more susceptible to wind and water forces1,2,3,4,5 (Figure 1). While these processes might be evident, they can be difficult to quantify.
The purpose of this study is first to provide an efficient method for quantifying and describing spatial patterns of erosion and deposition at the field scale using global positioning system (GPS) technology and geographic information systems (GIS) mapping tools. A simple conceptual model relating these patterns to near-surface soil carbonates (CaCO3) is also presented and tested by prescribed field and laboratory methods. These relationships provide indirect measures of erosion and deposition, while validating the results of the GPS method. The present paper emphasizes the methods used in Sherrod et al. so that they can be repeated, in part or whole, for similar research in other locations6.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/56064/56064fig1.jpg" /> 
Figure 1. Photos of (a) Erosion and (b) Deposition at the Study Site Following a Heavy Rainfall Event. A tractor tire track in the lower right corner of photo (b) indicates the depth of deposition at the wheat/fallow strip border.
Various direct methods for measuring soil erosion were reviewed by Stroosnijder7. Suggested methods vary with the measurement purpose and resources available, but a &#34;change in surface elevation&#34; method is recommended at the hillslope scale and provides the advantage of measuring both erosion and deposition. One way to apply this method is to install pins in the soil and monitor the change in height of the soil relative to the top of the pin7. With advances in land surveying technology, however, this labor-intensive approach can be replaced by other techniques, such as terrestrial laser scanning (TLS)8,9,10,11,12,13,14,15,16, airborne laser scanning (ALS)17,18,19,20,21, GPS6,22, advanced photogrammetry23,24, or combinations of these techniques25,26,27. While laser scanning, commonly referred to as LiDAR (Light Detection And Ranging), provides the most rapid acquisition of dense surface elevation data sets, corrections must be made to remove standing objects, such as vegetation. With millimeter-level vertical precision, TLS can detect the smallest elevation change, however Perroy et al. recommended ALS over TLS for gulley erosion estimates due to the larger scanning footprint and better instrument orientation (less topographic shadowing) for scanning into deeply incised gullies28. Real-time kinematic GPS (RTKGPS), providing centimeter-level precision without data post-processing, is used for this study. The spatial resolution and precision of RTKGPS-collected data are optimal for detecting the dominant erosional and depositional features in an agricultural field or other environments with substantial ground cover.
The pressure-calcimeter method for quantifying soil CaCO3 relies on the soil&#39;s reaction to acid in a closed system, resulting in the release of CO2. The increase in pressure within the reaction vessel at a constant temperature is linearly correlated to the amount of soil CaCO329. Modifications to the traditional pressure-calcimeter method, described by Sherrod et al., include changing the reaction vessel to serum bottles and using a pressure transducer wired to a digital voltmeter for the detection of pressure changes30. These modifications allow for lower detection limits and a higher capacity for daily soil sample runs. Gravimetric or simple titrimetric methods for soil CaCO3 measurement produced larger errors and detection limits than this modified pressure-calcimeter method30.
Conceptual Model
When direct measures of erosion and deposition are not feasible, indirect indicators of these processes may be used. Sherrod et al. hypothesized that soil surface layer CaCO3 concentration in a semi-arid climate is inversely correlated with the change in ground surface elevation (positively correlated with erosion, negatively correlated with deposition)6. The hypothesis should apply broadly, but specific relationships will depend upon site conditions (soil, vegetation, management, and climate). Soils at the test site (Table 1) typically contain a distinct calcareous layer 15-20 cm below the soil surface. Conceptually, erosion will remove the surface layer of relatively low CaCO3 concentration leaving this calcareous layer of high CaCO3 closer to the soil surface. The low CaCO3 soil is then transported to the depositional areas, causing the calcareous layer to be buried deeper below the soil surface (Figure 2). Sampling these soils over time at appropriate depth intervals, either erosion or deposition (or neither) may be inferred by CaCO3 concentration, according to this model.
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>Soil Series</td>
			<td>Slope</td>
			<td>Taxonomic classification</td>
			<td>Depth</td>
			<td>pH</td>
			<td>EC</td>
			<td>Total N</td>
			<td>SOC</td>
			<td>CaCO3</td>
		</tr>
		<tr>
			<td></td>
			<td>%</td>
			<td></td>
			<td>cm</td>
			<td>1:2</td>
			<td>dS m-1</td>
			<td>g kg-1</td>
			<td>g kg-1</td>
			<td>g kg-1</td>
		</tr>
		<tr>
			<td>Colby loam</td>
			<td>5-9</td>
			<td>fine-silty, mixed, superactive, calcareous, mesic Aridic Ustorthent</td>
			<td>0-15</td>
			<td>8.2</td>
			<td>0.24</td>
			<td>0.7</td>
			<td>6.1</td>
			<td>69.8</td>
		</tr>
		<tr>
			<td></td>
			<td></td>
			<td></td>
			<td>15-30</td>
			<td>8.3</td>
			<td>0.24</td>
			<td>0.5</td>
			<td>4.0</td>
			<td>84.3</td>
		</tr>
		<tr>
			<td>Kim sandy loam</td>
			<td>2-5</td>
			<td>fine-loamy, mixed, active, calcareous, mesic Ustic Torriorthent</td>
			<td>0-15</td>
			<td>7.8</td>
			<td>0.26</td>
			<td>0.8</td>
			<td>7.0</td>
			<td>29.8</td>
		</tr>
		<tr>
			<td></td>
			<td></td>
			<td></td>
			<td>15-30</td>
			<td>8.0</td>
			<td>0.27</td>
			<td>0.6</td>
			<td>5.0</td>
			<td>51.5</td>
		</tr>
		<tr>
			<td></td>
			<td>5-9</td>
			<td>fine-loamy, mixed, active, calcareous, mesic Ustic Torriorthent</td>
			<td>0-15</td>
			<td>8.1</td>
			<td>0.22</td>
			<td>0.6</td>
			<td>5.4</td>
			<td>26.7</td>
		</tr>
		<tr>
			<td></td>
			<td></td>
			<td></td>
			<td>15-30</td>
			<td>8.1</td>
			<td>0.19</td>
			<td>0.5</td>
			<td>4.1</td>
			<td>25.8</td>
		</tr>
		<tr>
			<td>Wagonwheel loam</td>
			<td>0-2</td>
			<td>coarse-silty, mixed, superactive, mesic Aridic Calciustept</td>
			<td>0-15</td>
			<td>8.2</td>
			<td>0.23</td>
			<td>0.7</td>
			<td>5.9</td>
			<td>66.2</td>
		</tr>
		<tr>
			<td></td>
			<td></td>
			<td></td>
			<td>15-30</td>
			<td>8.2</td>
			<td>0.23</td>
			<td>0.6</td>
			<td>3.7</td>
			<td>98.1</td>
		</tr>
		<tr>
			<td></td>
			<td>2-5</td>
			<td>coarse-silty, mixed, superactive, mesic Aridic Calciustept</td>
			<td>0-15</td>
			<td>8.3</td>
			<td>0.23</td>
			<td>0.8</td>
			<td>6.6</td>
			<td>52.0</td>
		</tr>
		<tr>
			<td></td>
			<td></td>
			<td></td>
			<td>15-30</td>
			<td>8.4</td>
			<td>0.26</td>
			<td>0.7</td>
			<td>5.4</td>
			<td>118.3</td>
		</tr>
	</tbody>
</table>
Table 1. Soils at the Test Site. Soil mapping units and taxonomic classification, with average soil pH, electrical conductivity (EC), total N, soil organic C (SOC), and CaCO3 concentrations in the 0- to 15- and 15- to 30-cm depth increments for the Scott field in 2012 (from Sherrod et al.)6.
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/56064/56064fig2.jpg" /> 
Figure 2. Conceptual Soil Profiles. Conceptual soil profiles for (a) a static soil matrix with CaCO3 leached from the surface layer and precipitated in a deeper layer, (b) moderate erosion of the surface layer, and (c) moderate deposition of material above the previous surface layer. Depth intervals (left) are approximate based on site data (from Sherrod et al.)6. <a href="//cloudflare.jove.com/files/ftp_upload/56064/56064fig2large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
Site Description and History
The 109-ha Scott field is part of the Drake Farm in northeastern Colorado (40.61oN, 104.84oW, Figure 3) and was monitored from 2001 to 2012 for this study. Average annual precipitation and evapotranspiration were approximately 350 and 1200 mm, respectively, in this semi-arid climate, where convective rain of short duration and high intensity were common during the summer. Elevations range from 1559 to 1588 m in this undulating terrain with distinct landscape positions: summit, sideslope north-facing (side-NF), sideslope south-facing (side-SF), and toeslope (Figure 4b). Alternating strips (~120 m wide) were typically managed in this rainfed winter wheat-fallow rotation such that every other strip was fallow for about 14 months out of every 24-month rotation cycle. Shallow tillage (~7 cm), typically v-blade sweeps, occurred 4 to 6 times through the fallow period for weed control. Soils at the site were classified to have a soil-loss tolerance, or T value, of 11 Mg ha-1 year-1, where erosion rates below this T value are considered acceptable for continued agricultural production4.
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/56064/56064fig3.jpg" /> 
Figure 3. Site Location is Shown on a Topographic Relief Image (1011 to 4401 m) of the State of Colorado, USA. Mean elevation of the site is 1577 m.
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/56064/56064fig4.jpg" /> 
Figure 4. Soils map and Land Surface Elevation of the Scott Field. (a) Soils map of the Scott Field showing point soil sample locations and the crop management strips. Soil unit abbreviations are: 1 = Wagonwheel loam 0-2% slope, 2 = Wagonwheel loam 2-5% slope, 3 = Colby loam 5 - 9% slope, 4 = Kim fine sandy loam 2 - 5% slope, 5 = Kim fine sandy loam 5-9% slope; and (b) land surface elevation of the field based on the 2001 5-m grid digital elevation model (DEM), with soil sample locations shown by land classification (from Sherrod et al.)6.
The first ground surface elevation survey was collected by RTKGPS in 2001 to produce a digital elevation model (DEM) for the site. In conjunction with McCutcheon et al., an intensive soil sample (Figure 4a) was also performed in 2001, from which surface soil CaCO3 were analyzed by a modified pressure-calcimeter method30,31. Visually evident erosion and deposition occurring over the subsequent decade due to wind, predominantly from the northwest, and rainfall-runoff events prompted a second RTKGPS elevation survey in 2009 (with a portion of the field completed in 2010). Comparison of the new DEM to the original 2001 DEM via a DEM of Difference map32 confirmed significant erosion and deposition, displaying patterns which suggested multiple controlling factors for these processes (Figure 5). Given the substantial surface soil redistribution at the site and the historical soil CaCO3 data, the 2001 soil sample was repeated in 2012 to test a conceptual model of hydropedological processes6, as described in the previous section.
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/56064/56064fig5.jpg" /> 
Figure 5. Map of Changes (2001-2009*) in Land Surface elevation (&#916;z) on a 5-m Grid within the Scott Field in Northeastern Colorado. Crop strip numbers are labeled over the alternating winter-wheat-fallow cropping system, and section A-A&#39; is shown (details given in Figure 11). *Strips 2, 4, 6, 8 surveyed in 2010 to complete the 2009 DEM (from Sherrod et al.)6. <a href="//cloudflare.jove.com/files/ftp_upload/56064/56064fig5large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>

<table><tbody><tr><td>Real-time kinematic GPS system</td><td>Trimble</td><td>Model 5800</td><td></td></tr><tr><td>GPS field data collector</td><td>Trimble</td><td>Model TSC2</td><td></td></tr><tr><td>GPS field software</td><td>Trimble</td><td></td><td>Trimble Access (Trimble Survey Controller used in 2001 for site calibration but this software is no longer supported)</td></tr><tr><td>Hydraulic soil coring machine</td><td>Giddings Machine Company</td><td></td><td></td></tr><tr><td>Utility vehicle</td><td>John Deere</td><td>Gator 6x4</td><td></td></tr><tr><td>GIS software</td><td>ESRI</td><td></td><td>ArcGIS for Desktop with Spatial Analyst and Geostatistical Analyst Extensions</td></tr><tr><td>Statistical software</td><td>SAS</td><td></td><td>SAS Institute Inc.</td></tr><tr><td>Pressure transducer 0-105 kPa</td><td>Serta</td><td>Model 280E</td><td>Setra Systems, In., Boxborough, MA</td></tr><tr><td>Volt meter</td><td>WaveTek</td><td>5XL</td><td>Digital meter set to read volts</td></tr><tr><td>Serum Bottles</td><td>Wheaton</td><td>223747</td><td>100 mL</td></tr><tr><td>Serum Bottles</td><td>Wheaton</td><td>223762</td><td>20 mL</td></tr><tr><td>Sealing Cap 20 mm Aluminum</td><td>Wheaton</td><td>224183-01</td><td>Case of 1,000</td></tr><tr><td>20 mm gray butyl stopper (2-prong)</td><td>Wheaton</td><td>224100-192</td><td>Septum; Case of 1,000</td></tr><tr><td>Hand crimper</td><td>Wheaton</td><td>W225303</td><td>20 mm size</td></tr><tr><td>Hand Decapper</td><td>Wheaton</td><td>W225353</td><td>20 mm size</td></tr><tr><td>Acid vials</td><td>Wheaton</td><td>224881</td><td>0.50 dram size (2-mL)</td></tr><tr><td>Power supply</td><td>SR Components</td><td>DDU240060</td><td>Class 2 Transformer AC adaptor; Input 120VAC , Output 24VDC</td></tr><tr><td>Calcium carbonate</td><td>Fisher</td><td>471-34-1</td><td>500 g of 100% w/w CaCO&lt;sub&gt;3&lt;/sub&gt;</td></tr></tbody></table>

measuring and mapping patterns of soil erosion and deposition related to soil carbonate concentrations under agricultural management

Spatial patterns of soil erosion and deposition can be inferred from differences in ground elevation mapped at appropriate time increments. Such changes in elevation are related to changes in near-surface soil carbonate (CaCO3) profiles. The objective is to describe a simple conceptual model and detailed protocol for repeatable field and laboratory measurements of these quantities. Here, accurate elevation is measured using a ground-based differential global positioning system (GPS); other data acquisition methods could be applied to the same basic method. Soil samples are collected from prescribed depth intervals and analyzed in the lab using an efficient and precise modified pressure-calcimeter method for quantitative analysis of inorganic carbon concentration. Standard statistical methods are applied to point data, and representative results show significant correlations between changes in soil surface layer CaCO3 and changes in elevation consistent with the conceptual model; CaCO3 generally decreased in depositional areas and increased in erosional areas. Maps are derived from point measurements of elevation and soil CaCO3 to aid analyses. A map of erosional and depositional patterns at the study site, a rain-fed winter wheat field cropped in alternating wheat-fallow strips, shows the interacting effects of water and wind erosion affected by management and topography. Alternative sampling methods and depth intervals are discussed and recommended for future work relating soil erosion and deposition to soil CaCO3.

Mapping DEM differences from 2001 and 2009 reveals erosion (red) and deposition (green) over that 8-year period, with decimeter-level changes in elevation over most areas (Figure 5). At the field-scale, erosion is dominant in the west and southwest, while deposition is seen along a northwest to southeast diagonal band on the eastern side of the field. Alternating bands of erosion and deposition are seen at the management-scale, often with abrupt changes at th...

Watch this Scientific Journal Video about Measuring and Mapping Patterns of Soil Erosion and Deposition Related to Soil Carbonate Concentrations Under Agricultural Management at JoVE.com

Measuring and Mapping Patterns of Soil Erosion and Deposition Related to Soil Carbonate Concentrations Under Agricultural Management

Spatial patterns of soil erosion and deposition can be inferred from differences in ground elevation mapped at appropriate time increments. Such changes in elevation are related to changes in near-surface soil carbonates. Repeatable methods for field and laboratory measurements of these quantities and data analysis methods are described here.

Spatial patterns of soil erosion and deposition can be inferred from differences in ground elevation mapped at appropriate time increments. Such changes ...

measuring-mapping-patterns-soil-erosion-deposition-related-to-soil

Research

Education

JoVE Core

JoVE Journal

Civil Engineering

Environment

Unit 1

Geographic Information Systems (GIS)

SCIENTISTS IN ACTION

11.7K Views.  USDA Agricultural Research Service (ARS). Ruimtelijke patronen van bodemerosie en depositie kunnen van verschillen in hoogte van de grond toegewezen op tijdig stappen worden afgeleid. Dergelijke wijzigingen in de hoogte hebben betrekking op wijzigingen in in de buurt van-oppervlakte bodem carbonaten. Herhaalbare methoden voor veld en laboratorium metingen van deze hoeveelheden en data analysemethoden worden hier beschreven.

Video: Measuring and Mapping Patterns of Soil Erosion and Deposition Related to Soil Carbonate Concentrations Under Agricultural Management

Design and Operation of a Continuous 13C and 15N Labeling Chamber for Uniform or Differential, Metabolic and Structural, Plant Isotope Labeling

Tracing rare stable isotopes from plant material through the ecosystem provides the most sensitive information about ecosystem processes; from CO2 fluxes and soil organic matter formation to small-scale stable-isotope biomarker probing. Coupling multiple stable isotopes such as 13C with 15N, 18O or&#160;2H has the potential to reveal even more information about complex stoichiometric relationships during biogeochemical transformations. Isotope labeled plant material has been used in various studies of litter decomposition and soil organic matter formation1-4. From these and other studies, however, it has become apparent that structural components of plant material behave differently than metabolic components (i.e. leachable low molecular weight compounds) in terms of microbial utilization and long-term carbon storage5-7. The ability to study structural and metabolic components separately provides a powerful new tool for advancing the forefront of ecosystem biogeochemical studies. Here we describe a method for producing 13C and 15N labeled plant material that is either uniformly labeled throughout the plant or differentially labeled in structural and metabolic plant components.
Here, we present the construction and operation of a continuous 13C and 15N labeling chamber that can be modified to meet various research needs. Uniformly labeled plant material is produced by continuous labeling from seedling to harvest, while differential labeling is achieved by removing the growing plants from the chamber weeks prior to harvest. Representative results from growing Andropogon gerardii Kaw demonstrate the system&#39;s ability to efficiently label plant material at the targeted levels. Through this method we have produced plant material with a 4.4 atom%13C and 6.7 atom%15N uniform plant label, or material that is differentially labeled by up to 1.29 atom%13C and 0.56 atom%15N in its metabolic and structural components (hot water extractable and hot water residual components, respectively). Challenges lie in maintaining proper temperature, humidity, CO2 concentration,&#160;and light levels in an airtight 13C-CO2 atmosphere for successful plant production. This chamber description represents a useful research tool to effectively produce uniformly or differentially multi-isotope labeled plant material for use in experiments on ecosystem biogeochemical cycling.

Design and Operation of a Continuous 13C and 15N Labeling Chamber for Uniform or Differential, Metabolic and Structural, Plant Isotope Labeling

This method explains how to build and operate a continuous 13C and 15N isotope labeling chamber for uniform or differential plant tissue labeling. Representative results from metabolic and structural labeling of Andropogon gerardii are discussed.

Tracing rare stable isotopes from plant material through the ecosystem provides the most sensitive information about ecosystem processes; from CO2 fluxes ...

A Protocol for Conducting Rainfall Simulation to Study Soil Runoff

Rainfall is a driving force for the transport of environmental contaminants from agricultural soils to surficial water bodies via surface runoff. The objective of this study was to characterize the effects of antecedent soil moisture content on the fate and transport of surface applied commercial urea, a common form of nitrogen (N) fertilizer, following a rainfall event that occurs within 24&#160;hr&#160;after fertilizer application. Although urea is assumed to be readily hydrolyzed to ammonium and therefore not often available for transport, recent studies suggest that urea can be transported from agricultural soils to coastal waters where it is implicated in harmful algal blooms. A rainfall simulator was used to apply a consistent rate of uniform rainfall across packed soil boxes that had been prewetted to different soil moisture contents. By controlling rainfall and soil physical characteristics, the effects of antecedent soil moisture on urea loss were isolated. Wetter soils exhibited shorter time from rainfall initiation to runoff initiation, greater total volume of runoff, higher urea concentrations in runoff, and greater mass loadings of urea in runoff. These results also demonstrate the importance of controlling for antecedent soil moisture content in studies designed to isolate other variables, such as soil physical or chemical characteristics, slope, soil cover, management, or rainfall characteristics. Because rainfall simulators are designed to deliver raindrops of similar size and velocity as natural rainfall, studies conducted under a standardized protocol can yield valuable data that, in turn, can be used to develop models for predicting the fate and transport of pollutants in runoff.

A rainfall simulator was used to apply a consistent rate of uniform rainfall to packed soil boxes in a study of the fate and transport of urea, a nonpoint source environmental contaminant. Under uniform soil and rainfall conditions, antecedent soil moisture content exerted strong control over urea loss in surface runoff.

Rainfall is a driving force for the transport of environmental contaminants from agricultural soils to surficial water bodies via surface runoff. The ...

Measurement of Greenhouse Gas Flux from Agricultural Soils Using Static Chambers

Measurement of greenhouse gas (GHG) fluxes between the soil and the atmosphere, in both managed and unmanaged ecosystems, is critical to understanding the biogeochemical drivers of climate change and to the development and evaluation of GHG mitigation strategies based on modulation of landscape management practices. The static chamber-based method described here is based on trapping gases emitted from the soil surface within a chamber and collecting samples from the chamber headspace at regular intervals for analysis by gas chromatography. Change in gas concentration over time is used to calculate flux. This method can be utilized to measure landscape-based flux of carbon dioxide, nitrous oxide, and methane, and to estimate differences between treatments or explore system dynamics over seasons or years. Infrastructure requirements are modest, but a comprehensive experimental design is essential. This method is easily deployed in the field, conforms to established guidelines, and produces data suitable to large-scale GHG emissions studies.

This article showcases the static chamber-based method for measurement of greenhouse gas flux from soil systems. With relatively modest infrastructure investments, measurements may be obtained from multiple treatments/locations and over timeframes ranging from hours to years.

Measurement of greenhouse gas (GHG) fluxes between the soil and the atmosphere, in both managed and unmanaged ecosystems, is critical to understanding the ...

A CO2 Concentration Gradient Facility for Testing CO2 Enrichment and Soil Effects on Grassland Ecosystem Function

Continuing increases in atmospheric carbon dioxide concentrations (CA) mandate techniques for examining impacts on terrestrial ecosystems. Most experiments examine only two or a few levels of CA concentration and a single soil type, but if CA can be varied as a gradient from subambient to superambient concentrations on multiple soils, we can discern whether past ecosystem responses may continue linearly in the future and whether responses may vary across the landscape. The Lysimeter Carbon Dioxide Gradient Facility applies a 250 to 500 &#181;l L-1 CA gradient to Blackland prairie plant communities established on lysimeters containing clay, silty clay, and sandy soils. The gradient is created as photosynthesis by vegetation enclosed in in temperature-controlled chambers progressively depletes carbon dioxide from air flowing directionally through the chambers. Maintaining proper air flow rate, adequate photosynthetic capacity, and temperature control are critical to overcome the main limitations of the system, which are declining photosynthetic rates and increased water stress during summer. The facility is an economical alternative to other techniques of CA enrichment, successfully discerns the shape of ecosystem responses to subambient to superambient CA enrichment, and can be adapted to test for interactions of carbon dioxide with other greenhouse gases such as methane or ozone.

A CO2 Concentration Gradient Facility for Testing CO2 Enrichment and Soil Effects on Grassland Ecosystem Function

The Lysimeter Carbon Dioxide Gradient Facility creates a 250 to 500 &#181;l L-1 linear carbon dioxide gradient in temperature-controlled chambers housing grassland plant communities on clay, silty clay, and sandy soil monoliths. The facility is used to determine how past and future carbon dioxide levels affect grassland carbon cycling.

Continuing increases in atmospheric carbon dioxide concentrations (CA) mandate techniques for examining impacts on terrestrial ecosystems. Most ...

Measuring Carbon-based Contaminant Mineralization Using Combined CO2 Flux and Radiocarbon Analyses

A method is described which uses the absence of radiocarbon in industrial chemicals and fuels made from petroleum feedstocks which frequently contaminate the environment. This radiocarbon signal &#8212; or rather the absence of signal &#8212; is evenly distributed throughout a contaminant source pool (unlike an added tracer) and is not impacted by biological, chemical or physical processes (e.g., the 14C radioactive decay rate is immutable). If the fossil-derived contaminant is fully degraded to CO2, a harmless end-product, that CO2 will contain no radiocarbon. CO2 derived from natural organic matter (NOM) degradation will reflect the NOM radiocarbon content (usually &#60;30,000 years old). Given a known radiocarbon content for NOM (a site background), a two end-member mixing model can be used to determine the CO2 derived from a fossil source in a given soil gas or groundwater sample. Coupling the percent CO2 derived from the contaminant with the CO2 respiration rate provides an estimate for the total amount of contaminant degraded per unit time. Finally, determining a zone of influence (ZOI) representing the volume from which site CO2 is collected allows determining the contaminant degradation per unit time and volume. Along with estimates for total contaminant mass, this can ultimately be used to calculate time-to-remediate or otherwise used by site managers for decision-making.

Measuring Carbon-based Contaminant Mineralization Using Combined CO2 Flux and Radiocarbon Analyses

A protocol is described wherein CO2 mineralized from organic contaminant (derived from petroleum feedstocks) biodegradation is trapped, quantified, and analyzed for 14C content. A model is developed to determine CO2 capture zone's spatial extent. Spatial and temporal measurements allow integrating contaminant mineralization rates for predicting remediation extent and time.

A method is described which uses the absence of radiocarbon in industrial chemicals and fuels made from petroleum feedstocks which frequently contaminate ...

Methods of Soil Resampling to Monitor Changes in the Chemical Concentrations of Forest Soils

Recent soils research has shown that important chemical soil characteristics can change in less than a decade, often the result of broad environmental changes. Repeated sampling to monitor these changes in forest soils is a relatively new practice that is not well documented in the literature and has only recently been broadly embraced by the scientific community. The objective of this protocol is therefore to synthesize the latest information on methods of soil resampling in a format that can be used to design and implement a soil monitoring program. Successful monitoring of forest soils requires that a study unit be defined within an area of forested land that can be characterized with replicate sampling locations. A resampling interval of 5 years is recommended, but if monitoring is done to evaluate a specific environmental driver, the rate of change expected in that driver should be taken into consideration. Here, we show that the sampling of the profile can be done by horizon where boundaries can be clearly identified and horizons are sufficiently thick to remove soil without contamination from horizons above or below. Otherwise, sampling can be done by depth interval. Archiving of sample for future reanalysis is a key step in avoiding analytical bias and providing the opportunity for additional analyses as new questions arise.

Repeated soil sampling has recently been shown to be an effective way to monitor forest soil change over years and decades. To support its use, a protocol is presented that synthesizes the latest information on soil resampling methods to aid in the design and implementation of successful soil monitoring programs.

Recent soils research has shown that important chemical soil characteristics can change in less than a decade, often the result of broad environmental ...

Laboratory and Field Protocol for Estimating Sheet Erosion Rates from Dendrogeomorphology

Sheet erosion is among the crucial drivers of soil degradation. Erosion is controlled by environmental factors and human activities, which often lead to severe environmental impacts. The understanding of sheet erosion is, consequently, a worldwide issue with implications for both environment and economies. However, the knowledge on how erosion evolves in space and time is still limited, as well as its effects on the environment. Below, we explain a new dendrogeomorphological protocol for deriving eroded soil thickness (Ex) by acquiring accurate microtopographic data using both terrestrial laser scanning (TLS) and microtopographic profile gauges. Additionally, standard dendrogeomorphic procedures, dependent on anatomical variations in root rings, are utilized to establish the timing of exposure. Both TLS and microtopographic profile gauges are used to obtain ground surface pro&#64257;les, from which Ex is estimated after the threshold distance (TD) is determined, i.e., the distance between the root and the sediment knickpoint, which allows de&#64257;ning the lowering of the ground surface caused by sheet erosion. For each profile, we measured the height between the topside of the root and a virtual plane tangential to the ground surface. In this way, we intended to avoid small-scale impacts of soil deformation, which may be due to pressures exerted by the root system, or by the arrangement of exposed roots. This may provoke small amounts of soil sedimentation or erosion depending on how they physically affect the surface runoff. We demonstrate that an adequate microtopographic characterization of exposed roots and their associated ground surface is very valuable to obtain accurate erosion rates. This finding could be utilized to develop the best management practices designed to eventually halt or perhaps, at least, lessen soil erosion, so that more sustainable management policies can be put into practice.

Characterizing erosion from dendrogeomorphology has usually focused on accurately finding the starting time of root exposure, by examining macroscopic or cell level changes caused by exposure. Here, we offer a detailed description of different novel techniques to obtain more precise erosion rates from highly accurate microtopographic data.

Sheet erosion is among the crucial drivers of soil degradation. Erosion is controlled by environmental factors and human activities, which often lead to ...

Use of Principal Components for Scaling Up Topographic Models to Map Soil Redistribution and Soil Organic Carbon

Landscape topography is a critical factor affecting soil formation and plays an important role in determining soil properties on the earth surface, as it regulates the gravity-driven soil movement induced by runoff and tillage activities. The recent application of Light Detection and Ranging (LiDAR) data holds promise for generating high spatial resolution topographic metrics that can be used to investigate soil property variability. In this study, fifteen topographic metrics derived from LiDAR data were used to investigate topographic impacts on redistribution of soil and spatial distribution of soil organic carbon (SOC). Specifically, we explored the use of topographic principal components (TPCs) for characterizing topography metrics and stepwise principal component regression (SPCR) to develop topography-based soil erosion and SOC models at site and watershed scales. Performance of SPCR models was evaluated against stepwise ordinary least square regression (SOLSR) models. Results showed that SPCR models outperformed SOLSR models in predicting soil redistribution rates and SOC density at different spatial scales. Use of TPCs removes potential collinearity between individual input variables, and dimensionality reduction by principal component analysis (PCA) diminishes the risk of overfitting the prediction models. This study proposes a new approach for modeling soil redistribution across various spatial scales. For one application, access to private lands is often limited, and the need to extrapolate findings from representative study sites to larger settings that include private lands can be important.

Landscape processes are critical components of soil formation and play important roles in determining soil properties and spatial structure in landscapes. We propose a new approach using stepwise principal component regression to predict soil redistribution and soil organic carbon across various spatial scales.

Landscape topography is a critical factor affecting soil formation and plays an important role in determining soil properties on the earth surface, as it ...

Monitoring Pedogenic Inorganic Carbon Accumulation Due to Weathering of Amended Silicate Minerals in Agricultural Soils.

The present study aims to demonstrate a systematic procedure for monitoring inorganic carbon induced by enhanced weathering of comminuted rocks in agricultural soils. To this end, the core soil samples taken at different depth (including 0-15 cm, 15-30 cm, and 30-60 cm profiles) are collected from an agriculture field, the topsoil of which has already been enriched with an alkaline earth metal silicate containing mineral (such as wollastonite). After transporting to the laboratory, the soil samples are air-dried and sieved. Then, the inorganic carbon content of the samples is determined by a volumetric method called calcimetry. The representative results presented herein showed five folded increments of inorganic carbon content in the soils amended with the Ca-silicate compared to control soils. This compositional change was accompanied by more than 1 unit of pH increase in the amended soils, implying high dissolution of the silicate. Mineralogical and morphological analyses, as well as elemental composition, further corroborate the increase in the inorganic carbon content of silicate-amended soils. The sampling and analysis methods presented in this study can be adopted by researchers and professionals looking to trace pedogenic inorganic carbon changes in soils and subsoils, including those amended with other suitable silicate rocks such as basalt and olivine. These methods can also be exploited as tools for verifying soil inorganic carbon sequestration by private and governmental entities to certify and award carbon credits.

The verification method described here is adaptable for monitoring pedogenic inorganic carbon sequestration in various agricultural soils amended with alkaline earth metal silicate-containing rocks, such as wollastonite, basalt, and olivine. This type of validation is essential for carbon credit programs, which can benefit farmers that sequester carbon in their fields.

The present study aims to demonstrate a systematic procedure for monitoring inorganic carbon induced by enhanced weathering of comminuted rocks in ...

Protocol for Soil Hydraulic Property Measurement via Evaporation

Measurements of Soil Water Potential and Conductivity based on a Simple Evaporation Experiment using a Hydraulic Property Analyzer

The measurement of soil hydraulic properties is critical in understanding the physical components of soil health as well as integrated knowledge of soil systems under various management practices. Collecting reliable data is imperative for informing decisions that affect agriculture and the environment. The simple evaporation experiment described here uses instrumentation in a laboratory setting to analyze soil samples collected in the field. The soil water tension of the sample is measured by the instrument, and tension data is modeled by software to return soil hydraulic properties. This method can be utilized to measure soil water retention and hydraulic conductivity and give insight into differences in treatments or environmental dynamics over time. Initial establishment requires a user, but data acquisition is automated with the instrument. Soil hydraulic properties are not easily measured with traditional experiments, and this protocol offers a simple and optimal alternative. Interpretation of results and options for extending the data range are discussed.

Author Spotlight: Advancing Agricultural Land Ecosystem Research with a Hydraulic Property Analyzer to Assess Soil Health

This article features a simple evaporation experiment using a hydraulic property instrument for a soil sample. Through efficient means, measurements can be taken over a series of days to generate high-quality data.

The measurement of soil hydraulic properties is critical in understanding the physical components of soil health as well as integrated knowledge of soil ...

Measuring and Mapping Patterns of Soil Erosion and Deposition Related to Soil Carbonate Concentrations Under Agricultural Management

Samenvatting

Meer video's verkennen