Comparisons between acoustic measurements and predictions of the bedload transport of marine gravels

Continuous, detailed records of marine gravel transport have been obtained acoustically and compared with bedload transport rates (q_b) predicted by five bedload transport equations using measurements of the near-bed turbulent current flow. When mean flow data are used in these equations, total q_b estimates are similar to those measured. However, when instantaneous flow data are used, total q_b is over-estimated by approximately one order-of-magnitude. Based on the acoustic measurements, an empirical equation has been obtained that gives accurate estimates of total q_b over a tidal cycle and simulates well the intermittent characteristics of marine bedload transport.     

The partitioning of the total sediment load of a river into suspended load and bedload: a review of empirical data

The partitioning of the total sediment load of a river into suspended load and bedload is an important problem in fluvial geomorphology, sedimentation engineering and sedimentology. Bedload transport rates are notoriously hard to measure and, at many sites, only suspended load data are available. Often the bedload fraction is estimated with ‘rule of thumb’ methods such as Maddock’s Table, which are inadequately field‐tested. Here, the partitioning of sediment load for the Pitzbach is discussed, an Austrian mountain stream for which high temporal resolution data on both bedload and suspended load are available. The available data show large scatter on all scales. The fraction of the total load transported in suspension may vary between zero and one at the Pitzbach, while its average decreases with rising discharge (i.e. bedload transport is more important during floods). Existing data on short‐term and long‐term partitioning is reviewed and an empirical equation to estimate bedload transport rates from measured suspended load transport rates is suggested. The partitioning averaged over a flood can vary strongly from event to event. Similar variations may occur in the year‐to‐year averages. Using published simultaneous short‐term field measurements of bedload and suspended load transport rates, Maddock’s Table is reviewed and updated. Long‐term average partitioning could be a function of the catchment geology, the fraction of the catchment covered by glaciers and the extent of forest, but the available data are insufficient to draw final conclusions. At a given drainage area, scatter is large, but the data show a minimal fraction of sediment transported in suspended load, which increases with increasing drainage area and with decreasing rock strength for gravel‐bed rivers, whereby in large catchments the bedload fraction is insignificant at ca 1%. For sand‐bed rivers, the bedload fraction may be substantial (30% to 50%) even for large catchments. However, available data are scarce and of varying quality. Long‐term partitioning varies widely among catchments and the available data are currently not sufficient to discriminate control parameters effectively.     

Modeling Bedload Transport of Juvenile Bivalves: Predicted Changes in Distribution and Scale of Postlarval Dispersal

The potential consequences of bedload transport of postlarvae for patterns of distribution of marine invertebrates were explored by developing a bedload transport model for juvenile bivalves in a small estuary in New Jersey, USA. A simple numerical model of tidal current hydrodynamics was developed based on field measurements of shear stresses near the bottom. Burrowing behavior of bivalves was incorporated into the model of bedload transport by using estimates of entrainment rates of Gemma gemma and Mya arenaria in a laboratory flume, and jump lengths of the bivalves were estimated by methods previously developed for noncohesive particles. Based on the flood domination and strong gradient of shear stresses in the Navesink estuary, our model predicted that juvenile bivalves would accumulate in the center of the estuary, traveling up to several kilometers over 30 days. Field distributions of juvenile bivalves were consistent with the model predictions for other species of bivalves but not for G. gemma, for which field distributions of both <500- and >500-μm individuals were concentrated in the eastern end of the estuary. Differences between the bedload model and G. gemma distributions suggest that spatial variation in burrowing behavior or biological interactions are playing an important role in maintaining distribution patterns of this species in spite of high levels of bedload transport. This modeling approach is applicable to other juvenile benthic invertebrates that disperse as bedload and is a useful model against which to compare field observations of rates of transport and patterns of distribution and abundance.     

Guidelines for sampling bedload transport with minimum uncertainty

