VZJ Journal of Natural Resources and Life Sciences Education
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Tokunaga, T. K.
Right arrow Articles by Wan, J.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Tokunaga, T. K.
Right arrow Articles by Wan, J.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Tokunaga, T. K.
Right arrow Articles by Wan, J.
Related Collections
Right arrow Capillary Barriers
Right arrow Water Retention/Capillary Pressure
Right arrow Macroporous/Aggregated Media
Published in Vadose Zone Journal 2:322-329 (2003)
© 2003 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA

ORIGINAL RESEARCH PAPER

Moisture Characteristics of Hanford Gravels

Bulk, Grain-Surface, and Intragranular Components

Tetsu K. Tokunaga*, Keith R. Olson and Jiamin Wan

Lawrence Berkeley National Laboratory, Berkeley, CA 94720
* Corresponding author (tktokunaga{at}lbl.gov).

Received 2 May 2003.



    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Although gravels comprise large portions of some vadose zones, their unsaturated hydraulic properties have received relatively little attention. This study examines moisture retention relations in the 2- and 6-mm size fractions of gravels from the Hanford formation vadose zone (Washington State). Understanding flow and transport within this formation is important because parts of it have become contaminated by leakage of radioactive wastes at the Hanford Site. Moisture retention relations were obtained for a very wide energy range, with attention to water retained in intragranular pores and along grain surfaces. External surfaces of these gravels have root mean-squared roughnesses (rmsr) in the micrometer range, with sparsely distributed deep (hundreds of micrometers) pits. Water films on these external surfaces are volumetrically insignificant at matric potentials less than about -2 kPa. Residual water in these gravels occurs in intragranular pores, accounts for about 10% of the total porosity, and is effectively hydraulically immobile. The intragranular domain in Hanford gravels also has a large specific surface area of about 11 m2 g-1. Thus, exchanges of solutes (including contaminants) between the intragranular domain of Hanford gravels and their immediate surrounding are significant and diffusion limited.

Abbreviations: AFM, atomic force microscope • BET, Brunauer–Emmett–Teller • RH, relative humidity • rmsr, root mean-squared roughness


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
GRAVELS CAN MAKE UP large fractions of the subsurface, including vadose zones. The extensiveness of some gravel deposits is evident from the fact that they serve as major aquifers that supply groundwater for agricultural, industrial, and municipal use. Some gravel deposits are mined for use in construction. Because some gravel deposits do occur in the vadose zone, understanding unsaturated flow and transport in such settings requires knowledge of their hydraulic properties. Some gravel deposits occur in heavily contaminated vadose zones such as the Hanford Site (Washington State), where radioactive wastes have leaked. Gravels are also an important component in engineered capillary barriers for subsurface waste isolation (Ross, 1990; Selker, 1997). However, relatively little information is available on the unsaturated hydraulic properties of gravels. Measurements of moisture characteristics on gravels have been reported by Miyazuki (1988), Khaleel and Relyea (1997), Conca et al. (1998), and Tokunaga et al. (2002). Unsaturated hydraulic conductivities of gravels have been reported by Miller and Bunger (1963), Conca et al. (1998), and Khaleel and Relyea (2001).

Because of their large size, water retention in unsaturated gravels must be mainly associated with grains rather than with the macroscopic pore network. Thus, water retention depends on intragranular and grain surface morphology. Some previous studies have examined aspects of water retention in gravel, especially as it relates to availability for plant uptake (e.g., Coile, 1953; Flint and Childs, 1984; Montagne et al., 1992; Brouwer and Anderson, 2000). We suggest several possible basic grain morphologies depicted in Fig. 1 for cross sections through rounded grains. Analogous varieties exist for angular grains. The simplest representation is that of a smooth-surfaced ellipsoid of revolution, with negligible intragranular porosity. Other possibilities include rough-surfaced grains and varying importance of internal porosity, as well as possible additional weathering-induced near-surface alterations.



