Imaging of water flow in porous media by magnetic resonance imaging microscopy

Magnetic resonance imaging (MRI) was used to study the flow of water in a column 14 mm in diameter packed with glass heads. The sample was fully saturated and water was pumped through the column using a peristaltic pump, at flow rates of 125 and 250 mL h(-1). This corresponds to mean velocities of 0.5 and 1 mm s(-1), given a porosity of 0.46 m3 m(-3). Nuclear magnetic resonance (NMR) images of the proton density and velocities within a 2-mm slice were taken at a spatial resolution of 0.15 x 0.15 x 2 mm3. At a mean pore water velocity of 1 mm s(-1) we approximated hydrodynamic dispersion using NMR-measured velocity distributions in a 2-mm slice through the sample. Additionally, we conducted a step pulse tracer experiment with chloride through the same column and at identical initial and boundary conditions. We fitted the convection-dispersion equation to the breakthrough curve and compared the column scale dispersion of the tracer experiment with the respective NMR estimate derived at the slice scale.     

Three-dimensional nickel ion transport through porous media using magnetic resonance imaging

The transport of Ni2+ ions in a column, filled with porous media, was observed in three dimensions and time by magnetic resonance imaging (MRI) in a clinical scanner. For porous media we used glass beads or quartz sand in a saturated continuous flow mode. The magnetic moment of Ni2+ decreased the T1 relaxation time of 1H in aqueous solution. This concentration-dependent effect was used by a fast low angle shot (FLASH) MRI sequence for imaging the concentration of the dissolved ions. Since Ni2+ behaves as a conservative tracer under the chosen conditions, the tracer motion was representative for the water flow in the porous medium. Currently, we can achieve an isotropic spatial resolution of 1.5 mm and a temporal resolution of 170 s. The transport observation gives direct access to hydraulic flow properties of the porous media. The fluid flow velocity field was calculated by a fronttracking method and the statistical properties of the velocities were investigated. We also compared the experimental data with the three-dimensional particle tracking model PARTRACE, which uses the experimental flow field as input.     

Satellite Rainfall Estimates Over South America – Possible Applicability to the Water Management of Large Watersheds1

Pereira Filho, Augusto J., Richard E. Carbone, John E. Janowiak, Phillip Arkin, Robert Joyce, Ricardo Hallak, and Camila G.M. Ramos, 2010. Satellite Rainfall Estimates Over South America – Possible Applicability to the Water Management of Large Watersheds. Journal of the American Water Resources Association (JAWRA) 46(2):344-360. DOI: 10.1111/j.1752-1688.2009.00406.x Abstract: This work analyzes high-resolution precipitation data from satellite-derived rainfall estimates over South America, especially over the Amazon Basin. The goal is to examine whether satellite-derived precipitation estimates can be used in hydrology and in the management of larger watersheds of South America. High spatial-temporal resolution precipitation estimates obtained with the CMORPH method serve this purpose while providing an additional hydrometeorological perspective on the convective regime over South America and its predictability. CMORPH rainfall estimates at 8-km spatial resolution for 2003 and 2004 were compared with available rain gauge measurements at daily, monthly, and yearly accumulation time scales. The results show the correlation between satellite-derived and gauge-measured precipitation increases with accumulation period from daily to monthly, especially during the rainy season. Time-longitude diagrams of CMORPH hourly rainfall show the genesis, strength, longevity, and phase speed of convective systems. Hourly rainfall analyses indicate that convection over the Amazon region is often more organized than previously thought, thus inferring that basin scale predictions of rainfall for hydrological and water management purposes have the potential to become more skillful. Flow estimates based on CMORPH and the rain gauge network are compared to long-term observed average flow. The results suggest this satellite-based rainfall estimation technique has considerable utility. Other statistics for monthly accumulations also suggest CMORPH can be an important source of rainfall information at smaller spatial scales where in situ observations are lacking.     

Modelling deposition and air concentration of reduced nitrogen in Poland and sensitivity to variability in annual meteorology