The stability of river channels and their suitability as habitat for aqueous organisms is strongly controlled by the rate of bedload transport. Quantification of bedload transport rates in rivers is difficult, not only because of the temporal variation in transport, but also because of the cross‐channel variation in transport. The objectives of this study were: (i) to determine the effect of cross‐channel variation in bedload transport on the uncertainty of width‐integrated transport rates; and to use this knowledge (ii) to improve guidelines for bedload sampling. This was done through a thorough statistical evaluation of stochastic and systematic uncertainties involved in bed‐load transport measurements. Based on this evaluation, new guidelines are presented for determination of the number of samples and sampling positions across the channel that are required for bedload measurements in several types of sand‐bed rivers and gravel‐bed rivers. The guidelines relate to bedload measurements made with pressure‐difference (Helley‐Smith type) samplers that require numerous bedload samples of short duration at several positions across the channel. The results show that generally more sampling positions across the channel are required in gravel‐bed rivers than in sand‐bed rivers. For gravel‐bed rivers with unknown cross‐channel distribution of transport, at least 10 sampling positions are recommended, whereas for most sand‐bed rivers five positions suffice. In addition, at least 12 short‐duration samples are required at each position to obtain bedload estimates with uncertainties below 20%. If the same level of uncertainty is desired in the case of high spatial and temporal variation in transport rates, the number of short‐duration samples needed per sampling position increases to 40.     

Calibration of piezoelectric bedload impact sensors in the Pitzbach mountain stream

Non-invasive techniques to continuously measure the rates or intensities of bedload transport in gravel bed rivers would be quite useful for both research and practice. Since 1986 measurements have been made with piezoelectric bedload impact sensors (PBIS) in several mountain torrents. In 1994 and 1995, a calibration campaign was carried out at a water intake in the Pitzbach mountain stream in Austria, where weighing cells and a water pressure sensor in the settling basin allow sediment transport loads to be determined in 15 minute intervals. Daily averages of PBIS impulses correlate reasonably well with the measured sediment load. Using a linear calibration relation, the relative error between predicted and measured sediment loads is within a factor of about 1.5 to 2. The scatter between PBIS impulses and bedload volumes is much larger if the measurements are averaged over one hour intervals or when the original 15 minute values are considered. A comparison with measurements from the Swiss Erlenbach torrent and from laboratory experiments using sediment particles from the Erlenbach torrent indicates that the proportionality coefficient between PBIS impulses and sediment load appears to depend on site specific conditions, including factors such as flow intensity and sediment properties.     

Bedload transport of aggregated muddy alluvium from Cooper Creek, central Australia: a flume study

The ability of mud aggregates to form depositional bedforms is of considerable sedimentological importance for explaining the geomorphology of the Channel Country of central Australia as well as for understanding the depositional environment of certain argillaceous fluvial sequences in the rock record. The sediment transport and bedform development of mud aggregates from the floodplain of Cooper Creek, central Australia, was examined in a laboratory flume over a range of flow conditions. The aggregates were found to be clay-rich (>60% clay), nonsaline (<0·02%), fine sand-sized (mean d₅₀=0·13 mm), low density (2300 kg m^?3) and water-stable. Three wetting rates were applied to the sediment in the laboratory prior to wet sieving to replicate various field conditions and results in three mean aggregate sizes. Immersion wetting (no tension) represents inundation of the sediment by overland flow and results in aggregates of 0·13 mm. Tension wetting at 20 and 50 mm corresponds to high- and low-intensity rainfall and results in mean d₅₀ sizes of 0·75 and 0·70 mm, respectively. Immersion wetting is the most applicable wetting mode for hydraulic transport of aggregated sediment on the Cooper Creek floodplain. Considerable variability in sediment transport rates in the field could result from differences in pre-wetting of the aggregated sediment. The dominance of smectite in the clay mineralogy of the sediment is an important factor in the development of the aggregates; disaggregated sediment reaggregated in a laboratory after 2–3 wetting/drying cycles. In flume experiments, bedforms of aggregated mud ranging from lower-regime plane beds to upper-regime antidunes were observed. The aggregates moved predominantly as bedload with measured peak bedload concentrations being high compared with other flume studies. The highly mobile nature of this sediment in the field is due to the ready entrainment of low-density aggregates in the form of self-mulching vertisols across extensive floodplains. The occurrence of low-sinuosity braid-like channels on this extensive low-gradient semi-arid floodplain can be attributed to: (a) the passage of floodwaters across a floodplain with steeper gradients than adjacent more sinuous anastomosing channels; (b) the highly mobile nature of the low-density sediment aggregates; (c) the ability of the aggregates to be transported as bedload; and (d) their durable nature during transport.     

On Bagnold's sediment transport equation in tidal marine environments and the practical definition of bedload

Bagnold's sediment transport equation has proved to be important in studying tidal marine environments. This paper discusses three problems concerning Bagnold's transport equation and its practical application:

