Here we present a methodology which matches a soil sample size and a hydraulic conductivity measurement device to prevent the so-called wall flow along the inside of the soil container from being erroneously included in water flow measurements. Its use is demonstrated with samples collected from a wastewater irrigation site.
Since the early 1960s, an alternative wastewater discharge practice at The Pennsylvania State University has been researched and its impacts monitored. Rather than discharging treated wastewater to a stream, and thereby directly impacting the stream quality, the effluent is applied to forested and cropped land managed by the University. Concerns related to reductions in soil hydraulic conductivity occur when considering wastewater reuse. The methodology described in this manuscript, matching soil sample size with the size of the laboratory-based hydraulic conductivity measurement apparatus, provides the benefits of a relatively rapid collection of samples with the benefits of controlled laboratory boundary conditions. The results suggest that there may have been some impact of wastewater reuse on the soil's ability to transmit water at deeper depths in the depressional areas of the site. Most of the reductions in the soil hydraulic conductivity in the depressions appear to be related to the depth from which the sample was collected, and by inference, associated with the soil structural and textural differences.
Discharge of treated wastewater from municipalities into streams has been a standard practice for decades. Such wastewater is treated primarily for the purpose of reducing the potential for biological oxygen consumption by microorganisms in the receiving waters, as a result of the discharged wastewater effluent. Oxygen consumption by microorganisms degrades organic materials in the wastewater reducing the oxygen levels in the water body into which the effluent is discharged and thereby harm aquatic organisms, including fish.
In recent decades concerns have developed related to inorganic nutrients, some metals, and other chemicals within wastewater which create harm. Due to a study published by Kolpin et al.1, a greater focus on a range of chemicals not previously considered has evolved. This study, published by the United States Geological Society, raised awareness regarding the wide range of personal care products and other chemicals in rivers and streams across the US due to discharge from wastewater treatment facilities.
Since the early 1960s, researchers at Penn State University have investigated and developed an alternative wastewater discharge practice somewhat unique in a humid region. Rather than discharging treated wastewater to a stream, and thereby directly impacting the stream quality, the effluent is applied to the forested and the cropped land managed by the University. This application area, nicknamed "The Living Filter", presently accepts all wastewater effluent generated from the campus plus some from the municipality. This reduces the likelihood for excess nutrients to enter streams which deliver water to the Chesapeake Bay, protects the local cold-water fishery from discharges of warm wastewater which is harmful to the fish, and prevents the delivery of other chemicals contained in the wastewater from directly contacting aquatic ecosystems.
However, there are always consequences of behavior changes, and this alternative use facility is not immune to such. Questions have arisen regarding whether the application of the wastewater effluent has negatively impacted the soil's ability to allow water to infiltrate the soil surface2,3,4,5 and caused greater runoff, whether there is a possible contamination of the local wells with chemicals (nutrients, antibiotics or other pharmaceutical compounds, personal care products) contained in the wastewater effluent, and whether those chemicals are creating negative environmental impacts, such as through the uptake of chemicals into plants6 grown on the site, or the development of antibiotic resistance in soil organisms7 at the site.
As a result of some of these concerns, this study is conducted to determine the impacts of the irrigation of wastewater effluent on soil hydraulic conductivity at saturation. The approach used involves collecting soils from selected sites either within or outside the irrigated area and matching the soil sample container size with the laboratory setup. It is important for the soil sample container to fit into the laboratory apparatus and for the water that moves downward through the soil matrix in the sample to be separated from the water that moves downward between the soil and the soil sample container. The protocol describes how the laboratory apparatus is constructed to ensure this occurred.
Soil samples are collected using a hydraulic core sampler attached to a tractor. Soil cores are collected from selected areas in the undulating landscape and retained in a plastic sleeve fitted into the soil core sampler. These cores are collected from a Hagerstown silt loam, located either in a summit landscape position or in a depressional area. Six representative summits and six depressional sites are sampled from the irrigated area (a total of 12 irrigated area sampling sites). In addition, three summits and three depressional sites are sampled from an adjacent, non-irrigated area (a total of six non-irrigated sites). A maximum of six cores is collected at each site to a maximum depth of approximately 1,200 mm, with each core sample being approximately 150 mm long (100 mm of the sample being contained in the plastic sleeve and 50 mm being contained in the cutting head of the metal sampler). After removal from the metal sampler, the plastic sleeves containing the collected soil cores are fitted with end caps, transported upright to the laboratory, and stored upright until they are used to determine the saturated hydraulic conductivity. Concurrently, soil samples are collected at each depth for the determination of soil and soil solution concentrations of Calcium (Ca), Magnesium (Mg), and Sodium (Na) using a Mehlich 3 extraction for estimates of soil concentrations8 and deionized water extracts at a 1:2 ratio of soil mass:water mass. The chemical analyses of the water extracts were obtained from Inductively Coupled Plasma Atomic Emission Spectroscopy (ICP-AES) and were used to calculate the Sodium Adsorption Ratio (SAR).
The determination of the saturated hydraulic conductivity is carried out primarily using a constant head method9. A solution containing Ca and Na salts to mimic the effluent electrical conductivity (EC) and SAR of the effluent is created so the soil will be exposed to water quality variables similar to the wastewater applied in the field. In this case, the EC is 1.3 dS/m and the SAR is 3, reflecting the EC and SAR of the effluent in recent years prior to the sample period. [Technically, the units for SAR are (milliequivalents/liter)½ and are not usually identified in the literature.]
The modification to the constant head method of Klute and Dirksen9 is the development of a flow separator by Walker8 to prevent flow through the column which occurred outside the soil matrix from being included in the estimate of the soil hydraulic conductivity. The flow separator is built using polyvinyl chloride (PVC) tubing selected and machined to match the soil sample size. A screen supports the soil sample and allows the water that has moved through the soil matrix to flow out the bottom of the sample. A second outlet emits the water that has flowed down the inside of the plastic sleeve, thereby eliminating so-called "wall flow" from being incorrectly included in the estimate of the amount of water that moves through the soil matrix.
1. Selecting Soil Sampling Locations
2. Collecting Soil Samples
Supplemental Figure 1: Drill rig used for sampling.
3. Building a Constant Head, Multiple Column, Soil Hydraulic Conductivity Setup
NOTE: The hydraulic conductivity laboratory apparatus is based on work by Walker10. It involves the use of a permeameter which is constructed to separate the flow in between the outer edge of the sample and the cylinder containing the ring from the flow through the soil matrix. The ID of any PVC pipe referred to below is not a strict tolerance. Some may fit well, and others may require some work (light sanding).
Supplemental Figure 2: Front view of saturated hydraulic conductivity apparatus. Please click here to view a larger version of this figure.
Supplemental Figure 3: End view of water supply gutter. Please click here to view a larger version of this figure.
4. Obtaining Soil Hydraulic Conductivity Values
To investigate the question of whether the application of wastewater effluent at the Living Filter site has impacted the ability of the soil to transmit water, we conducted experiments to measure the saturated hydraulic conductivity of the soils. We compared hydraulic conductivity of soils from irrigated areas of the site with those in non-irrigated areas of the site. The impact of wastewater effluent on the soil hydraulic conductivity is a question of concern, as there have been some reports in less humid regions of reductions in the ability of the soil to transmit water resulting from (for example) a buildup of sodium in the soil or from the development of a surface biological crust. For each of the samples collected for hydraulic conductivity measurements, samples from adjacent locations within a meter of the hydraulic conductivity samples were collected for a measurement of the major cations of Ca, Mg, and Na in the soil solution. The values of Ca, Mg, and Na in the soil solution were estimated from 1:2 soil mass:water mass ratio deionized water extracts of the soils, with values back-calculated based on the soil water content at the time of sampling. These soil solution concentrations were used to calculate the sodium adsorption ratio of the soil solution, as defined by:
where the values for Na, Ca, and Mg are given in milliequivalents (meq)/liter. (For Na, the number of meq/l = mg/l divided by 23; for Ca, meq/l = mg/l divided by 20; for Mg, meq/l = mg/l divided by 12.)
It was observed that the hydraulic conductivity decreased with depth in the depressional areas (Figure 1), with the strongest relationship of Ksat by depth existing in the non-irrigated areas.
When considering the relationship of Ksat with SAR, for soil samples below 20 cm, there was a strong positive relationship between Ksat and the SAR of the soil solution for the non-irrigated depressions (Figure 2), but a strong negative relationship between Ksat and SAR for the irrigated depressions. It should be noted that the values of SAR in the non-irrigated area were significantly below those in the irrigated area and were in a range in which it is not expected that the hydraulic conductivity of the soil would be impacted. The SAR values in the irrigated area were much higher and although they were in a range not expected to create any reductions in the hydraulic conductivity, it is observed that there was a strong negative relationship between Ksat and SAR (Figure 2).
Looking at the distribution of SAR by depth for the non-irrigated depressions (Figure 3), there was little relationship between soil solution SAR and soil depth. Thus, it is likely that the relationship between the decreasing Ksat values by depth is dominated by changes in the soil characteristics (structure, texture), which tend to change with depth. In contrast, Figure 4 demonstrates that there is a strong relationship between soil depth and soil solution SAR for irrigated depressional areas. The strong relationship between Ksat and SAR for the irrigated area is likely due to both a strong correlation between soil depth and soil solution SAR and to some extent to a decrease in Ksat with an increase in depth. It does appear, however, that there may be some decrease in Ksat simply due to the increased SAR.
Figure 1: Relationship between saturated hydraulic conductivity and soil depth for both irrigated and non-irrigated depressions.
Figure 2: Relationship between saturated hydraulic conductivity and sodium adsorption ratio (SAR) of the soil solution for both irrigated and non-irrigated depressions.
Figure 3: Value of the average SAR in soil solution by sampling depth for sites located in non-irrigated depressions.
Figure 4: Value of the soil solution SAR by depth for sites located in irrigated depressions.
The ability to collect field-based, undisturbed soil samples and obtain their hydraulic conductivity values is important in obtaining data representative of a site. In order to best represent field conditions, it is important to use soil samples which remain in a physical state representative of their environment in the field. Soil samples collected from a field site which are then disturbed by subsampling or by handling induced compaction, for example, will experience structural changes which impact the saturated hydraulic conductivity.
It is also important to have a means of measuring the soil hydraulic conductivity in a controlled laboratory setting. However, using a laboratory method for the saturated hydraulic conductivity which fails to account for a flow along the inside of the sample container (so-called "wall flow") will result in a poor reproduction of results, and high variability due to methodology rather than due to natural variability.
Since the development of the Disc Permeameter in 1982 by Perroux11, based on work by Clothier and White12, much effort has been expended to more satisfactorily obtain field-based soil hydraulic conductivity measurements13. Having field-based measurements is highly desirable, as "disturbed" soil samples (for example, dried and ground samples repacked into a column) do not reflect the natural soil conditions.
However, there are also drawbacks to field-based methods to obtain soil hydraulic conductivities. One assumption made with field-based methods is that the soil underlying the equipment used to measure the rate of the water intake into the soil is uniform13. However, most soils are not uniform but are comprised of layers of soil materials which differ in their conductivities.
Another drawback is the fact that an extensive sampling campaign can require from 1–4 h (or more) measurement time per sample, in addition to the site preparation time. Earlier work at this site14 required several weeks to complete using the method of Ankeny et al13. The consequence is that the collection of a large number of samples will require such a substantial period of time that the field conditions will change (e.g., plant growth, water content, etc.), and the sampling may also interfere with field operations (e.g., in this case, wastewater irrigation applications and harvest). Differences in environmental conditions (e.g., rainfall) can result in changes in soil chemical properties. In the case of this experiment, soil Ca, Mg, and Na concentrations changed because of the leaching of Ca, Mg, and Na through the soil due to rainfall infiltration and downward water movement.
Because of the physical labor involved in and the length of time required to prepare the site when it is covered with vegetation14, and the length of time required to collect field-based saturated hydraulic conductivity values, the ability to obtain representative values over a range of soil depths and field locations can require up to a half day per sample per depth. Necessary plant production field operations, including irrigation, may constrain the length of time for sample collection.
In addition, even if the time required to conduct field measurements of saturated hydraulic conductivity at many sites is available, the length of time to obtain samples across a field site and at several depths will necessarily result in many samples being collected under different environmental conditions, as those change from day to day (or more frequently).
Samples collected from the field with a hydraulic soil sampler can be collected in a much shorter time, thereby reducing the changes which may take place at a field site over time. However, the laboratory procedures to obtain the soil hydraulic conductivity from such samples have the distinct disadvantage of being subject to so-called "wall flow"10. Wall flow is the flow of water along the interior of the sample container when the sample is placed on a constant head device which is typically used to obtain soil hydraulic conductivity estimates. Such flow, if it is included in the measure of the rate of water movement through the soil, results in an erroneously high estimate of hydraulic conductivity. This manuscript describes the use of a laboratory set up to eliminate wall flow from estimates of soil hydraulic conductivity, and of a sample collection method which matches the soil sample size to the size of the laboratory equipment.
A critical step is to collect soil samples which have not been compacted. Although the soil moisture status will impact the resistance to the insertion of the sampler and thus the compatibility of the soil, the recommendation is that the length of the sample collected should be somewhat shorter than the sample liner which is inserted into the metal sample tube.
Transporting samples from the field to the lab should be done in a manner which minimizes disturbance to them. Keeping them upright and ensuring they are not tightly packed against each other will help reduce handling disturbances.
The most critical step in the protocol is to build the laboratory apparatus to match the sample size collected from the field so that wall flow is not included in the water collected from the soil matrix10. Although the description of the laboratory apparatus presented herein is for a particular size of sample container, other sized containers may be used if the sample holders in the laboratory apparatus are similarly matched in size.
Once a prototype is assembled, test samples intentionally created to result in wall flow should be utilized to ascertain that the construction of the apparatus does truly separate wall flow from soil matrix flow. Another important observation is whether the final design permits the establishment of a constant head of water on top of the soil sample without overtopping the plastic soil container. The top of the soil container must be above the water level in the supply gutter. This is critical. If the water overtops the plastic soil container, then the dimensions were likely not correctly measured. This can be overcome by fitting a rubber ring to the top of the plastic soil container, being careful not to disturb the soil sample.
The sample collection time required will be dependent upon both the criteria set for the experiment's precision as well as the rate of water movement through the soil. For example, samples may need to be collected for 10–20 min every hour for a period of 12 h until a relatively constant amount of water moves through the soil core and into the sampling container each time a sample is taken. In other cases, the samples may only need to be collected for 8–10 min for a period of 3 or 4 h before a constant amount of water moves through the soil for a given sample period. Constant volumes of water over the same period of sample time would indicate a "steady state" condition had been reached.
Soil cores were prepared for a saturated hydraulic conductivity analysis by placing cheesecloth on the bottom of each core and then placing the core in a tub of water to slowly saturate the samples from the bottom up, for a period of at least 24 h.
After presaturation, the cores were removed from the tub and the bottom of each of the cores was set on a flow separator designed to separate any flow along the sides of the plastic sleeve containing the soil, from flow through the soil itself. Six soil cores at a time were placed on this device which included a bench with a gutter system10 containing a weir for the delivery of water to the top of the soil samples at a constant head through siphon tubes. Water was pumped to the gutter system from a reservoir using a submersible pump.
The flow separator is essentially a 100 mm long and 100 mm diameter PVC tube which serves as a holder in which the plastic sleeve holding the soil sample sits. A second PVC tube (approximately 75 mm in diameter and 75 mm long) is sharpened so that the soil sample contacts the sharp edge of this PVC tube, and fits inside the outside PVC tube, with the plastic sleeve holding the soil sample fitting outside the smaller PVC tube. A screen in the smaller PVC tube supports the soil sample and allows the water that has moved through the soil to flow out the bottom of the sample. A second outlet emits the water that has flowed down the inside of the plastic sleeve, thereby eliminating so-called "bypass flow" from being incorrectly included in the estimate of the amount of water that moved through the soil matrix.
One important limitation to the technique occurs with soils, such as those with a high clay content, that have a low saturated hydraulic conductivity. Soils with very low saturated hydraulic conductivity typically must have their conductivity determined with a "falling head" approach7 instead of the constant head approach used here. The apparatus described here would need to be significantly modified to allow for a falling head approach to be utilized.
The design has been found to provide more consistent results of saturated soil hydraulic conductivity10 than traditional laboratory methods9. Use of the design should help reduce the frequency of recording erroneously high estimates of saturated soil hydraulic conductivity due to erroneously including wall flow in estimates of the volume of water flowing through the soil over a period of time.
The authors have nothing to disclose.
The authors would like to thank the Pennsylvania State University Office of Physical Plant for providing partial funding to support this project. Partial funding was also provided by the USDA-Regional Research Project W-3170. We would like to express our gratitude to Ephraim Govere for his assistance with the analytical work. Our deepest gratitude is to Charles Walker, whose engineering design and construction skills made it possible for us to conduct this work.
Sampling equipment: | |||
Soil Sampler Drill Rig | Giddings Machine Co. Inc | #25-TS / Model HDGSRTS | * NOTE: This model is comparable to the model we utilized but which is no longer produced |
Kelly Bar | Giddings Machine Co. Inc | #KB-208 8 Ft. Kelly Bar | |
Soil Sample Collection Tube | Giddings Machine Co. Inc | #ZC-180 4-3/4” X 7-1/4” | |
Soil Collection Tube Bit | Giddings Machine Co. Inc | #ZC-190 4-3/4” Standard Relief | |
Plastic Liner for Soil Sample | Giddings Machine Co. Inc | #ZC-208 3-5/8” x 6” | Enough for the number of samples being collected |
Black end caps a for bottom of sample liners | Giddings Machine Co. Inc | To retain samples in liners | |
Red end caps a for top of sample liners | Giddings Machine Co. Inc | To retain samples in liners | |
Cooler Chest | Store & maintain samples upright in sample liners during transport from field to lab | ||
Protective gear: | |||
Hardhats, googles, and gloves | other items as needed for personal protection | ||
Saw | |||
Drill and bits | |||
PVC Cement | |||
6 to 8 – 19 mm x 184 mm x 2438 mm boards | |||
2 – barbed fittings; 13 mm HB x MGHT | to connect plastic tubing to supply gutter and to drainage gutter | ||
6 – barbed fitting | to connect plastic tubing to outer PVC cylinder to allow for water drainage | ||
3000 mm long – 19 mm OD / 13 mm ID plastic tubing | |||
6 – 85 mm diameter circular mesh pieces | Can be cut from (e.g.) a 600 mm long, 6 mm x 18 gauge wire mesh (e.g. galvanized steel gutter guard) | ||
Schedule 40 PVC pipe – 96 mm ID / 114 mm OD | |||
Schedule 40 PVC pipe – 73 mm ID / 89 mm OD | |||
Schedule 40 PVC pipe – 63 mm ID / 73 mm OD, OR 6 – 73 mm plastic shower drains | |||
Schedule 40 PVC pipe – 25 mm ID | |||
6 – 6 mm thick x 155 mm square sheets of PVC | Can purchase 2 – 6 mm x 300 mm (appx) sheets for about $20 each from: https://www.interstateplastics.com/Pvc-Gray-Sheet-PVCGE~~SH.php?vid=20180212222911-7p | ||
6 – 140 mm by 19 mm plastic funnels | To direct water flowing from soil sample into collection beaker | ||
Adhesive caulk | |||
1 – length of 150 mm x 1200 mm wire mesh cloth | 4 Mesh works well | ||
2 – 120 mm x 1219 mm plastic gutter with end caps | |||
4 – gutter hangers | |||
1 – additional gutter end cap | To be cut as described in procedures to create a constant head in the supply gutter | ||
1 – large plastic tub | Appx 65 L in volume, for example, to serve as water source for the hydraulic conductivity procedure | ||
1 – large plastic tub | To serve for wetting up soil samples | ||
1 – Submersible pump | e.g. Beckett M400 AUL or M400 AS | ||
Plastic tubing | Various sized drainage tubes, water supply tube, and drain from drainage gutter | ||
Container of Cheese Cloth | To place at bottom of soil sample help retain soil in plastic sample container during hydraulic conductivity and wetting up | ||
Rubber bands | Large enough to fit around plastic sample liners tightly | ||
Scale which measures to at least 0.1 gram | |||
Beaker or other container to collect water from each sample | |||
Sodium Chloride | For creating a water quality similar to that which is typically applied to the soil | ||
Calcium Chloride | For creating a water quality similar to that which is typically applied to the soil |