View larger version (72K):
[in this window]
[in a new window]
 
Fig. 1. Some possible combinations of intragranular porosity and grain surface roughness, illustrated as cross sections through grains.

 
Whether crystalline or amorphous, the larger grain size of gravels offers the potential for significant intragranular porosity and intragranular surface area. Diffusion coefficients within this intragranular domain can be very low because of low porosity and constrictive and tortuous pathways of microporous networks of boundaries between microcrystals. Diffusive transport times within grains scale with the square of the diameter, such that equilibration times within 4-mm gravel will be 400 times longer than in 200-µm fine sand of otherwise equal nature. Potassium ion diffusion coefficients ranging from 2.2 x 10-14 to 1.7 x 10-11 m2 s-1 were measured in a variety of gravels by Conca and Wright (1990). Thus, porous gravels can behave as sources and sinks of solutes, including contaminants, with slow exchange times. Water retained in intragranular pores will likewise exchange with the surrounding environment only over long times, thus behaving hydraulically as residual moisture. The distribution coefficient (Kd) values measured by Kaplan et al. (2000) for Sr2+ equilibration provided strong evidence for significant amounts of intragranular surface in Hanford gravels. The Kd values in their gravels were about 50% that of the <2-mm size fraction, hence still quite significant, and much greater than expected from external surface area alone.

The external surfaces of gravels can also be important since their unsaturated hydraulic properties can be strongly influenced by water films coating grain surfaces. Grain surface wettability and microtopography are then important in controlling water film hydraulic relations. The various possible configurations illustrated in Fig. 1 suggest that the unsaturated hydraulic properties of gravels can be quite diverse. Measurements of intragranular and surface characteristics will help identify which of the characteristics illustrated in Fig. 1 apply to specific types of gravels.

In a recent study on Hanford formation gravels, large saturation changes occurring within the near-zero matric potential region were reported (Tokunaga et al., 2002). Since that study was restricted to a narrow energy range, an understanding of the nature of residual saturation remains to be developed. This was especially warranted because the levels of residual saturation were relatively high, in the range of 0.1 to 0.2. Thus, the present work addresses a much broader range of matric potentials and saturations in Hanford gravels. We used suction plate, pressure plate, and vapor pressure methods to characterize samples. Water vapor adsorption–desorption measurements were obtained to test the hypothesis that the intragranular porosity of these gravels is substantial.

In addition to these bulk measurements, external grain surface properties were characterized microscopically (laser profilometry, atomic force microscopy, and scanning electron microscopy) and hydraulically (surface water film measurements at controlled matric potentials). We define the gravel grain surface moisture characteristic as the relation between the average (external surface) water film thickness and the matric potential (Tokunaga and Wan, 1997). The grain surface moisture characteristic is the surface analog of the familiar volumetric moisture characteristic relation. This relation is important to quantify because it can strongly influence unsaturated flow and transport in open framework (well-sorted, with minor fine fraction) gravels.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Gravels
Gravel samples were obtained from the Hanford Site in south-central Washington State. The Hanford formation gravels and sandy gravels were deposited during cataclysmic floods released from Pleistocene ice-dammed lakes. The less weathered gravel fraction retains most features of the parent rocks, predominantly basalt. Hydraulically important features imparted by the basalt include interconnected, very fine microporosity due to boundaries between individual mineral grains and much larger vesicles due to gas bubbles trapped during cooling of the rock. Influences of these features on moisture retention will be discussed below. Our samples were excavated from Site 218-E-12B, in the northeast corner of the 200-East Area. We report primarily on the 2.0- to 2.4- and 5.6- to 6.3-mm size fractions, referred to here as 2- and 6-mm Hanford gravels, respectively. As shown below, the Hanford gravels have significant intragranular porosity and surface area. To obtain comparisons with gravels that have less intragranular influences, additional samples were obtained from deposits along the Maurice River, New Jersey. These gravels are primarily quartz (98–99%). The 2.0- to 2.4- and 8.0- to 9.5-mm size fractions of these gravels were only used in a few water vapor sorption measurements and will be referred to as 2- and 9-mm quartz gravel, respectively. The overall shapes of the Hanford and quartz gravels were similar, with both having sphericities of about 0.7 and roundness of about 0.6 (Krumbein and Sloss, 1963).

