About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Adsorption of fluoride on goethite surfaces-implications on dental epidemiology

Fluoride ion interaction with synthetically prepared goethite has been investigated over a range of pH values (4–9) and F^– concentrations (10^–3–10^–5 M). The amount of F^– retained by goethite suspensions was found to be a function of pH, media ionic strength, F^– concentration, and goethite concentration. The lowest ionic strength (0.001 M KNO₃) gave the highest adsorption medium. Uptake was minimal at pH >7 and increased with decreasing pH. Thermodynamic properties for fluoride adsorption at 298 K and 323 K were investigated. The isosteric heat of adsorption (H _r) was calculated and the heterogeneity and homogeneity of the surface examined for goethite. In view of the importance of fluoride in dental health, the interaction of fluoride on goethite in the physical environment has important implications on dental epidemiology.     

Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

An examination of directional discontinuities and magnetic polarity changes around interplanetary sector boundaries usingE > 2 keV electrons

Past studies of interplanetary magnetic sector boundaries have been based on the assumption that one can determine the field polarities by comparing the field directions with those of the nominal Parker spiral angles. Previous investigators have found evidence for decreases of B, the magnitude of the magnetic fieldB, and increases of , the angle betweenB and the ecliptic plane, at sector boundaries. Others have argued that the characteristic thickness of sector boundaries exceeds that of tangential discontinuities, making sector boundaries a separate class of structures.We use a simple technique for inferring the polarities of interplanetary magnetic fields based on the assumption thatE > 2 keV electrons are always flowing along the magnetic field away from the Sun. Electron data from the UC Berkeley experiment on the ISEE-3 spacecraft are used to examine periods around several apparent sector boundaries in 1978 and 1979. We compare properties of (a) boundaries with field polarity changes and (b) large-angle ( > 60°) directional discontinuities with no field polarity changes. We find no significant differences between the sector boundaries and the directional discontinuities in terms of associated decreases in B or of values of . These results suggest no significant difference between sector boundaries and directional discontinuities other than the change in field polarities. Within limited statistics we find that about half the polarity changes would not have been identified using a requirement that > 90° and that half of the > 120° discontinuities would have been misidentified as polarity changes.     

Doppler velocity measurements made with a scanning photoelectric magnetograph

We discuss the problems connected with the measurements and evaluation of line-of-sight velocities, obtained with a scanning photoelectric magnetograph using a line-shifter with enhanced sensitivity. We bring arguments for the validity of the results of our photoelectric Doppler velocity recordings. We have found a network of cellularly shaped patterns in the distribution of photo-electrically measured line-of-sight motions, upflowing in the magnetically quiet (blue-shifted) and downflowing in magnetically active (red-shifted) areas of the photosphere, if the mean velocity level is estimated for a sufficiently large measured area. The features of both directions are mutually complementary. We demonstrate the effect of the shift of the reference zero velocity level on the topology of the line-of-sight velocity maps, and the dependence of this level on the size of the area from which it is estimated.     

Dynamical evolution of NEAs: Close encounters, secular perturbations and resonances

We discuss the main mechanisms affecting the dynamical evolution of Near-Earth Asteroids (NEAs) by analyzing the results of three numerical integrations over 1 Myr of the NEA (4179) Toutatis. In the first integration the only perturbing planet is the Earth. So the evolution is dominated by close encounters and looks like a random walk in semimajor axis and a correlated random walk in eccentricity, keeping almost constant the perihelion distance and the Tisserand invariant. In the second integration Jupiter and Saturn are present instead of the Earth, and the 3/1 (mean motion) and v ₆ (secular) resonances substantially change the eccentricity but not the semimajor axis. The third, most realistic, integration including all the three planets together shows a complex interplay of effects, with close encounters switching the orbit between different resonant states and no approximate conservation of the Tisserand invariant. This shows that simplified 3-body or 4-body models cannot be used to predict the typical evolution patterns and time scales of NEAs, and in particular that resonances provide some fast-track dynamical routes from low-eccentricity to very eccentric, planet-crossing orbits.On leave from the Department of Mathematics, University of Pisa, Via Buonarroti 2, 56127 Pisa, Italy, thanks to the G. Colombo fellowships of the European Space Agency.     

Comparison of structural and least-squares lines for estimating geologic relations

Two different goals in fitting straight lines to data are to estimate a true linear relation (physical law) and to predict values of the dependent variable with the smallest possible error. Regarding the first goal, a Monte Carlo study indicated that the structural-analysis (SA) method of fitting straight lines to data is superior to the ordinary least-squares (OLS) method for estimating true straight-line relations. Number of data points, slope and intercept of the true relation, and variances of the errors associated with the independent (X) and dependent (Y) variables influence the degree of agreement. For example, differences between the two line-fitting methods decrease as error in X becomes small relative to error in Y. Regarding the second goal—predicting the dependent variable—OLS is better than SA. Again, the difference diminishes as X takes on less error relative to Y. With respect to estimation of slope and intercept and prediction of Y, agreement between Monte Carlo results and large-sample theory was very good for sample sizes of 100, and fair to good for sample sizes of 20. The procedures and error measures are illustrated with two geologic examples.