Estimation of groundwater evaporation and salt flux from Owens Lake, California, USA

Groundwater evaporation and subsequent precipitation of soluble salts at Owens Lake in eastern California have created one of the single largest sources of airborne dust in the USA, yet the evaporation and salt flux have not been fully quantified. In this study, we compare eddy correlation, microlysimeters and solute profiling methods to determine their validity and sensitivity in playa environments. These techniques are often used to estimate evaporative losses, yet have not been critically compared at one field site to judge their relative effectiveness and accuracy. Results suggest that eddy correlation methods are the most widely applicable for the variety of conditions found on large playa lakes. Chloride profiling is shown to be highly sensitive to thermal and density-driven fluxes in the near surface and, as a result, appears to underestimate yearly groundwater evaporation. Yearly mean groundwater evaporation from the playa surface estimated from the three study areas was found to range from 88 to 104 mm year⁻¹, whereas mean evaporation from the brine-covered areas was 872 mm year⁻¹. Uncertainties on these mean rates were estimated to be ±25%, based on comparisons between eddy correlation and lysimeter estimates. On a yearly basis, evaporation accounts for approximately 47 × 10⁶ m³ of water loss from the playa surface and open-water areas of the lake. Over the playa area, as much as 7.5 × 10⁸ kg (7.5 × 10⁵ t) of salt are annually concentrated by evaporation at or near the playa surface, much of which appears to be lost during dust storms in area.     

Optimal weighting in the finite difference solution of the convection-dispersion equation

Oscillation and numerical dispersion limit the reliability of numerical solutions of the convection-dispersion equation when finite difference methods are used. To eliminate oscillation and reduce the numerical dispersion, an optimal upstream weighting with finite differences is proposed. The optimal values of upstream weighting coefficients numerically obtained are a function of the mesh Peclet number used. The accuracy of the proposed numerical method is tested against two classical problems for which analytical solutions exist. The comparison of the numerical results obtained with different numerical schemes and those obtained by the analytical solutions demonstrates the possibility of a real gain in precision using the proposed optimal weighting method. This gain in precision is verified by interpreting a tracer experiment performed in a laboratory column.     

Extreme value theory based on the r largest annual events: a robust approach

In Smith (1986, J. Hydrol. 86, 27–43), a family of statistical distributions and estimators for extreme values based on a fixed number r > = 1 of the largest annual events are presented. The method of estimation was numerical maximum likelihood. In this paper, we consider the robust estimation of parameters in such families of distributions. The estimation technique, which is based on optimal B-robust estimates, will assign weights to each observation and give estimates of the parameters based on the data which are well modeled by the distribution. Thus, observations which are not consistent with the proposed distribution can be identified and the validity of the model can be assessed. The method is illustrated on Venice sea level data.     

The identifiability of parameters in a water quality model of the Biebrza River, Poland

The identifiability of model parameters of a steady state water quality model of the Biebrza River and the resulting variation in model results was examined by applying the Monte Carlo method which combines calibration, identifiability analysis, uncertainty analysis, and sensitivity analysis. The water quality model simulates the steady state concentration profiles of chloride, phosphate, ammonium, and nitrate as a function of distance along a river. The water quality model with the best combination of parameter values simulates the observed concentrations very well. However, the range of possible modelled concentrations obtained for other more or less equally eligible combinations of parameter values is rather wide. This range in model outcomes reflects possible errors in the model parameters. Discrepancies between the range in model outcomes and the validation data set are only caused by errors in model structure, or (measurement) errors in boundary conditions or input variables. In this sense the validation procedure is a test of model capability, where the effects of calibration errors are filtered out. It is concluded that, despite some slight deviations between model outcome and observations, the model is successful in simulating the spatial pattern of nutrient concentrations in the Biebrza River.     

Observations of particle movement in a monitoring well using the colloidal borescope

The colloidal borescope consists of a set of lenses and miniature video cameras capable of observing natural particles in monitoring wells. Based on field observations of these particles, it appears possible to measure in situ groundwater velocity in a well bore. Field observations have shown that directional measurements using the colloidal borescope are generally in good agreement with expected flow directions. However, the magnitude of flow velocity is higher compared with values based on conventional test methods. High relative flow velocities, even after correction factors have been applied to compensate for well bore effects, are believed to be due to preferential flow zones in the surrounding aquifer. Low flow zones exhibit swirling multidirectional flow that does not allow for a linear velocity measurement. Consequently, groundwater flow velocities measured by the colloidal borescope in heterogeneous aquifers will be biased toward the maximum velocity values present in the aquifer. A series of laboratory experiments was conducted to assess the reliability of the instrument. Based on this work, a seepage velocity correction factor ( ) of 1–4 was found for quantifying groundwater seepage velocity in the adjacent aquifer from observations in a well bore. Laboratory measurements also indicate that preferential flow in the surrounding aquifer dominates flow in the well. Results of this work suggest the possibility of quantifying higher-flow velocities associated with preferential flow zones in the subsurface.     

Methodical approach to distribution of epilithic and drifting algae communities in a French subalpine river: Inferences on water quality assessment