Intragranular Properties
Water vapor adsorption and desorption isotherms were measured on the Hanford and quartz gravels to characterize moisture retention at low (very negative) water potentials and to obtain estimates of specific surface area. The oven-dry (110°C, 24 h) Hanford sample masses were about 34 and 66 g for the 2- and 6-mm fractions, respectively. The 2- and 9-mm fractions of the quartz gravel were tested in nominally 80-g samples. The initial motivation for testing two different size fractions was based on the assumption that the measurements would be strongly correlated to external surface area; hence, water retention would be inversely proportional to the grain size. Large sample masses were used to partially compensate for the lower specific external surface area. After oven drying and cooling to room temperature, samples were placed into weighing jars and equilibrated over (aqueous) saturated salt solutions (Table 1) kept in a controlled temperature (25 ± 0.5°C) incubator. These salt solutions provided controlled relative humidities (RH) ranging from 0.07 to 0.97 and water potentials from -364 to -4.6 MPa. The RH is the ratio of the water vapor pressure to the saturation vapor pressure (3.17 kPa at 25°C). The quartz gravels were tested only up to a RH = 0.33 because only adsorption in the low RH range can be used to infer surface areas and because the moisture characteristics of these gravels were beyond the scope of this study. The jars containing gravel samples were periodically removed from constant humidity chambers and weighed to monitor progression toward adsorption equilibrium. Vapor adsorption equilibrium was reached in 10 to 20 d, after which gravels were again oven dried. For desorption isotherms, gravel samples were immersed in water overnight, quickly drained, placed into weighing jars, and equilibrated at controlled RH with the same saturated salt solutions. Desorption equilibrium was reached in about 30 d for all samples except those in the chamber with RH = 0.967. These latter samples exhibited very slow drying because of weak vapor pressure gradients acting between the gravels and salt solution reservoir.


View this table:
[in this window]
[in a new window]
 
Table 1. Saturated salt solutions used in water vapor adsorption–desorption measurements.

 
Estimates of surface areas were obtained in two ways. In the first approach, water vapor adsorption data obtained with RH less than 0.30 were used in the Brunauer–Emmett–Teller (BET; Brunauer et al., 1938) method to calculate surface areas. This approach amounted to applying the standard BET analysis briefly described below, but with water as the adsorbate. Conventional surface area measurements were also obtained on these samples by BET analyses using Kr and N2 adsorbates (Micromeritics, Norcross, GA).

Bulk Gravel Moisture Characteristics
Moisture characteristics were only obtained for the Hanford gravel. The determination of moisture characteristic relations for gravels presents an experimental challenge because significant changes in saturation occur over matric head intervals that are shorter than typical column heights. At matric potentials close to zero, the upper portion of an equilibrium gravel column can be significantly less saturated than the bottom portion. Methods to obtain correct saturation–potential relations under these conditions have been presented by Dane et al. (1992), Liu and Dane (1995), Schroth et al. (1996), Jalbert and Dane (2001), and Tokunaga et al. (2002). The suction plate approach described in Tokunaga et al. (2002) was used here, with gravel packed to a height of 30 mm in large (82.6-mm i.d.), modified Tempe cells. Boundary matric potentials were controlled to within ± 15 Pa in this region of the drainage curve. Measurements at lower (more negative) matric potentials were obtained using the pressure plate method (Dane and Hopmans, 2002), with two-stage pressure regulation (final regulated potentials controlled to within ± 100 Pa), and with a 0.05-MPa (0.5 bar), high-flow plate and 0.5- MPa (5 bar) standard plate (Soilmoisture Equipment Corp., Goleta, CA). The pressure plate measurements were obtained down to a matric potential of -150 kPa. A filter paper interface was placed between the gravel and ceramic to facilitate hydraulic contact.