• 1 Bagnold's suspended-load transport equation and the total-load transport equation with are incorrect from the viewpoint of energy conservation. In these equations the energy loss due to bedload transport has been counted twice. The correct form should be for suspended-load transport and for total-load transport with
• 2 The commonly used Bagnold's transport coefficient K varies as a non-linear function of the dimensionless excess shear stress, which can be represented best by the power law , where the coefficient A and exponent B depend on sediment grain size D. The empirical values of A and B for fine to medium grained sands are determined using Guy et al.'s (1966) flume-experiment data.
• 3 The sediment transport rates predicted from this equation are compared with bedform migration measurements in the flume and the field. This comparison shows that the sediment transport rates measured from bedform migrations are higher than the predicted bedload transport rates, but comparable to the calculated total-load (bedload plus intermittent suspended-load) transport rates. This indicates that bedform migration involves both bedload and intermittent suspended-load transport. As a logical conclusion, bedform migration data should be compared with Bagnold's total-load transport equation rather than with his bedload transport equation. In this respect the term ‘bed material’ might be more appropriate than the term ‘bedload’ for estimating sediment transport rate from bedform migration data.

The sediment transport rates predicted from this modified Bagnold transport equation are in good agreement with field measurements of bedform migration rates in four individual tidal marine environments, which cover a wide range of sediment grain size, flow velocity and bedform conditions (ranging from small ripples, megaripples to sandwaves).     

The dynamics of a river bend: a study in flow and sedimentary processes

Comprehensive field measurements of flow and sedimentary processes have been made with the aid of stable scaffolding bridges spaced along the length of a bend of the River South Esk, Scotland. At river stages between about two-thirds full and bankfull, channel width, mean depth and mean flow velocity at a cross-section vary little in the streamwise direction. Flow resistance reaches a maximum at these stages, and the bed topography is stable and in equilibrium with flow and bedload transport. Stable flow geometry is thus related in some way to energy conservation, and to maximization of flow resistance. Detailed observations over a large range of river stages of mean velocity distributions, secondary circulation, water surface configuration, bed shear stress and resistance to flow, bed configurations and bed load transport rates agree with much (but not all) of the comparable published experimental studies and selected theoretical work. Generalized physical models of flow and sediment transport in natural curved channels (Engelund, 1974; Bridge, 1977) are demonstrated to be sound in basis and can simulate the bend studied very well. Although there is a pressing need for further development of these models, the results lend confidence to their use in simulating ancient river sedimentation. Sediment deposited on point bars is the result mainly of bedload transport over a range of near-bankfull stages. The areal distribution of grain-size characteristics and bed configurations at these stages give rise, with lateral deposition, to vertical facies sequences that vary substantially in the streamwise direction.     

Direct numerical simulation of bedload transport using a local, dynamic boundary condition

ABSTRACT Temporally and spatially averaged models of bedload transport are inadequate to describe the highly variable nature of particle motion at low transport stages. The primary sources of this variability are the resisting forces to downstream motion resulting from the geometrical relation (pocket friction angle) of a bed grain to the grains that it rests upon, variability of the near‐bed turbulent velocity field and the local modification of this velocity field by upstream, protruding grains. A model of bedload transport is presented that captures these sources of variability by directly integrating the equations of motion of each particle of a simulated mixed grain‐size sediment bed. Experimental data from the velocity field downstream and below the tops of upstream, protruding grains are presented. From these data, an empirical relation for the velocity modification resulting from upstream grains is provided to the bedload model. The temporal variability of near‐bed turbulence is provided by a measured near‐bed time series of velocity over a gravel bed. The distribution of pocket friction angles results as a consequence of directly calculating the initiation and cessation of motion of each particle as a result of the combination of fluid forcing and interaction with other particles. Calculations of bedload flux in a uniform boundary and simulated pocket friction angles agree favourably with previous studies.     

Geometry and kinematics of dunes during steady and unsteady flows in the Calamus River, Nebraska, USA

The geometry and kinematics of river dunes were studied in a reach of the Calamus River, Nebraska. During day-long surveys, dune height, length, steepness, migration rate, creation and destruction were measured concurrently with bedload transport rate, flow depth, flow velocity and bed shear stress. Within a survey, individual dune heights, lengths and migration rates were highly variable, associated with their three-dimensional geometry and changes in their shape through time. Notwithstanding this variability, there were discernible changes in mean dune height, length and migration rate in response to changing discharge over several days. Changes in mean dune height and length lagged only slightly behind changes in discharge. Therefore, during periods of both steady and unsteady flow, mean dune lengths were quite close to equilibrium values predicted by theoretical models. Mean dune steepnesses were also similar to predicted equilibrium values, except during high, falling flows when the steepness was above that predicted. Variations in mean dune height and length with discharge are similar to those predicted by theory under conditions of low mean dune excursion and discharge variation with a short high water period and long low water period. However, the calculated rates of change of height of individual dunes vary considerably from those measured. Rates of dune creation and destruction were unrelated to discharge variations, contrary to previous results. Instead, creations and destructions were apparently the result of local variations in bed shear stress and sediment transport rate. Observed changes in dune height during unsteady flows agree with theory fairly well at low bed shear stresses, but not at higher bed shear stresses when suspended sediment transport is significant.     