The relative contribution of reduced nitrogen to acid and eutrophic deposition in Europe has increased recently as a result of European policies which have been successful in reducing SO(2) and NO(x) emissions but have had smaller impacts on ammonia (NH(3)) emissions. In this paper the Fine Resolution Atmospheric Multi-pollutant Exchange (FRAME) model was used to calculate the spatial patterns of annual average ammonia and ammonium (NH(4)(+)) air concentrations and reduced nitrogen (NH(x)) dry and wet deposition with a 5 km × 5 km grid for years 2002-2005. The modelled air concentrations of NH(3) and dry deposition of NH(x) show similar spatial patterns for all years considered. The largest year to year changes were found for wet deposition, which vary considerably with precipitation amount. The FRAME modelled air concentrations and wet deposition are in reasonable agreement with available measurements (Pearson's correlation coefficients above 0.6 for years 2002-2005), and with spatial patterns of concentrations and deposition of NH(x) reported with the EMEP results, but show larger spatial gradients. The error statistics show that the FRAME model results are in better agreement with measurements if compared with EMEP estimates. The differences in deposition budgets calculated with FRAME and EMEP do not exceed 17% for wet and 6% for dry deposition, with FRAME estimates higher than for EMEP wet deposition for modelled period and lower or equal for dry deposition. The FRAME estimates of wet deposition budget are lower than the measurement-based values reported by the Chief Inspectorate of Environmental Protection of Poland, with the differences by approximately 3%. Up to 93% of dry and 53% of wet deposition of NH(x) in Poland originates from national sources. Over the western part of Poland and mountainous areas in the south, transboundary transport can contribute over 80% of total (dry + wet) NH(x) deposition. The spatial pattern of the relative contribution of national sources to total deposition of NH(x) may change significantly due to the general circulation of air.     

WATER QUALITY DURING STORM EVENTS FROM TWO CONSTRUCTED WETLANDS RECEWING MINE DRAINAGE1

ABSTRACT: Flow rates, pH, iron concentration, and manganese concentration were measured during several storm events at two constructed wetlands receiving mine water. During a substantial rain event, flow rates at both the wetland outlets surpassed flow rates at the wetland inlets, reflecting incident rainfall and differences in wetland area at the two sites. A significant positive correlation existed between local rainfall and outflow rates at the larger wetland, but not between rainfall and inflow rates. During storm events, outlet pH, relative to inlet pH, was slightly elevated at the larger wetland, and depressed at the smaller wetland. However, over the course of one year, rainfall was uncorrelated to outlet pH in the larger wetland. A substantial rain event at the smaller wetland resulted in a temporary elevation in outlet iron concentrations, with treatment efficiency reduced to near zero. However, in the larger wetland, outlet iron concentrations were not significantly affected by storm events. Although rainfall and outlet iron concentration were not significant correlates at the larger wetland, flow rate was positively correlated to outlet iron concentration. A normal manganese treatment efficiency of 50 percent at the smaller wetland was reduced to zero during a heavy rain.     

Comparison of Measured and MM5 Modeled Meteorology Data for Simulating Flow in a Mountain Watershed1

Herr, Joel W., Krish Vijayaraghavan, and Eladio Knipping, 2010. Comparison of Measured and MM5 Modeled Meteorology Data for Simulating Flow in a Mountain Watershed. Journal of the American Water Resources Association (JAWRA) 46(6):1255–1263. DOI: 10.1111/j.1752-1688.2010.00489.x Abstract: Accurate simulation of time-varying flow in a river system depends on the quality of meteorology inputs. The density of meteorology measurement stations can be insufficient to capture spatial heterogeneity of precipitation, especially in mountainous areas. The Watershed Analysis Risk Management Framework (WARMF) model was applied to the Catawba River watershed of North and South Carolina to simulate flow and water quality in rivers and a series of 11 reservoirs. WARMF was linked with the AMSTERDAM air model to analyze the water quality benefit from reduced atmospheric emissions. The linkage requires accurate simulation of meteorology for all seasons and for all types of precipitation events. WARMF was driven by the mesoscale meteorology model MM5 processed by the Meteorology Chemistry Interface Processor, which provides greater spatial density but less accuracy than meteorology stations. WARMF was also run with measured data from the National Climatic Data Center (NCDC) to compare the performance of the watershed model using measured data vs. modeled meteorology as input. A one year simulation using MM5 modeled meteorology performed better overall than the simulation using NCDC data for the volumetric water balance measure used for calibration, but MM5 represented precipitation from a dissipated hurricane poorly, which propagated into errors of simulated flow.     

Evaluation of a Watershed Model for Estimating Daily Flow Using Limited Flow Measurements1

Abstract: The Soil and Water Assessment Tool (SWAT) model was evaluated for estimation of continuous daily flow based on limited flow measurements in the Upper Oyster Creek (UOC) watershed. SWAT was calibrated against limited measured flow data and then validated. The Nash‐Sutcliffe model Efficiency (NSE) and mean relative error values of daily flow estimations were 0.66 and 15% for calibration, and 0.56 and 4% for validation, respectively. Also, further evaluation of the model’s estimation of flow at multiple locations was conducted with parametric paired t‐test and nonparametric sign test at a 95% confidence level. Among the five main stem stations, four stations were statistically shown to have good agreement between predicted and measured flows. SWAT underestimated the flow of the fifth main stem station possibly because of the existence of complex flood control measures near to the station. SWAT estimated the daily flow at one tributary station well, but with relatively large errors for the other two tributaries. The spatial pattern of predicted flows matched the measured ones well. Overall, it was concluded from the graphical comparisons and statistical analyses of the model results that SWAT was capable of reproducing continuous daily flows based on limited flow data as is the case in the UOC watershed.     

MODELING STREAMFLOW FROM ARTIFICIALLY DRAINED AGRICULTURAL WATERSHEDS IN ILLINOIS1

ABSTRACT: Evaluation of the applicability and validity of hydrologic simulation models for various cropping systems in different hydrogeologic and soil conditions is needed for a range of spatial scales. We calibrated and tested the ADAPT model for simulating streamflow from 552 to 1,985 km² watersheds in central Illinois, where more than 79 percent of the land is used for maize‐soybean production and tile drainage is common. Model calibration was performed with a seven year period (1987 through1993) of measured streamflow from one of the watersheds, and model testing was done using independent weather and measured streamflow data from the two neighboring watersheds for the same seven year period. Simulations of annual streamflow were accurate with a coefficient of determination and Willmott's index of agreement of 0.98 and 0.99, respectively. For simulation of monthly streamflow, Willmott's index of agreement ranged from 0.93 to 0.95. For simulation of daily streamflow, Willmott's index of agreement ranged from 0.84 to 0.85. The daily simulations challenged the temporal and spatial resolution of our measured precipitation data. Discrepancies between simulated and measured data may result from the model's inability to effectively address frozen soils and snowmelt runoff processes and in accurately representing evapotranspiration.     

Fluorescence imaging applied to tracer distributions in variably saturated fractured clayey till

The study of mechanisms controlling preferential flow and transport in variably saturated fractured clayey till is often hindered by insufficient spatial resolution or unknown measuring volume. With the objective to study these mechanisms while circumventing the obstacles, tracer experiments with two fluorescent tracers Acid Yellow 7 (AY7) and Sulforhodamine B (SB) were performed at three different rain events for a fall and a summer season. Irrigated areas were excavated down to depths of 2.8 m and the movement of both tracers in the exposed profiles was delineated simultaneously by high spatial resolution apparent concentration maps (pixel approximately 1 mm(2)) obtained with an imaging device. The device consists of a light source and a CCD camera, both equipped with tracer-specific-filters for fluorescent light. The fluorescence images were corrected for nonuniform lighting, changing surface roughness, and varying optical properties of the soil profile. The resulting two-dimensional apparent concentration distribution profiles of the tracers showed that: (i) relative low water content in the upper 10 cm of the irrigated till in summer had a pronounced retardation effect on the AY7-migration and no effect on the SB-migration; (ii) the dead-end biopores were not activated in the fall season; (iii) only 3D fracture-plans connected to hydraulically active 1D-biopores contributed to the leaching; (iv) the tracer migration primary followed macropores during both seasons, though AY7 also followed a topsoil piston transport in summer; (v) the highest tracer pixel apparent concentrations were often found in macropores and most pronounced in the summer season; and (vi) 3D-dilution in fractures seems to play a dominating role in AY7-migration in the fall season.     

Using heat to characterize streambed water flux variability in four stream reaches

Estimates of streambed water flux are needed for the interpretation of streambed chemistry and reactions. Continuous temperature and head monitoring in stream reaches within four agricultural watersheds (Leary Weber Ditch, IN; Maple Creek, NE; DR2 Drain, WA; and Merced River, CA) allowed heat to be used as a tracer to study the temporal and spatial variability of fluxes through the streambed. Synoptic methods (seepage meter and differential discharge measurements) were compared with estimates obtained by using heat as a tracer. Water flux was estimated by modeling one-dimensional vertical flow of water and heat using the model VS2DH. Flux was influenced by physical heterogeneity of the stream channel and temporal variability in stream and ground-water levels. During most of the study period (April-December 2004), flux was upward through the streambeds. At the IN, NE, and CA sites, high-stage events resulted in rapid reversal of flow direction inducing short-term surface-water flow into the streambed. During late summer at the IN site, regional ground-water levels dropped, leading to surface-water loss to ground water that resulted in drying of the ditch. Synoptic measurements of flux generally supported the model flux estimates. Water flow through the streambed was roughly an order of magnitude larger in the humid basins (IN and NE) than in the arid basins (WA and CA). Downward flux, in response to sudden high streamflows, and seasonal variability in flux was most pronounced in the humid basins and in high conductivity zones in the streambed.     

长江三峡库区横向扩散系数研究初探

进行长江长寿段示踪剂岸边有限时段连续投放试验,求得横向扩散系数及其与流速、水深、流量、污染物浓度等关系。以长寿湖大沙坝河段为类比试验,预测得到三峡库区在90％保证率的最枯月流量下典型江段的横向扩散系数,并可用于库区的污染带计算。     

Phosphorus transport in overland flow in response to position of manure application

Phosphorus (P) loss in overland flow varies with spatial distribution of soil P, management, and hydrological pathways. The effect of flow time, flowpath length, and manure position on P loss in overland flow from two central Pennsylvania soils packed in boxes of varying length (0.5, 1.0, 1.5, 2.75, and 4.0 m long x 15 cm wide x 5 cm deep) were examined by collecting flow samples at 5-min intervals for 30 min (50 mm h(-1) rainfall) without and with 75 kg P ha(-1) applied as swine (Sus scrofa) manure over 0.5 m of the box slope length at distances of 0 to 3.5 m from the downslope collection point. Dissolved reactive P concentration was more closely related to the proportion of clay in sediment of overland flow before (r = 0.98) than after (r = 0.56) manure application. This was attributed to the transport of larger, low-density particles after applying manure. The concentration of dissolved and particulate P fractions decreased with increasing flowpath length, due to dilution rather than sorption of P by surface soil during overland flow. Total P loss (mainly as particulate P) from the Watson channery silt loam (fine-loamy, mixed, active, mesic Typic Fragiudult) was more than from Berks channery silt loam (loamy-skeletal, mixed, active, mesic Typic Dystrudept), even with manure applied. Thus, while P loss in overland flow is affected by where manure is applied relative to flowpath length, initial soil P concentration should not be discounted when looking at areas of potential P loss within a watershed.     

Determining long-term (decadal) deep drainage rate using multiple tracers

The deep drainage rate is a critical hydrological parameter in understanding contamination mechanisms of soil and groundwater. Little research has been conducted on the temporal variations in deep drainage rate during the last century. The objective of this study was to determine the long-term deep drainage rate on a cultivated loamy soil in the Canadian Prairies. Three tracers were used: KCl applied in 1971, fallout tritium in 1963, and NO3* released during the initial cultivation of the field (1923). Two soil cores to a depth of 3.6 m were taken along a flat portion of the field, and soil Cl(-), 3H, and NO3* concentrations were measured as a function of depth. An additional four cores were taken for soil water content measurements between 2000 and 2003. Distinct peaks in the depth distribution of these three tracers were located at 1.27 m for Cl(-), 1.31 m for 3H, and 1.52 m for NO3*, 32, 40, and 80 yr after the application of Cl(-), 3H, and NO3*, respectively. The average deep drainage rates, calculated as the product of the estimated tracer velocity and volumetric soil water content below the active root zone, were 2.0 mm yr(-1) from the Cl(-) tracer, 2.2 mm yr(-1) from 3H, and 2.5 mm yr(-1) from the NO3* tracer. Therefore, there was little temporal variability in the groundwater recharge over the eight decades that the field has been cultivated. The recharge rates are less than 1% of the mean annual precipitation (333 mm).     

Influence of prairie restoration on CT-measured soil pore characteristics

Restored prairies are expected to improve soil physical properties, yet little is known about the extent of change to soil properties and how rapidly these changes take place. The objective of this study was to compare effects of prairie restoration on computed tomography (CT)-measured pore parameters. Undisturbed soil cores (76 mm diam. by 76 mm long) from native prairie (NP), restored prairie (RP), conservation reserve program (CRP), and no-till corn (Zea mays L.)-soybean (Glycine max (L.) Merr.; CS) sites were collected with six replicates from the 0- to 40-cm depth in 10-cm increments. Five CT images were acquired from each soil core using a medical CT scanner with 0.2 by 0.2 mm pixel resolution with 0.5 mm slice thickness, and then images were analyzed. Soil bulk density and hydraulic conductivity (K(sat)) were also measured. Soils under NP, RP, CRP, and CS areas had 83, 43, 48, and 26 pores on a 2500 mm(2) area, respectively, for the 0- to 40-cm depth. The number of pores, number of macropores (>1000 microm diam.), macroporosity, mesoporosity (200-1000 microm diam.), and fractal dimension were significantly higher and pore circularity was lower for NP, RP, and CRP than the CS treatment. The CT-measured mesoporosity and macroporosity of the CS treatment were 20 and 18% of the values for the NP site. CT-measured number of pores and macropores explained 43 and 40% of the variation for K(sat). The study showed that prairie restoration improves CT-measured soil pore parameters and decreases bulk density which are related to soil water infiltration.     

Modeling macropore flow effects on pesticide leaching: inverse parameter estimation using microlysimeters

Macropore flow is a key factor determining pesticide fate, but models accounting for this process need parameters that cannot be easily measured. This study was conducted to investigate the use of inverse techniques to estimate parameters controlling macropore flow and pesticide fate in the dual-permeability model MACRO. Undisturbed columns were sampled at three landscape positions (hilltop, slope, hollow) with contrasting texture and organic carbon content. Transient leaching experiments were performed for an anionic tracer and the herbicide MCPA (4-chloro-2methylphenoxy acetic acid) during a 4-mo period, first under natural rainfall, and then under controlled irrigation in the laboratory. The tracer breakthrough for the liner-textured soil from the hilltop showed strong evidence of macropore flow, resulting in a rapid leaching of MCPA, while leaching was minimal from the organic-rich hollow soil, since macropore flow was weaker and adsorption stronger. The MACRO model was linked to the inverse modeling program SUFI (Sequential Uncertainty Fitting) to enable calibration of nine key model parameters. Based on calculated model efficiencies, MACRO-SUFI gave generally good predictions of water movement and tracer and pesticide transport, although some errors were attributed to difficulties in simulating the effects of soil moisture on degradation and the timing of water outflows. Even after calibration, significant uncertainties remained for some key parameters controlling macropore flow. Nevertheless, the parameter estimates were significantly different between landscape positions and could also be related to basic soil properties. The posterior uncertainty ranges could probably be reduced with a more exhaustive sampling of the parameter space and improved experimental designs.     

Effects of soil data resolution on SWAT model stream flow and water quality predictions

The prediction accuracy of agricultural nonpoint source pollution models such as Soil and Water Assessment Tool (SWAT) depends on how well model input spatial parameters describe the characteristics of the watershed. The objective of this study was to assess the effects of different soil data resolutions on stream flow, sediment and nutrient predictions when used as input for SWAT. SWAT model predictions were compared for the two US Department of Agriculture soil databases with different resolution, namely the State Soil Geographic database (STATSGO) and the Soil Survey Geographic database (SSURGO). Same number of sub-basins was used in the watershed delineation. However, the number of HRUs generated when STATSGO and SSURGO soil data were used is 261 and 1301, respectively. SSURGO, with the highest spatial resolution, has 51 unique soil types in the watershed distributed in 1301 HRUs, while STATSGO has only three distributed in 261 HRUS. As a result of low resolution STATSGO assigns a single classification to areas that may have different soil types if SSURGO were used. SSURGO included Hydrologic Response Units (HRUs) with soil types that were generalized to one soil group in STATSGO. The difference in the number and size of HRUs also has an effect on sediment yield parameters (slope and slope length). Thus, as a result of the discrepancies in soil type and size of HRUs stream flow predicted was higher when SSURGO was used compared to STATSGO. SSURGO predicted less stream loading than STATSGO in terms of sediment and sediment-attached nutrients components, and vice versa for dissolved nutrients. When compared to mean daily measured flow, STATSGO performed better relative to SSURGO before calibration. SSURGO provided better results after calibration as evaluated by R(2) value (0.74 compared to 0.61 for STATSGO) and the Nash-Sutcliffe coefficient of Efficiency (NSE) values (0.70 and 0.61 for SSURGO and STATSGO, respectively) although both are in the same satisfactory range. Modelers need to weigh the benefits before selecting the type of data resolution they are going to use depending on the watershed size and level of accuracy required because more effort is required to prepare and calibrate the model when a fine resolution soil data is used.     

Research on distribution performance of desert sand for heat storage in downcomer

In this research, desert sand is used as the sensible heat storage medium, which exchanges heat with air in the downcomer to realize heat storage and heat release. The desert sand distribution uniformity has a significant impact on the heat exchange performance and efficiency between desert sand and air for the process of convection in the downcomer. Given the superiority of sensible heat storage in convective heat transfer between desert sand and air, distributors with cylinder or conical bore solid particles and homogeneity performance testing device are designed and manufactured on the basis of convection system equipped with solid particle–air downcomer. Then, the convection experiment between solid sand and air is researched. The greater the desert sand flow rate and higher the volume density, the larger the variance of regional mass flow rate and the worse the homogeneity performance. For the cylinder bore distributor, the smaller the sand particle size is, the greater affected the sand groups can be. The sand homogeneity performance is preferable with the two particle size ranges: 0.18-0.25 mm and 0.15-0.18 mm. The total sand flow rate decreases, but the uniformity improves with the increase of the air flow velocity, and the best distribution performance is achieved at an air velocity of 0.6 m/s. However, the distribution performance declines with the air flow velocity persistently increasing because the sand groups are pushed to one pipe side close to the wall. The sand groups deflect seriously with the air flow velocity increasing.     

Enhancing Hyrdological Simulation Program – FORTRAN Model Channel Hydraulic Representation1

Abstract: The hydrological simulation program – FORTRAN (HSPF) is a comprehensive watershed model that employs depth‐area‐volume‐flow relationships known as the hydraulic function table (FTABLE) to represent the hydraulic characteristics of stream channel cross‐sections and reservoirs. An accurate FTABLE determination for a stream cross‐section site requires an accurate determination of mean flow depth, mean flow width, roughness coefficient, longitudinal bed slope, and length of stream reach. A method that uses regional regression equations to estimate mean flow depth, mean flow width, and roughness coefficient is presented herein. FTABLES generated by the proposed method (Alternative Method) and FTABLES generated by Better Assessment Science Integrating Point and Nonpoint Sources (BASINS) were compared. As a result, the Alternative Method was judged to be an enhancement over the BASINS method. First, the Alternative Method employs a spatially variable roughness coefficient, whereas BASINS employs an arbitrarily selected spatially uniform roughness coefficient. Second, the Alternative Method uses mean flow width and mean flow depth estimated from regional regression equations whereas BASINS uses mean flow width and depth extracted from the National Hydrography Dataset (NHD). Third, the Alternative Method offers an option to use separate roughness coefficients for the in‐channel and floodplain sections of compound channels. Fourth, the Alternative Method has higher resolution in the sense that area, volume, and flow data are calculated at smaller depth intervals than the BASINS method. To test whether the Alternative Method enhances channel hydraulic representation over the BASINS method, comparisons of observed and simulated streamflow, flow velocity, and suspended sediment were made for four test watersheds. These comparisons revealed that the method used to estimate the FTABLE has little influence on hydrologic calibration, but greatly influences hydraulic and suspended sediment calibration. The hydrologic calibration results showed that observed versus simulated daily streamflow comparisons had Nash‐Sutcliffe efficiencies ranging from 0.50 to 0.61 and monthly comparisons had efficiencies ranging from 0.61 to 0.84. Comparisons of observed and simulated suspended sediments concentrations had model efficiencies ranging from 0.48 to 0.56 for the daily, and 0.28 to 0.70 for the monthly comparisons. The overall results of the hydrological, hydraulic, and suspended sediment concentration comparisons show that the Alternative Method yielded a relatively more accurate FTABLE than the BASINS method. This study concludes that hydraulic calibration enhances suspended sediment simulation performance, but even greater improvement in suspended sediment calibration can be achieved when hydrological simulation performance is improved. Any improvements in hydrological simulation performance are subject to improvements in the temporal and spatial representation of the precipitation data.     

Interactions between soil texture and placement of dairy slurry application: I. Flow characteristics and leaching of nonreactive components

Land application of manure can exacerbate nutrient and contaminant transfers to the aquatic environment. This study examined the effect of injecting a dairy cattle (Bostaurus L.) manure slurry on mobilization and leaching of dissolved, nonreactive slurry components across a range of agricultural soils. We compared leaching of slurry-applied bromide through intact soil columns (20 cm diam., 20 cm high) of differing textures following surface application or injection of slurry. The volumetric fraction of soil pores >30 microm ranged from 43% in a loamy sand to 28% in a sandy loam and 15% in a loam-textured soil. Smaller active flow volumes and higher proportions of preferential flow were observed with increasing soil clay content. Injection of slurry in the loam soil significantly enhanced diffusion of applied bromide into the large fraction of small pores compared with surface application. The resulting physical protection against leaching of bromide was reflected by 60.2% of the bromide tracer was recovered in the effluent after injection, compared with 80.6% recovery after surface application. No effect of slurry injection was observed in the loamy sand and sandy loam soils. Our findings point to soil texture as an important factor influencing leaching of dissolved, nonreactive slurry components in soils amended with manure slurry.     

A HYDRODYNAMIC MODELING STUDY TO DETERMINE THE OPTIMUM WATER INTAKE LOCATION IN LAKE PALDANG,KOREA1

A three‐dimensional hydrodynamic model was applied to Lake Paldang, South Korea. The lake has three inflows, of which Kyoungan Stream has the smallest flow rate and poorest water quality. Since all drinking water intake stations are located near the confluence of Kyoungan Stream within the lake, this contaminated tributary may have a significant impact on the quality of drinking water sources. The optimum drinking water intake location was determined from the applied model. The model was calibrated and verified using the data measured under different hydrological conditions. The model results were in reasonable agreement with the field measurements in both calibration and verification. The circulation and spreading patterns of the incoming flows in the lake, as well as their composition ratios to the drinking water intakes were determined from the model, and three alternative intake locations were proposed. The simulation results suggested that the horizontal and vertical relocations of the intake aqueduct could significantly decrease the composition ratio of the contaminated water. From this study, it was concluded that the three‐dimensional hydrodynamic model could successfully simulate the temporal and spatial mixing patterns of incoming flows and become a useful tool in determining the optimum water intake location in Lake Paldang.