PERCOLATE WATER AND BROMIDE MOVEMENT IN THE ROOT ZONE OF EFFLUENT IRRIGATION SITES1

ABSTRACT: A bromide tracer was used to evaluate percolate water and ion movement in the upper 1.2 m of soil at a proposed sewage effluent irrigation site located in the Missouri Ozarks. Two plots representing Doniphan silt loam and Crider silt loam soils were sprinkler irrigated with local ground water at a rate of 7.62 cm/week from June through August 1976. Soil water potential, percent soil moisture by volume, and background levels of bromide in soil water, ground water, and precipitation were measured at the study plots. Bromide exchange properties and saturated hydraulic conductivity of the soils were determined in the laboratory. During two selected time periods, irrigation water, was spiked with NaBr (5.0 mg/l Br). Bromide movement through the upper profile was quantified by soil water samples and post-sampling neutron activation analysis. Soil moisture was near saturatin in both soils when the Br tracer was applied. Bromide concentrations above background levels (0.023 mg/l Br, Doniphan silt loam and 0.016 mg/l Br, Crider silt loam) were detected within 2.60 hours at 0.9 m in the Doniphan soil and within 3.75 hours at that depth in the Crider soil. The rate of Br movement in the profile was greater in both soils than the measured saturated hydraulic conductivity, Bromide concentrations above background levels were present in soil water from the study plots for a minimum of 21 days after irrigation with the Br tracer.     

MODELING PERCOLATION LOSSES FROM A PONDED FIELD UNDER VARIABLE WATER-TABLE CONDITIONS1

ABSTRACT: A numerical simulation model was developed to predict the vertical and lateral percolation losses from a ponded agricultural field. The two-dimensional steady-state unsaturated/ saturated flow equation was solved using the finite-difference technique. A constant ponding depth was maintained at the soil surface with different water table conditions in an application of the model for rice fields bordered by bunds. Field experiments were conducted for two different water table depths to collect data on the spatial distribution of volumetric soil-moisture content for model verification. The measured soil-moisture content values were found to be in close agreement with those predicted by the model. The sensitivity analysis of the model with selected hydrologic conditions shows that the model is most sensitive to the values of saturated hydraulic conductivity, but relatively less sensitive to water table depth, ponding depth, and evaporation rate from the soil surface. It implies that, in a ponded rice field condition, the lateral and vertical percolation losses are mostly governed by the hydraulic conductivity of the soil. The vertical percolation losses were almost equal to the saturated hydraulic conductivity values and, in most cases, these losses increased with deeper water table depths. The lateral percolation losses also increased with deeper water table depths; however, these losses were relatively small in comparison to the vertical percolation losses. The vertical and lateral percolation losses increased with the increase in ponding depths. The lateral percolation losses through the bund decreased when the evaporation losses increased from the soil surface. The results of this study indicate that the percolation losses from a ponded field may be predicted accurately for a wide range of soil and hydrological conditions when the values of hydraulic conductivity, evaporation rate, depth of ponding, and water table depth are accurately known.     

Productivity equation for reclaiming surface mines

This article addresses the development of an agricultural productivity equation for predicting new soil (neo-sol) plant growth potential in Clay County, Minnesota, USA. Soil factors examined in the study include percent organic matter, percent slope, percent rock fragments, hydraulic conductivity, electrical conductivity, pH, topographic position, available water-holding capacity, bulk density, and percent clay. Squared terms and two-factor interaction terms were also examined as possible regressors. A best equation was selected that had a multiple coefficient of determination of 0.7399 and has five significant regressors and intercept withP.0001. The regressors are hydraulic conductivity, percent slope squared, bulk density times percent rock fragments, electrical conductivity times percent rock fragments, and electrical conductivity times percent organic matter. The regressors predict soil suitability for a general crop model. The crops included in the model are wheat, oats, barley, soybeans, sugar beets, sunflowers, and grasses/legumes.     

VARIABILITY OF SOIL WATER PROPERTIES AND CROP YIELD IN A SLOPED WATERSHED1