Validation and implications of an energy‐based bedload transport equation

A recently developed bedload equation (Abrahams & Gao, 2006) has the form i_b = ωG^3˙4, where i_b is the immersed bedload transport rate, ω is the stream power per unit area, G = 1?θ_c/θ, θ is the dimensionless shear stress and θ_c is the associated threshold value for the incipient motion of bed grains. This equation has a parsimonious form and provides good predictions of transport rate in both the saltation and sheetflow regimes (i.e. flows with low and high θ values, respectively). In this study, the equation was validated using data independent of those used for developing it. The data represent bedload of identical sizes transported in various steady, uniform, fully rough and turbulent flows over plane, mobile beds. The equation predicted i_b quite well over five orders of magnitude. This equation was further compared with six classic bedload equations and showed the best performance. Its theoretical significance was subsequently examined in two ways. First, based on collision theory, the parameter G was related to the ratio of grain‐to‐grain collisions to the total collisions including both grain‐to‐grain and grain‐to‐bed collisions, P_g by P_g = G², suggesting that G characterizes the dynamic processes of bedload transport from the perspective of granular flow, which partly accounts for the good performance of the equation. Moreover, examining the ability of two common equations to predict bedload in gravel‐bed rivers revealed that G can also be used to simplify equations for predicting transport capacities in such rivers. Second, a simple dimensionless form of the equation was created by introducing B = i_b/ω. The theoretical nature of the term B was subsequently revealed by comparing this equation with both the Bagnold model and two commonly used parameters representing dimensionless bedload transport rates.     

An experimental study of flow, bedload transport and bed topography under conditions of erosion and deposition and comparison with theoretical models

The nature of flow, sediment transport and bed texture and topography was studied in a laboratory flume using a mixed size-density sediment under equilibrium and non-equilibrium (aggradational, degradational) conditions and compared with theoretical models. During each experiment, water depth, bed and water surface elevation, flow velocity, bed shear stress, bedload transport and bed state were continuously monitored. Equilibrium, uniform flow was established with a discharge of about 0.05 m³ s^?1, a flow depth of about 0.01 m, a flow velocity of about 0.81–0.88 m s^?1, a spatially averaged bed shear stress of about 1.7–2.2 Pa and a sediment transport rate of about 0.005–0.013 kg m^?1 s^?1 (i.e. close to the threshold of sediment transport). Such equilibrium flow conditions were established prior to and at the end of each aggradation or degradation experiment. Pebble clusters, bedload sheets and low-lying bars were ubiquitous in the experiments. Heavy minerals were relatively immobile and occurred locally in high concentrations on the bed surface as lag deposits. Aggradation was induced by (1) increasing the downstream flow depth (flume tilting) and (2) sediment overloading. Tilt-induced aggradation resulted in rapid deposition in the downstream half of the flume of a cross-stratified deposit with downstream dipping pebbles (pseudo-imbricated). and caused a slight decrease in the equilibrium mean water surface slope and total bedload transport rate. These differences between pre- and post-aggradation equilibrium flow conditions are due to a decrease in the local grain roughness of the bed. Sediment overloading produced a downstream fining and thinning wedge of sediment with upstream dipping pebbles (imbricated), whereas the equilibrium flow and sediment transport conditions remained relatively unchanged. Degradation was induced by (1) decreasing the downstream flow depth (flume tilting) and (2) cutting off the sediment feed. Tilt-induced degradation produced rapid downstream erosion and upstream deposition due to flow convergence with little change to the equilibrium flow and sediment transport conditions. The cessation of sediment feed produced degradation and armour development, a reduction in the mean water surface slope and flow velocity, an increase in flow depth, and an exponential decrease in bedload transport rate as erosion proceeded. A bedload transport model predicted total and fractional transport rates extremely well when the coarse-grained (or bedform trough) areas of the bed are used to define the sediment available to be transported. A sediment routing model, MIDAS, also reproduced the equilibrium and non-equilibrium flow conditions, total and fractional bedload transport rates and changes in bed topography and texture very well.     