In a sub-alpine river, the Asse, with an unpredictable discharge regime, chlorophyll, density and taxonomic diversity of both drifting and periphytic communities were studied from a methodical approach. The investigations based on stamping two close cross-sections to know their spatial distributions in a heterogeneous bed substrate showed a heterogeneous pattern of colonization mainly influenced by hydrological variables as the current velocity, the size of the substrata (gravel-pebble-boulder). For the drifting algae, mainly composed with benthic species temporarily unhooked to the substrata, the heterogeneity of the distribution was reduced in fast flow which induced a mixing of the whole mass of water. For the periphyton, some results indicated that the greatest diversity and the highest chlorophyll concentration were found on pebble substratum where the current velocity was moderate. These results have also shown that the classical methods used to day for investigations were inadequate for most Mediterranean rivers. It was really difficult to estimate the water quality from classical index created for use in regular large plain rivers.     

Predisposition to Induced Seismicity in Some Czech Coal Mines

—Mining-induced seismicity associated with longwall face operations in the Ostrava- Karviná coal mines, Czech Republic, has been investigated in order to establish the conditions leading to a focal zone generation. The study, based on macroseismic and instrumental observations, proved that seismicity is influenced by natural as well as mining conditions. The first group includes the influence of faults, washouts and red beds, while the second one is represented by shaft and/or crosscut safety pillars and various types of remnant pillars. All the cases discussed show that many focal zones are generated in overstressed strata as a consequence of interaction of natural conditions and/or old workings with the active coal face. Received June 18, 1996, accepted March 10, 1997     

Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift

Recent technological advances in current measuring devices has resulted in a large observational database related to wind-driven motions in the upper ocean mixed layer. This has served to highlight the fact that transient motions make up a substantial contribution of the resulting Ekman currents. At the same time, certain discrepancies have emerged between the observed angular deflections of the steady-state currents from the surface wind stress, both at the surface and at sub-surface depths, which cannot be reconciled using the classical Ekman model. This paper seeks to tackle these two issues.First a general analytical method is presented for solving the time dependent horizontal momentum Ekman equations. Analysis of the unsteady terms that arise from simple special cases shows how the evolution proceeds through three stages. At early times, the Coriolis acceleration is insignificant, and the current is unidirectional and deepens through downward diffusion of momentum. Later Coriolis acceleration deflects the current vectors in the upper layers, whilst downward diffusion of momentum continues to deepen the layer. Finally, once diffusion has penetrated down to the depth of the steady-state current, then the transients decay on the inertial or diffusive timescale, depending upon the boundary conditions of the particular problem.In the second half of the paper, a new steady-state model is developed that includes the effects of wind-generated waves, through the action of their Stokes drift on the planetary vorticity. Comparisons between observations and the theoretical predictions, demonstrate that inclusion of the Stokes drift is the key to reconciling the discrepancies in the angular deflections of the steady-state currents. This leads to the conclusion that Ekman layer currents are significantly influenced by the surface waves.     

华南石炭纪岩关-大塘界线期牙形石地层分带

系统采集不同沉积相区和连续沉积的剖面，即浅海台地相广西柳江龙殿山剖面、台缘缓坡相云南施甸鱼硐剖面和台内盆地相广西柳州碰冲剖面，运用数理统计方法，研究了华南岩关一大塘界线过渡期牙形石的产出丰度和相分布，在恢复生物演化序列的基础上，对各剖面首先详细分带，再进一步综合成系统的5个牙形石带，即Gnathodus typicus—G．cuneiformis带、Gnathodus typicus—Protognathodus cordiformis带、Scaliognathus anchoralis—Gnathodus pseu-dosemiglaber带、Gnathodus praebilineatus带和Paragnathodus horrtopunctatus带。以生物系统发生为标准，严格限定各带的界线，并与北美、西欧同期地层进行精确对比，为建立杜内阶一维宪阶界线的全球层型(GSSP)提供可靠依据。     

Crystal morphologies in pillow basalts: implications for mid-ocean ridge processes

We present the results of a detailed petrological study of a sparsely phyric basalt (MAPCO CH98-DR11) dredged along the Mid-Atlantic Ridge (30°41′N). The sample contains microphenocrysts of olivine that display four different rapid-growth morphologies. Comparison of these morphologies with those obtained in dynamic crystallization experiments allows us to constrain the thermal history of the sample. The dendritic morphology (swallowtail, chain and lattice olivine) is directly related to the final quenching during magma–seawater interaction. In contrast, the three other morphologies, namely the complex polyhedral crystal, the closed hopper and the complex swallowtail morphology result from several cycles of cooling–heating (corresponding to a maximum degree of undercooling of 20–25°C) during crystal growth. These thermal variations occurred before eruption and are interpreted to be the result of turbulent convection in a small magmatic body beneath the ridge. The results suggest that the Mid-Atlantic Ridge is underlain by a mush zone that releases batches of liquid during tectonic segregation. Aphyric basalts are emitted during eruptions controlled by the tectonic activity, whereas phyric basalts correspond to small fractions of magma from the mush zone mobilized by reinjections of primitive magmas.