ABSTRACT: Spatial distribution of soil and water properties and the correlations between them and crop yield were determined for a natural rainfall environment. Hydraulic conductivity, soil texture, water retention, and soil-water flux were variables used to investigate their relationship to crop yield using multiple regression techniques. Variations in crop yields on a watershed with a 3 to 4 percent slope and moderately erosive soils were related to soil-water characteristics and soil properties along slope and with depth. Climatic conditions to sustain crop growth and yield ranged from inadequate soil water in 1983 to adequate soil water in 1984. Crop yield was predicted with models using both available and measured soil-water content. Available water content provided a better model for the prediction of water yield and does not require field measurements of actual soil-water content. Soil water holding capacity was more significant for predicting crop yield in soils with moderate to high silt content than infiltrability of water into the soil.     

Modeling the Spatially Varying Water Balance Processes in a Semiarid Mountainous Watershed of Idaho1

Stratton, Benjamin T., Venakataramana Sridhar, Molly M. Gribb, James P. McNamara, and Balaji Narasimhan, 2009. Modeling the Spatially Varying Water Balance Processes in a Semiarid Mountainous Watershed of Idaho. Journal of the American Water Resources Association (JAWRA) 45(6):1390‐1408. Abstract: The distributed Soil Water Assessment Tool (SWAT) hydrologic model was applied to a research watershed, the Dry Creek Experimental Watershed, near Boise Idaho to investigate its water balance components both temporally and spatially. Calibrating and validating SWAT is necessary to enable our understanding of the water balance components in this semiarid watershed. Daily streamflow data from four streamflow gages were used for calibration and validation of the model. Monthly estimates of streamflow during the calibration phase by SWAT produced satisfactory results with a Nash Sutcliffe coefficient of model efficiency 0.79. Since it is a continuous simulation model, as opposed to an event‐based model, it demonstrated the limited ability in capturing both streamflow and soil moisture for selected rain‐on‐snow (ROS) events during the validation period between 2005 and 2007. Especially, soil moisture was generally underestimated compared with observations from two monitoring pits. However, our implementation of SWAT showed that seasonal and annual water balance partitioning of precipitation into evapotranspiration, streamflow, soil moisture, and drainage was not only possible but closely followed the trends of a typical semiarid watershed in the intermountain west. This study highlights the necessity for better techniques to precisely identify and drive the model with commonly observed climatic inversion‐related snowmelt or ROS weather events. Estimation of key parameters pertaining to soil (e.g., available water content and saturated hydraulic conductivity), snow (e.g., lapse rates, melting), and vegetation (e.g., leaf area index and maximum canopy index) using additional field observations in the watershed is critical for better prediction.     

流域非点源分布式模型AnnAGNPS参数的不确定性研究

流域非点源分布式模型的不确定性研究成为目前的热点问题之一。本研究以宁波市章溪河流域为研究区,通过GIS软件AreView和AnnAGNPS模型,集成利用莫里斯分类筛选法,选取临界源面积、SCS曲线系数、水土保持因子、饱和导水率、氮磷吸收率、施肥量等6个因子进行敏感性分析。模拟结果表明,土壤饱和导水率、水土保持因子和径流曲线CN值对模型输出结果影响较大;水土保持因子对泥沙、总磷和总有机碳负荷的模拟结果影响最大,均呈显著负效应;径流曲线CN值对总氮负荷计算结果影响最大,呈显著正效应。     

Characterizing scale- and location-dependent correlation of water retention parameters with soil physical properties using wavelet techniques

Understanding the correlation between soil hydraulic parameters and soil physical properties is a prerequisite for the prediction of soil hydraulic properties from soil physical properties. The objective of this study was to examine the scale- and location-dependent correlation between two water retention parameters (alpha and n) in the van Genuchten (1980) function and soil physical properties (sand content, bulk density [Bd], and organic carbon content) using wavelet techniques. Soil samples were collected from a transect from Fuxin, China. Soil water retention curves were measured, and the van Genuchten parameters were obtained through curve fitting. Wavelet coherency analysis was used to elucidate the location- and scale-dependent relationships between these parameters and soil physical properties. Results showed that the wavelet coherence between alpha and sand content was significantly different from red noise at small scales (8-20 m) and from a distance of 30 to 470 m. Their wavelet phase spectrum was predominantly out of phase, indicating negative correlation between these two variables. The strong negative correlation between alpha and Bd existed mainly at medium scales (30-80 m). However, parameter n had a strong positive correlation only with Bd at scales between 20 and 80 m. Neither of the two retention parameters had significant wavelet coherency with organic carbon content. These results suggested that location-dependent scale analyses are necessary to improve the performance for soil water retention characteristic predictions.     