Grain-size distributions of bed load: Inferences from flume experiments using heterogeneous sediment beds

The grain-size fractions in the bedload transported over the five heterogeneous sediment beds of different values of bed roughness were computed from the flume experiments. The existence of an entrapment factor associated with the sorting observed from the bed to active layer was modeled based on the modified critical shear stress to estimate the grain-size fractions in the transport layer under given hydraulic conditions. The efficiency of these models was tested with the observed data. Subsequently, the patterns of observed grain-size distributions in the transport layer were tested to identify the distributions developed in the active layer due to sorting using three probability density functions (pdf), such as, log-normal, log-hyperbolic and log-skew-Laplace. Tests indicated that a log-skew-Laplace distribution fitted best for 49%, a log-hyperbolic for 31%, and a log-normal for 20% out of forty-five bedload samples collected under unidirectional flow with changes in flow discharge and bed roughness. The results of this study would be useful to specify the grain-size distributions in the bedload formed under different hydrodynamic conditions in various sedimentary environments.     

长江中游城汉河段推移质输沙率计算公式

三峡水库蓄水运用后,城汉河段悬移质输沙量大幅度减少,推移质输沙量所占比重增加,造床作用日益凸显,故开展其推移质输沙率计算公式研究具有理论与实际意义。利用螺山和汉口水文站实测资料建立了推移质输沙率与流量之间的指数关系式,并据此推算了两站逐日推移质输沙率,结果表明:螺山站、汉口站输沙率均与流量的0.912 78次方成正比,多年（2009—2015年）平均推移质输沙量分别为137万t、152万t,主要集中在汛期。利用推移质实测资料对Engelund、Einstein、Yalin 3个公式进行了检验与修正,结果表明:修正前各公式计算结果比较分散,且与实测值偏差很大;修正后各公式计算精度显著提高,综合比较发现,修正后的Yalin公式精度最高,Engelund公式、指数关系式次之,Einstein公式精度相对较低。因此,修正后的Yalin公式更适合于城汉河段推移质输沙率计算,可用于该河段的演变分析与数学模型计算。     

The relationship of sediment transport rates and paths to sandbanks in a tidally dominated area off the coast of East Anglia, U.K.

ABSTRACT Sand transport measurements of bedload and suspended load in the Sizewell-Dunwich Banks area, East Anglia have shown that the suspended mode is dominant. The depth-integrated spring tidal residual is 5.66 g cm⁻¹ sec⁻¹, although the neap rate is only one-fifth of this. The calculated bedload transport rates also vary, from 0.012 to 0.040 g cm⁻¹ sec⁻¹, correlating with changing meteorological conditions.
In order to predict the bedload sediment circulation pattern from midwater current meter measurements, five sediment transport equations were calibrated, using fluorescent dyed sand. Yalin's relationship gave the best estimates. The bed shear stress was determined by extrapolating the velocity profile as a power law relationship, with an exponent equal to 0.1, from midwater down to 2 m and as a lognormal profile from 2 m to the sea-bed. Roughness length values appropriate to the substrate were used.
Although bedload transport residuals are mainly to the south, the banks trend northwards from the coast and have also elongated in this direction. This is thought to be in response to the dominance of the suspended sediment transport. It is suggested that a tidal residual eddy mechanism is responsible for the banks'maintenance, similar to the process operating in Start Bay, Devon. The well-documented westward movement of the banks is likely to be related to wave processes.     

The response of subaqueous dunes to floods in sand and gravel bed reaches of the Dutch Rhine

Abstract The branches of the River Rhine in the Netherlands, characterized by a sand–gravel bed in the upstream part and a sand bed in the downstream part of the river system, show migrating dunes, especially during floods. In the last 20 years, these dunes have been studied extensively. High-resolution echo-sounding measurements of these dunes, made with single and multibeam equipment, were analysed for three different sections of the Rhine river system during several floods. This analysis was done to quantify the growth, decay and migration rates of the dunes during floods. In addition, the migrating dunes were used to calculate bedload transport rates with dune tracking. The results of dune growth and decay and migration rate are shown to be very different for the various sections during the various floods, and these differences are related to differences in grain size of the bed and to differences in the distribution of discharge over the main channel and the floodplain. The relations are used to show that the growth and migration rate of dunes, and the calculated bedload transport rates during the rising stage of a flood wave can be predicted from the mobility of the bed material with simple power relations.     

In situ acoustic measurements of marine gravel threshold and transport