External Grain Surface Properties
To evaluate the role of external grain surfaces on water retention, measurements were made to characterize grain surface topography and to quantify the relation between average water film thickness and matric potential. Grain surface topography was characterized using a laser profilometer (UMB, Sunnyvale, CA) and an atomic force microscope (AFM; Park Scientific Autoprobe M5, Veeco Instruments, Woodbury, NY) as previously described (Tokunaga et al., 2000; Wan and Tokunaga, 2000). A scanning electron microscope was used to obtain images of gravel grain surfaces.

Grain surface moisture characteristics could in principle be obtained through very accurately determined adsorption isotherms and moisture characteristic measurements on bulk sediments, combined with information (or assumptions) on moisture within intragranular pores and in pendular rings. Instead, we chose to obtain average film thickness measurements directly at various points on a single grain surface. Microscopic inspection of external surfaces of 1- to 30-mm grain sizes of Hanford gravels did not reveal any size-dependent topography. Thus measurements on a single grain are expected to be representative of other grain sizes. Using a presaturated single grain provided the advantage of effectively constant grain matrix saturation, and avoiding the complicating influence of pendular rings. The method used to determine area-averaged film thickness is based on synchrotron X-ray fluorescence of a nonreactive tracer [selenium(VI) as SeO2-4, which does not undergo reduction because the system is kept equilibrated with atmospheric O2] dissolved in the aqueous phase, with the matric potential regulated by a modified suction plate device (Tokunaga et al., 2000). This method, previously tested on roughened quartz glass surfaces, was used on a 13-mm Hanford gravel specimen with its bottom surface cut flat to provide good hydraulic contact with the underlying porous plate. The X-ray microprobe beamline X26A at the National Synchrotron Light Source was used to obtain these measurements. Artifacts due to evaporation by local X-ray beam heating were insignificant (above our film thickness detection limit of 0.3 µm) because we used a monochromatic beam, short exposure times (5–7 s), and kept the sample in a humidified (pinhole-vented) chamber.


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Intragranular Properties
The 2- and 6-mm Hanford gravels had similar water vapor sorption–desorption isotherms, with both grain sizes retaining substantial mass fractions of water (Fig. 2). In contrast, the quartz gravels had low water contents, barely at the limit of quantification ({approx}0.0002 g g-1). These results provide evidence for significant intragranular porosity and intragranular surface area in Hanford gravels for two general reasons. First, the water contents are far too large to have been associated solely with films on external surfaces of solid (nonporous) grains and pendular rings in equilibrium with very low water potentials. A uniform water film coating a solid (nonporous) spherical grain of diameter D has a thickness f of

[1]
where w is the water/solid mass ratio, and {rho}s and {rho}w are the densities of the solid and water, respectively. This estimate leads to hypothetical film thicknesses ranging from 10 to 60 µm, for the matric potential range where water actually occurs in monomolecular to several molecular layer films (i.e., nanometer films). We show later that external surface roughness, although significant, cannot reconcile this approximately 104–fold discrepancy in hypothetical film thickness. Water retention in pendular rings (Fig. 3) also cannot account for much of this discrepancy since the energy range under consideration is associated with radii of curvature r < 5 nm, based on the Kelvin equation:

[2]
where M is the molecular weight of water (18.0 g mol-1), {sigma} is the air–water surface tension (7.3 x 10-2 N m-1), R is the gas constant (8.314 J K-1 mol-1), and T is the temperature (298 K). The factor of 2 appears in the numerator for the familiar case of hemispherical and spherical interfaces, but was removed in Eq. [2] since the concave radius of curvature bridging grains (r in Eq. [2] is designated r1 in Fig. 3) is much smaller than the convex radius of curvature that circumscribes the axis between grains (r2 in Fig. 3). Since the Laplace relation involves sums of inverse radii, the much larger radius has only minor influence. Whether or not the Kelvin equation is valid at r < 5 nm and whether or not the macroscopic value of {sigma} is appropriate at this scale are outside the scope of this study. However, it is worth noting that Fisher and Israelachvili (1979) demonstrated applicability of the Kelvin equation down to r = 4 nm for cyclohexane wetting mica (with its macroscopic {sigma}).



View larger version (39K):
[in this window]
[in a new window]
 
Fig. 2. Water vapor adsorption and desorption isotherms for 2- and 6-mm Hanford gravels, and for 2- and 9-mm quartz gravels.

 


View larger version (84K):
[in this window]
[in a new window]
 
Fig. 3. Geometry of pendular ring surface at a grain–grain contact. At low water contents r1 << r2, so that r-11 >> r-12.

 
The second reason why these macroscopic results suggest substantial intragranular water retention becomes evident on considering the fact that practically identical isotherms were obtained for the 2- and 6-mm Hanford gravels. If only external grain surfaces were responsible for water retention, then the 2-mm gravel should have three times greater water contents than the 6-mm gravel at any given potential. This scaling of surface area per unit volume, also implicit in Eq. [1], is not supported by the data.

Calculations of effective surface areas can lead to a better understanding of residual or immobile water in the Hanford gravels. Two approaches were taken, both relying on BET calculations (Brunauer et al., 1938). In the first approach, data from the water vapor sorption isotherms were expressed in terms of the BET variable:

[3]

When B is linear with respect to RH (up to RH ~0.3), multilayer adsorption conforms to the BET model. Such conformity was obtained with our measurements on Hanford gravels (Fig. 4). The value of w associated with monomolecular coverage is obtained from the inverse of the sum of the slope and intercept of B(RH). This monolayer w was combined with the effective area of water molecules, 0.106 nm2 per molecule (Gregg and Sing, 1982), to obtain the water-based specific surface area measurements of 27 and 22 m2 g-1 for the 2- and 6-mm gravels, respectively. Measurements obtained using less reactive gases (N2 and Kr) yielded specific surface areas of 11.7 and 11.3 m2 g-1 for the 2- and 6-mm gravels, respectively. It is also worth mentioning that N2 BET measurements of the <2-mm (median grain size ~500 µm) fraction from the Hanford formation yielded a value of 15.5 m2 g-1, that is, not much greater than the gravel fraction. All of these surface area values are large and in the range of values for kaolinite, a low surface area clay. Measurements using water as the adsorbate commonly do yield higher BET specific surface areas than N2- and Kr-based measurements, with differences attributed to solvation of adsorbed ions and hydration reactions (Quirk, 1955; Gregg and Sing, 1982). Differences in values obtained with different adsorbates further reflect the operationally dependent nature of specific surface area measurements (Sposito, 1984). Nevertheless, the collective evidence obtained from measurements of water uptake at low potentials indicates that the Hanford gravels have very significant internal surface area. Such high intragranular surface areas combined with large grain sizes are expected to cause strongly diffusion-limited water and solute exchange.



View larger version (35K):
[in this window]
[in a new window]
 