SURFACE WATER QUALITY IMPACTS OF TILLAGE PRACTICES UNDER LIQUID SWINE MANIJRE APPLICATION1

ABSTRACT: Field studies were conducted to investigate the effects of tillage practices on the saturated hydraulic conductivity, and quantity and quality of surface runoff water resulting from the application of the liquid swine manure as a fertilizer. As part of the study, infiltration experiments were conducted on silt-loam soil with no-tillage (NT) and disk tillage (DT) practices. Liquid swine manure was applied on test plots, and the rainfall was applied by the portable rainfall simulator. The infiltration data was analyzed for the saturated hydraulic conductivity (K8) and runoff volume determinations. The surface runoff water was analyzed for total N, total P, ammonia, and nitrate concentration determinations. The study indicated that the tillage had significant effects on K_s, and quantity and quality of runoff water. The K_s values of the NT plots were found to vary from 0.693 to 1.734 mm/min, with a mean of 1.494 mm/min, while they varied from 1.056 to 2.543 mm/min, with a mean of 2231 mm/mm in the DT plots. The total N, total P, ammonium nitrogen and nitrate nitrogen concentrations were lower in runoff generated from DT plots, compared to that from the NT plots. The chemical concentration levels were significantly different in runoff waters collected one-day after manure application than in those collected 40-days after the manure application. Study suggested that the DT practice must be preferred over the NT practice if liquid swine manure is used as the fertilizer. The study is further continued to assess the long-term impacts of swine manure application and tillage on the quantity and quality of surface runoff water.     

Long-term effect of a single application of organic refuse on carbon sequestration and soil physical properties

Restoration of degraded lands could be a way to reverse soil degradation and desertification in semiarid areas and mitigate greenhouse gases (GHG). Our objective was to evaluate the long-term effects of a single addition of organic refuse on soil physical properties and measure its carbon sequestration potential. In 1988, a set of five plots (87 m(2) each) was established in an open desert-like scrubland (2-4% cover) in Murcia, Spain, to which urban solid refuse (USR) was added in a single treatment at different rates. Soil properties were monitored over a 5-yr period. Sixteen years after the addition, three of the plots were monitored again (P0: control, P1: 13 kg m(-2), P2: 26 kg m(-2) of USR added) to assess the lasting effect of the organic addition on the soil organic carbon (SOC) pools and on the physical characteristics of the soil. The SOC content was higher in P2 (16.4 g kg(-1)) and in P1 (11.8 g kg(-1)) than in P0 (7.9 g kg(-1)). Likewise, aerial biomass increased from 0.18 kg m(-2) in P0 up to 0.27 kg m(-2) in P1 and 0.46 kg m(-2) in P2. This represents a total C sequestration of 9.5 Mg ha(-1) in P2 and 3.4 Mg ha(-1) in P1, most of the sequestered C remaining in the recalcitrant soil pool. Additionally, higher saturated hydraulic conductivity, aggregate stability, and available water content values and lower bulk density values were measured in the restored plots. Clearly, a single addition of organic refuse to the degraded soils to increase the potential for C sequestration was effective.     

MEASUREMENT OF HYDRAULIC CONDUCTIVITY AND DIFFUSIVITY FOR PREDICTING THE PROCESS OF SOIL WATER INFILTRATION FROM A TRICKLE SOURCE1

ABSTRACT: Two soil water functions, hydraulic conductivity K(θ) and diffusivity D(θ), were estimated by two methods In one method D(θ) was estimated according to Bruce and Klute (1956), and K(θ) was calculated from D(θ) and the retention curve. In the second, K(θ) was obtained by field estimation, with D(θ) being calculated from K(θ) and the retention curve. The criterion of reliability for both methods was agreement between experimental and predicted distribution of soil water content. The prediction was made using the functions K(θ) and D(θ) as soil water parameters in both methods. Theoretical and experimental agreement was generally good. The first method, however, was found to be best for high soil water content and the second for low soil water content. In addition, the water content at the end of the monotonic increase of function D(θ) (estimated according to Bruce and Klute 1956) was found to be about the upper limit of field soil water content. It can be used as a boundary condition in the numerical solution of a cylindrical model of infiltration from a trickle source. It was concluded that the best agreement between theory and experiment can be found when the combined values of D(θ) and K(θ) from both methods of estimation are used.     