Measurements of the nearbed turbulent current flow and the bedload transport of marine gravel have been made over three tidal cycles. The turbulence in the bottom boundary layer was measured using two electromagnetic current meters, and the gravel transport was measured using a passive acoustic system which monitored the interparticle collision noise of locally mobile material. Visual estimates of bedload were also obtained with an underwater TV camera. The acoustic technique, unlike a conventional bedload sampler, has allowed estimates of transport to be obtained with a temporal resolution comparable with the turbulence data collected. This has enabled a detailed comparison to be made between the turbulent flow and the sediment response to the instantaneous flow conditions. The results of the study show that of the turbulent bursting events which contribute towards the Reynolds stress, only the sweeps and outward interactions play a significant role in the transport of coarse sedimentary material. The measurements show that it is the instantaneous increases in the horizontal turbulent velocity fluctuations that generate excess shear stresses which drive the transport process.     

Sediment pathways and transport rates through a tide-dominated entrance, Rangaunu Harbour, New Zealand

Bed load sediment traps were deployed at two sections across channels in Rangaunu Harbour entrance. Traps were inspected and emptied by divers at hourly intervals through both spring and neap tidal cycles for a total of 292 trap deployments. Current velocities were measured simultaneously with the trap inspections. Transport is concentrated in sandy megaripple fields on the channel banks and sub-tidal platforms flanking the channels. There, transport is almost continuous throughout the tidal cycle, increasing with flow velocity but lagging by approximately one hour. The channel floors are lined with shell-gravel lag across which bedload transport rates are low and discontinuous. Tidal asymmetry produces a net seaward transport through the channel troughs and a net landward transport across the channel banks and flanking sub-tidal platforms. Sediment leaving the harbour recirculates in anticlockwise gyres across the ebb-tide delta to re-enter the harbou and maintain the supply of sand to the megaripple field. Transport during spring tides is typically 25–30 times that during neaps. Predictions of transport rates, from a method developed by Black & Healy utilizing the Yalin bedload equation, produced transport rates similar to the traps over sand beds. Transport over shell lag surfaces appears independent of near-bed velocity and more dependent on the passage of ribbons of sand across the lag surface.     

A method for the measurement of bedload sediment transport and passive faunal transport on intertidal sandflats

A simple and inexpensive sampler to measure bedload sediment transport in shallow subtidal or intertidal areas is described. The cylindrical sub-sediment trap with an aspect ratio of 20 (height: diameter) is an improvement over conventional bedload samplers which are difficult to use in shallow areas or fail to collect the biological material associated with bedload. Traps deployed on a low-energy intertidal sandflat for six months provided daily estimates of bedload transport (quartz grains: 0.001–40 kg m^?1 d^?1), passive infaunal transport (e.g., the bivalveMya arenaria, max: 800 ind m^?1 d^?1), and organic detrital flux (e.g., macrophyte fragments, max: 400 g dry wt m^?1 d^?1). Bedload rates estimated with traps were compared to predictions from a numerical bedload model to evaluate the trap’s collection and retention efficiency. A significant linear regression between observed (trap) and predicted (model) rates (r²=0.65, p<0.001, n=97) indicated that the traps were useful for the measurement of high- and low-frequency variability in bedload transport. Potential applications of the traps in benthic oceanography include recruitment and recolonization studies.     

Integration of Tracer Test Data to Refine Geostatistical Hydraulic Conductivity Fields Using Sequential Self-Calibration Method

On the basis of local measurements of hydraulic conductivity,geostatistical methods have been found to be useful in heterogeneity characterization of a hydraulic conductivity field on a regional scale. However,the methods are not suited to directly integrate dynamic production data,such as,hydraulic head and solute concentration,into the study of conductivity distribution. These data,which record the flow and transport processes in the medium,are closely related to the spatial distribution of hydraulic conductivity. In this study,a three-dimensional gradient-based inverse method-the sequential self-calibration (SSC) method-is developed to calibrate a hydraulic conductivity field,initially generated by a geostatistical simulation method,conditioned on tracer test results. The SSC method can honor both local hydraulic conductivity measurements and tracer test data. The mismatch between the simulated hydraulic conductivity field and the reference true one,measured by its mean square error (MSE),is reduced through the SSC conditional study. In comparison with the unconditional results,the SSC conditional study creates the mean breakthrough curve much closer to the reference true curve,and significantly reduces the prediction uncertainty of the solute transport in the observed locations. Further,the reduction of uncertainty is spatially dependent,which indicates that good locations,geological structure,and boundary conditions will affect the efficiency of the SSC study results.