Fig. 4. Brunauer–Emmett–Teller (BET) plots for estimating gravel surface areas from water vapor adsorption. Surface areas determined by N2 and Kr adsorption are shown for comparison.

 
Bulk Gravel Moisture Characteristics
Drainage curves obtained on the 2- and 6-mm Hanford gravels are shown in Fig. 5. The near-zero matric potential region (Region c in Fig. 5) within which large decreases in saturation accompanied small energy decreases is gravity-corrected for variable saturation within individual equilibrium states (Tokunaga et al., 2002). The intermediate potential range results obtained with a pressure plate are shown in Region b of Fig. 5. Up to 15 d were allowed for individual equilibration steps. Even with this precaution and the use of a filter paper interface, contact impedance can make it difficult to achieve equilibrium in such coarse media. However, the fact that the moisture contents obtained with the pressure plate are consistent with measurements obtained with vapor adsorption support the reliability of both methods. The lowest energy range (Region a in Fig. 5) is based on the water vapor adsorption isotherms described above. The water potentials obtained in the regulated relative humidity chambers were assumed to be approximately equal to matric potentials for this purpose. This amounts to assuming that the water potential equals the matric potential plus an additive but insignificant solute potential. This assumption is justified in nearly all of these vapor pressure equilibrations since the highest water contents obtained by vapor adsorption (region a in Fig. 5) are about equal to the lowest water contents obtained with the pressure plate (Region b of Fig. 5). The solute potential at this water content (saturation {approx}0.1) has a magnitude of about -102 kPa (estimated from the van't Hoff relation and measured major ion concentrations in Hanford pore water). Thus, the solute potential contributes only about 2% to the highest water potential imposed by vapor pressure regulation (-4.6 x 103 kPa in Table 1). When the saturation is reduced to 0.05, the solute potential is about -2 x 102 kPa, or about 0.1% of the corresponding water potential of -2 x 105 kPa (Fig. 5). Three main features of Hanford gravel moisture characteristics are evident from these results. First, practically all of the saturation changes occur within 1 to 2 kPa of zero matric potential. Changes in saturation within this near-zero potential region are large enough that gravity corrections were required for the 6-mm gravel. Gravity corrections for the 2-mm gravel drainage curve are minor (Tokunaga et al., 2002). Second, the fairly high residual saturation region spans a very broad energy range and associated with intragranular porosity (as discussed in the previous section). Intragranular moisture retention within gravel fractions in soils has been reported previously (e.g., Coile, 1953; Flint and Childs, 1984; Montagne et al., 1992; Brouwer and Anderson, 2000), and on tuff gravel by Conca et al. (1998). The results presented in Fig. 5 show that moisture capacities (the change in saturation per unit change in matric potential) for these gravels are very low. At most, there is only about 5% change in saturation for the interval from -1 to -104 kPa. Third, within this broad residual saturation region, the 2- and 6-mm gravel curves are practically the same. This too results from water being retained primarily in intragranular pores, so that retention scales with bulk density rather than with external grain surface area. The nearly identical residual saturations for the two grain sizes also indicate that even if weathering rinds are present they have no measurable impact on moisture retention. These results support the morphologies shown in Fig. 1b and 1e for the Hanford gravels.



View larger version (31K):
[in this window]
[in a new window]
 
Fig. 5. Moisture retention in 2- and 6-mm Hanford gravels, for a wide range of matric (water) potentials. Measurement methods used are indicated along the x axis as vapor pressure regulation (a), pressure plate (b), and suction plate (c). The vapor pressure equilibration region data were obtained in adsorption mode. The suction plate and pressure plate data were obtained in drainage mode.

 
Grain Surface Properties
The microtopography of grain surfaces reflect the combined influences of parent rock structure and weathering. In the case of Hanford gravels, many of these grains are basalt in origin and have been mechanically weathered through glacial fluvial processes. Thus, the vesicular imprint from trapped gas bubbles in the basalt remains, as shown in the photographs of the external surface and intragranular domain (Fig. 6). Together with the previously described results, we conclude that the Hanford gravels are represented by morphology depicted in Fig. 1e. Scanning electron micrographs of the external grain surface qualitatively show surface roughness over a range of scales (Fig. 7). Note that the surface pits, originally vesicles in the bulk basalt rock, are isolated and substantially deeper than the more pervasive roughness features. Isolated pits are typically in the range of 100 to 400 µm in depth, whereas the majority of the external surface roughness is about two orders of magnitude finer. The latter was quantified at intermediate scale (several hundred micrometer traverses, with 0.1-µm vertical resolution) through laser profilometry, with typical profiles shown in Fig. 8. Note that the vertical scale on this figure is expanded by a factor of 3. The rmsr at this scale is about 2 µm. Atomic force microscope scans were difficult to obtain because large-scale topography sometimes exceeded the instruments working range (8 µm). Nevertheless, some small areas (45–900 µm2) were scanned successfully, with typical rmsr values in the range of 0.19 to 1.92 µm. The average rmsr obtained from 22 area scans was 0.83 µm. Typical AFM surface profiles are shown in Fig. 9. The greater rmsr value obtained on the laser profiles probably resulted from scale dependence of surface roughness, increasing with measurement length or area. Brown and Scholz (1985) found that rock fracture surfaces exhibit scale-dependent roughness. In the case of granular materials, the overall curvature of the surface itself imparts greater apparent roughness with increased measurement scale, even on an ideally smooth surface, since greater elevation variations are encompassed. The AFM average rmsr of 0.83 µm is a good indicator of the local topographic constraint on water film thickness since it was obtained for distances ranging from 7 to 30 µm. Greater distances typically intercept major ridges and pits that support only very thin (nanometer) and very thick (up to hundreds of micrometers) films, respectively.



View larger version (96K):
[in this window]
[in a new window]
 
Fig. 6. Photographs of Hanford formation basalt gravel (a) external surface, and (b) cross section ({approx}300-µm-thick, polished slice). Note that vesicles appear as isolated circular voids in the interior region and are pits when exposed at external surfaces.

 


View larger version (99K):
[in this window]
[in a new window]
 
Fig. 7. Scanning electron micrographs of a Hanford gravel surface.

 


View larger version (24K):
[in this window]
[in a new window]
 
Fig. 8. Laser profiles of Hanford gravel grain surfaces (rmsr = root mean-squared roughness).

 


View larger version (33K):
[in this window]
[in a new window]
 
Fig. 9. Representative atomic force microscope (AFM) images of Hanford gravel grain surfaces.

 
Grain surface moisture characteristics were obtained at eight randomly selected spots on the gravel sample for the matric potential range of -0.1 to -6.1 kPa (-10 to -620 mm H2O). Representative drainage and wetting curves are shown in Fig. 10. Although there is some apparent variability, three main features are evident. First, very high values of average film thickness (7–10 µm) were obtained at very near-zero (-0.05 to -0.10 kPa) matric potentials. These high values of average film thickness result from filling isolated pits in the near-zero potential limit. Such pits are practically empty at lower potentials because of their large cross sections and shallow depths. Second, average film thicknesses are equal to or less than the quantification limit of 0.3 µm at matric potentials less than about -2 kPa. This result supports the argument that moisture retention at even lower energies (pressure plate and vapor pressure controlled measurements) largely occurs in the intragranular domain. Even when the film thickness is 1 µm, the gravimetric water content associated with such a coating on the exterior of a 2-mm-diameter grain is only about 10-3 g g-1 (from Eq. [1]). The third feature shown in these results is that surface film moisture characteristics do not exhibit hysteresis. This lack of hysteresis in film characteristics was demonstrated previously on a fractured rock surface (Tokunaga and Wan, 1997).



View larger version (23K):
[in this window]
[in a new window]
 
Fig. 10. Synchrotron X-ray fluorescence microprobe measurements of draining and wetting film moisture characteristics at three different locations (a, b, c). The spot size is about 500 by 500 µm.

 
Qualitative aspects of film flow along external surfaces of Hanford gravel can be deduced from these measurements. Film transmissivity, the unit gradient flux per surface length transverse to the flow direction, depends on surface topography and matric potential (Tokunaga and Wan, 1997). Since the vesicular pits are sparsely distributed rather than interconnected (Fig. 6 and 7), these deep surface features do not contribute significantly to film flow. Recall that the grain surfaces have been smoothed by glacial and fluvial weathering, resulting in rmsr values in the range of 1 µm. Under finite matric potentials, such smooth surfaces cannot support thick (e.g., >2 µm) water films. Therefore, fast, unsaturated flow is not possible in these gravels below the matric potential at which pendular rings are no longer interconnected, since water films then control flow (Wan and Tokunaga, 1997).


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Various techniques were used to quantify moisture characteristics of the 2- and 6-mm fractions of Hanford gravels. Differences between moisture characteristics for these two grain-size fractions were only important in the near-zero matric potential region associated with drainage of the primary pore network. Across a broad range of matric potentials, water content remains practically unchanged and nearly identical for the two grain sizes. Water retention in this range is associated with intragranular pores and amounts to about 10% of the total porosity. The specific surface area of these gravels is fairly high (about 11 m2 g-1) and largely associated with the intragranular domain. The combination of large grain size, significant intragranular porosity, and significant intragranular surface area make these gravels sinks and sources for solutes, with potentially long, diffusion-limited equilibration times. External grain surfaces have rmsr of about 1 µm for most locations, with isolated deep pits. Such surface topography does not support thick water films at matric potentials less than about -1 kPa and therefore will not support fast film flow at lower energies. Although measurements were done on specific size fractions of Hanford gravels, the conclusions concerning the intragranular domain and grain-surface characteristics remain valid for these gravels in well-graded mixtures of grain size.


    ACKNOWLEDGMENTS
 
We thank Andrew Mei, Eduardo Saiz, and Robert Conners of LBNL, and Tony Lanzirotti of University of Chicago for technical support, and John Zachara, Robert Lenhard, Steve Smith, and Bruce Bjornstad of PNNL for samples of Hanford formation sediment. We also thank Sam Ricci, Jr., and Ricci Brothers Sand Company, Port Norris, NJ for providing quartz gravel samples. Helpful internal review comments by Tim Kneafsey (LBNL), and final constructive comments by two anonymous reviewers are gratefully acknowledged. This work was carried out under U.S. Dep. of Energy (DOE) Contract no. DE-AC03-76SF00098, with funding provided by the DOE, Environmental Management Science Program. Research carried out (in part) at the National Synchrotron Light Source, Brookhaven National Laboratory, which is supported by the U.S. Department of Energy, Division of Materials Sciences and Division of Chemical Sciences under Contract no. DE-AC02-98CH10886.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 




This article has been cited by other articles:


Home page
Soil Sci.Home page
M. C. Duniway, J. E. Herrick, and H. C. Monger
The High Water-Holding Capacity of Petrocalcic Horizons
Soil Sci. Soc. Am. J., April 5, 2007; 71(3): 812 - 819.
[Abstract] [Full Text] [PDF]


Home page
Vadose Zone JHome page
K. Hansson and L.-C. Lundin
Water Content Reflectometer Application to Construction Materials and its Relation to Time Domain Reflectometry
Vadose Zone J., March 8, 2006; 5(1): 459 - 468.
[Abstract] [Full Text] [PDF]


Home page
Vadose Zone JHome page
Q. Hu, T. J. Kneafsey, J. J. Roberts, L. Tomutsa, and J. S. Y. Wang
Characterizing Unsaturated Diffusion in Porous Tuff Gravel
Vadose Zone J., November 1, 2004; 3(4): 1425 - 1438.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Tokunaga, T. K.
Right arrow Articles by Wan, J.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Tokunaga, T. K.
Right arrow Articles by Wan, J.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Tokunaga, T. K.
Right arrow Articles by Wan, J.
Related Collections
Right arrow Capillary Barriers
Right arrow Water Retention/Capillary Pressure
Right arrow Macroporous/Aggregated Media


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
The SCI Journals Agronomy Journal Crop Science
Journal of Natural Resources
and Life Sciences Education
Soil Science Society of America Journal
Journal of Plant Registrations Journal of
Environmental Quality
The Plant Genome