USE OF SOIL MOISTURE PROBES TO ESTIMATE GROUND WATER RECHARGE AT AN OIL SPILL SITE1

Soil moisture data collected using an automated data logging system were used to estimate ground water recharge at a crude oil spill research site near Bemidji, Minnesota. Three different soil moisture probes were tested in the laboratory as well as the field conditions of limited power supply and extreme weather typical of northern Minnesota: a self‐contained reflectometer probe, and two time domain reflectometry (TDR) probes, 30 and 50 cm long. Recharge was estimated using an unsaturated zone water balance method. Recharge estimates for 1999 using the laboratory calibrations were 13 to 30 percent greater than estimates based on the factory calibrations. Recharge indicated by the self‐contained probes was 170 percent to 210 percent greater than the estimates for the TDR probes regardless of calibration method. Results indicate that the anomalously large recharge estimates for the self‐contained probes are not the result of inaccurate measurements of volumetric moisture content, but result from the presence of crude oil, or borehole leakage. Of the probes tested, the 50 cm long TDR probe yielded recharge estimates that compared most favorably to estimates based on a method utilizing water table fluctuations. Recharge rates for this probe represented 24 to 27 percent of 1999 precipitation. Recharge based on the 30 cm long horizontal TDR probes was 29 to 37 percent of 1999 precipitation. By comparison, recharge based on the water table fluctuation method represented about 29 percent of precipitation.     

Tillage and field scale controls on greenhouse gas emissions

There is a lack of understanding of how associations among soil properties and management-induced changes control the variability of greenhouse gas (GHG) emissions from soil. We performed a laboratory investigation to quantify relationships between GHG emissions and soil indicators in an irrigated agricultural field under standard tillage (ST) and a field recently converted (2 yr) to no-tillage (NT). Soil cores (15-cm depth) were incubated at 25 degrees C at field moisture content and 75% water holding capacity. Principal component analysis (PCA) identified that most of the variation of the measured soil properties was related to differences in soil C and N and soil water conditions under ST, but soil texture and bulk density under NT. This trend became more apparent after irrigation. However, principal component regression (PCR) suggested that soil physical properties or total C and N were less important in controlling GHG emissions across tillage systems. The CO2 flux was more strongly determined by microbial biomass under ST and inorganic N content under NT than soil physical properties. Similarly, N2O and CH4 fluxes were predominantly controlled by NO3- content and labile C and N availability in both ST and NT soils at field moisture content, and NH4+ content after irrigation. Our study indicates that the field-scale variability of GHG emissions is controlled primarily by biochemical parameters rather than physical parameters. Differences in the availability and type of C and N sources for microbial activity as affected by tillage and irrigation develop different levels and combinations of field-scale controls on GHG emissions.     

STREAM DISCHARGE PREDICTION VIA A GRID BASED SOIL WATER ROUTING WITH PADDY FIELDS1

ABSTRACT: A grid based daily hydrologic model for a watershed with paddy fields was developed to predict the stream discharge. ASCII formatted elevation, soil, and land use data supported by the GRASS Geographic Information System are used to generate distributed results such as surface runoff and subsurface flow, soil water content, and evapotranspiration. The model uses a single flow path algorithm and simulates a water balance at each grid element. A linear reservoir assumption was used to predict subsurface runoff components. The model was applied to a 75.6 km² watershed located in the middle of South Korea, and observed stream flow hydrographs from 1995 and 1996 were compared to model predictions. The stream flow predictions of 1995 and 1996 generally agreed with the observed flow, resulting in a Nash‐Sutcliffe efficiency R² of 0.60 and 0.62, respectively. The hydraulic conductivity for percolating water through the saturated layer affected baseflow generation. The levee height of the paddy influenced the time and magnitude of the surface runoff, depending on irrigation management. The model will be used for making low flow management decisions by evaluating the role of each land use to stream flow, especially in case of paddy decrease by gradual urbanization of a watershed.     

Field Trial Assessment of Biological,Chemical, and Physical Responses of Soil to Tillage Intensity,Fertilization, and Grazing

Soil microbial populations can fluctuate in response to environmental changes and, therefore, are often used as biological indicators of soil quality. Soil chemical and physical parameters can also be used as indicators because they can vary in response to different management strategies. A long-term field trial was conducted to study the effects of different tillage systems (NT: no tillage, DH: disc harrow, and MP: moldboard plough), P fertilization (diammonium phosphate), and cattle grazing (in terms of crop residue consumption) in maize (Zea mays L.), sunflower (Heliantus annuus L.), and soybean (Glycine max L.) on soil biological, chemical, and physical parameters. The field trial was conducted for four crop years (2000/2001, 2001/2002, 2002/2003, and 2003/2004). Soil populations of Actinomycetes, Trichoderma spp., and Gliocladium spp. were 49% higher under conservation tillage systems, in soil amended with diammonium phosphate (DAP) and not previously grazed. Management practices also influenced soil chemical parameters, especially organic matter content and total N, which were 10% and 55% higher under NT than under MP. Aggregate stability was 61% higher in NT than in MP, 15% higher in P-fertilized soil, and also 9% higher in not grazed strips, bulk density being 12% lower in NT systems compared with MP. DAP application and the absence of grazing also reduced bulk density (3%). Using conservation tillage systems, fertilizing crops with DAP, and avoiding grazing contribute to soil health preservation and enhanced crop production.     

Effects of a waste paper product on soil phosphorus, carbon, and bulk density

Applications of manures to agricultural fields have increased soil test values for P to high levels in parts of the USA and thus increased the likelihood that P will be transported to surface water and degrade its quality. Waste paper applications to soils with high STP (soil test P) may decrease the risk of P transport to surface water by decreasing DRP (dissolved reactive P) by the formation of insoluble Al-P complexes and providing organic matter to improve infiltration. A field experiment was conducted near Booneville, AR (USA) to assess the effects of different rates of a waste paper product addition on STP, soil bulk density, and total soil C with a soil with approximately 45 mg Bray1-P kg-1 soil (dry weight). A Leadvale silt loam soil (fine-silty, siliceous, thermic Typic Fragiudult) was amended with 0, 22, 44, or 88 Mg waste paper product ha-1 to supply approximately 90, 170, or 350 kg Al ha-1, respectively. One year after additions, there was a strong negative correlation between waste paper product application rates and soil bulk density, and a strong positive correlation between rates and total soil C content. Soil bulk density and total C 2 yr after additions, and soil DRP and Bray1-P were not affected by waste paper additions. These results support the hypothesis that decreases in DRP in runoff from soils receiving waste paper additions were probably due to changes in soil organic matter and bulk density, rather than changes in the chemical forms of soil P.     

MODELING RUNOFF FROM VARIABLE SOURCE AREAS IN HUMID,SHALLOW WATER TABLE ENVIRONMENTS1

ABSTRACT: Variable Source Areas (VSAs) are zones with water saturated soils in forested wetlands fringing streams and creeks. Runoff from these areas is generated by saturation excess after a shallow water table rises and inundates the ground surface. In humid regions, like Florida and the Southeast, VSAs are believed to produce most of the runoff in shallow water table environments. Modeling the spatial extent and temporal fluctuation of a VSA is difficult because the formation of a VSA depends on a number of hydrological and morphological factors like rainfall intensity, soil texture, water table depth, and topographic attributes of the terrain. In this paper, we couple a digital elevation model with a two‐dimensional variable saturation model to illustrate the formation of a VSA at the hillside scale. The topography derived from the digital elevation model forms the upper domain geometry for the two‐dimensional finite element simulations of variable saturated flow. The objectives are: (1) to model the spatial and dynamic fluctuation of a VSA, and (2) to understand the roles of rainfall variability and terrain attributes on the formation of a VSA. Results show that hillsides with shallow water table depths, low saturated hydraulic conductivity, mild slopes, and concave slope curvature were more susceptible to runoff from a variable source. Runoff from a variable source showed little sensitivity to rainfall intensity. In general, landscapes with steep slopes generated a small VSA and a seepage face that vanished rapidly with time. In contrast, flat terrains are more amenable to VSA and retain ground surface inundation for longer periods of time.     

Developing risk models of Cryptosporidium transport in soils from vegetated, tilted soilbox experiments

Transport of Cryptosporidium parvum through macroporous soils is poorly understood yet critical for assessing the risk of groundwater contamination. We developed a conceptual model of the physics of flow and transport in packed, tilted, and vegetated soilboxes during and immediately after a simulated rainfall event and applied it to 54 experiments implemented with different soils, slopes, and rainfall rates. Using a parsimonious inverse modeling procedure, we show that a significant amount of subsurface outflow from the soilboxes is due to macropore flow. The effective hydraulic properties of the macropore space were obtained by calibration of a simple two-domain flow and transport model that accounts for coupled flow in the matrix and in the macropores of the soils. Using linear mixed-effects analysis, macropore hydraulic properties and oocyst attenuation were shown to be associated with soil bulk density and rainfall rate. Macropore flow was shown to be responsible for bromide and C. parvum transport through the soil into the underlying pore space observed during the 4-h experiments. We confirmed this finding by conducting a pair of saturated soil column studies under homogeneously repacked conditions with no macropores in which no C. parvum transport was observed in the effluent. The linear mixed-effects and logistic regression models developed from the soilbox experiments provide a basis for estimating macropore hydraulic properties and the risk of C. parvum transport through shallow soils from bulk density, precipitation, and total shallow subsurface flow rate. The risk assessment is consistent with the reported occurrence of oocysts in springs or groundwater from fractured or karstic rocks protected only by shallow overlying soils.     

INFLUENCES OF SOIL DATASET RESOLUTION ON HYDROLOGIC MODELING1

ABSTRACT: Soils represent a fundamental abiotic parameter in defining the characteristics of an ecosystem. The U.S. Department of Agriculture (USDA) Natural Resources Conservation Service (NRCS) produces the most detailed digital spatial soil datasets that are publicly available. The Soil Survey Geographic (SSURGO) database contains basic attributes for the continuous coverage of soils across the United States. In its standard format, the SSURGO database is incompatible for use within the ArcView Soil and Water Assessment Tool (SWAT). A modified version of the State Soil and Geographic (STATSGO) database is the template soils dataset used by ArcView SWAT. This paper presents the methodology and development of a SSURGO database preprocessor extension for the ArcView SWAT model. A case study for the Upper Sabinal River Watershed near Uvalde, Texas, is given. Results indicate that hydro‐logic output parameter differences occur when comparing the STATSGO and SSURGO database information in the ArcView SWAT model under identical modeling conditions. Specifically, the SSURGO model produced a greater daily mean water yield with evapotranspiration and surface runoff being found consistently lower across the watershed. The most likely causes assigned to this phenomenon were higher percolation and resulting ground water return flow values due to significantly larger saturated hydraulic conductivity values associated with the SSURGO 2.x database.     

USING SOIL TEXTURE TO ESTIMATE SATURATED HYDRAULIC CONDUCTIVITY AND THE IMPACT ON RAINFALL-RUNOFF SIMULATIONS1

ABSTRACT: In this paper a new set of soil texture data is used to estimate the spatial distribution of saturated hydraulic conductivity values for a small rangeland catchment. The estimates of conductivity are used to re-excite and re-evaluate a quasi-physically based rainfall-runoff model. The performance of the model is significantly reduced with conductivity estimates gleaned from soil texture data rather than the infiltration data used in our previous efforts.     

Soil chemical changes resulting from irrigation with water co-produced with coalbed natural gas

Land application of coalbed natural gas (CBNG) co-produced water is a popular management option within northwestern Powder River Basin (PRB) of Wyoming. This study evaluated the impacts of land application of CBNG waters on soil chemical properties at five sites. Soil samples were collected from different depths (0-5, 5-15, 15-30, 30-60, 60-90, and 90-120 cm) from sites that were irrigated with CBNG water for 2 to 3 yr and control sites. Chemical properties of CBNG water used for irrigation on the study sites indicate that electrical conductivity of CBNG water (EC(w)) and sodium adsorption ratio of CBNG water (SAR(w)) values were greater than those recommended for irrigation use on the soils at the study sites. Soil chemical analyses indicated that electrical conductivity of soil saturated paste extracts (EC(e)) and sodium adsorption ratio of soil saturated paste extracts (SAR(e)) values for irrigated sites were significantly greater (P < 0.05) than control plots in the upper 30-cm soil depths. Mass balance calculations suggested that there has been significant buildup of Na in irrigated soils due to CBNG irrigation water as well as Na mobilization within the soil profiles. Results indicate that irrigation with CBNG water significantly impacts certain soil properties, particularly if amendments are not properly utilized. This study provides information for better understanding changes in soil properties due to land application of CBNG water. These changes must be considered in developing possible criteria for preserving fragile PRB ecosystems.