A modified light transmission visualization method for DNAPL saturation measurements in 2-D models

In this research, a light transmission visualization (LTV) method was used to quantify dense non-aqueous phase liquids (DNAPL) saturation in two-dimensional (2-D), two fluid phase systems. The method is an expansion of earlier LTV methods and takes into account both absorption and refraction light theories. Based on this method, DNAPL and water saturations can rapidly be obtained point wise across sand-packed 2-D flow chambers without the need to develop a calibration curve. A single point calibration step is, however, needed when dyed DNAPL is used to account for the change in the transmission factor at the dyed DNAPL–water interface. The method was applied to measure, for the first time, undyed DNAPL saturation in small 2-D chambers. Known amounts of DNAPL, modeled by tetrachloroethylene (PCE), were added to the chamber and these amounts were compared to results obtained by this LTV method. Strong correlation existed between results obtained based on this method and the known PCE amounts with an R² value of 0.993. Similar experiments conducted using dyed PCE showed a stronger correlation between results obtained by this LTV method and the known amounts of dyed PCE added to the chamber with an R² value of 0.999. The method was also used to measure dyed PCE saturation in a large 2-D model following sparging experiments. Results obtained from image analyses following each sparging event were compared to results obtained by two independent techniques, namely gas chromatography–mass spectrometry (GC/MS) analyses and carbon column extraction. There was a good agreement between the results obtained by this LTV method and those obtained by the two independent techniques when experiments were carried out under stable light source conditions and errors in mass balance were minor. The method presented here can be expanded to measure fluid contents in three fluid phase systems and provide a non-destructive, non-intrusive tool to investigate changes in DNAPL architecture and flow characteristics in laboratory experiments.     

The fall of the St-Robert meteorite

Abstract The St-Robert (Québec, Canada) meteorite shower occurred on 1994 June 15 at 0h02m UT accompanied by detonations audible for >200 km from the fireball endpoint. The fireball was recorded by visual observers in Vermont, New York State, New Hampshire, Québec and Ontario as well as by optical and infrared sensors in Earth-orbit. Penetration to an altitude of 36 km occurred ～60 km to the northeast of Montreal, where the bolide experienced several episodes of fragmentation. A total of 20 fragments of this H5 chondrite, comprising a total mass of 25.4 kg, were recovered in an ellipse measuring 8 × 3.5 km. One fragment of the shower partially penetrated the aluminum roof of a shed. Interpretation of the visual and satellite data suggests that the fireball traveled from south-southwest to north-northeast, with a slope from the horizontal of 55°–61°. A statistical evaluation of the likely heliocentric orbits for the body prior to collision with the Earth, coupled with theoretical modeling of the entry, suggests an entry velocity in the range of 12.7–13.3 km/s; the meteoroid had moved in a low-inclination orbit, with orbital perihelion located extremely close to the Earth's orbit. From satellite optical data, it is found that the photometric mass consumed during the largest detonation is ～1200 kg. Estimation of the amplitude of the acoustic signal detected by the most distant observer yields a source energy near 0.5 kt TNT equivalent energy, which corresponds to a mass of order 10 metric tonnes. This measure is uncertain to approximately one order of magnitude. Theoretical modeling of the entry of the object suggests a mass near 1600 kg. Cosmogenic radionuclide activities constrain the lower initial mass to be ～700 kg with an upper limit near 4000 kg. Seismic data possibly associated with the fireball suggest extremely poor coupling between the airwave and the ground. The total mass estimated to have reached the ground is ～100 kg (in material comprising >55 g fragments), while the preatmospheric mass is found to be most probably in the range of 1200–2000 kg.     

Distribution and chemistry of manganese,iron, and suspended particulates in Orca Basin

The intense halocline and redoxcline in the Orca Basin, northwest Gulf of Mexico, induce dramatic water column profiles for manganese, iron, and suspended particulates. Within a 17 m interval, the salinity of the basin water increases from 66 to 260 & permil and dissolved oxygen decreases to zero. Midway through this transition zone, concentrations of suspended matter peak at 900 g/liter. Dissolved iron and manganese concentrations in the anoxic brine increase from oceanic values to maxima of 1.6 and 22 mg/liter, respectively. Upward migration of dissolved manganese from the brine leads to production of manganese-rich particles in the slightly oxygenated overlying water.     

Mapping and analysing virtual outcrops

Laser scanning is a very efficient way to generate realistic, high-resolution digital models of 3-D geological outcrops. This paper discusses the methodologies involved in the creation and analysis of virtual outcrops, based on laser scanner data. The visualisation of the laser scanner data as a photorealistic 3-D object is described. Geological features picked out on the virtual outcrop (e.g. fractures, faults or bedding planes) can be extrapolated outward, into space, and inward, into the subsurface, using tension surfaces.Electronic Supplementary Material Supplementary material is available in the online version of this article at Reviewed by: J.D. Clemens, D. Yuen     

The influence of autocorrelation on the ability to detect trend in hydrological series

This study investigated using Monte Carlo simulation the interaction between a linear trend and a lag‐one autoregressive (AR(1)) process when both exist in a time series. Simulation experiments demonstrated that the existence of serial correlation alters the variance of the estimate of the Mann–Kendall (MK) statistic; and the presence of a trend alters the estimate of the magnitude of serial correlation. Furthermore, it was shown that removal of a positive serial correlation component from time series by pre‐whitening resulted in a reduction in the magnitude of the existing trend; and the removal of a trend component from a time series as a first step prior to pre‐whitening eliminates the influence of the trend on the serial correlation and does not seriously affect the estimate of the true AR(1). These results indicate that the commonly used pre‐whitening procedure for eliminating the effect of serial correlation on the MK test leads to potentially inaccurate assessments of the significance of a trend; and certain procedures will be more appropriate for eliminating the impact of serial correlation on the MK test. In essence, it was advocated that a trend first be removed in a series prior to ascertaining the magnitude of serial correlation. This alternative approach and the previously existing approaches were employed to assess the significance of a trend in serially correlated annual mean and annual minimum streamflow data of some pristine river basins in Ontario, Canada. Results indicate that, with the previously existing procedures, researchers and practitioners may have incorrectly identified the possibility of significant trends. Copyright © Environment Canada. Published by John Wiley & Sons, Ltd.     

The use of geometric and gamma-related distributions for frequency analysis of water deficit

This paper presents an approach to perform statistical frequency analysis of water deficit duration and severity using respectively the geometric and exponential distributions. Monthly mean water discharges are compared to a given threshold and classified in two mutually exclusive ways. This leads to a two state random variable such that: a success represents the absence of a water deficit event (mean monthly discharge exceeds threshold), and a failure, a water deficit event (mean monthly discharge is below threshold). If we suppose that this random variable gives rise to a Markov process of order 1, then the duration of a water deficit event X (consecutive months in deficit) will have a geometric distribution. In turn, the summation of discharges in deficit will give the severity of a water deficit event which can be represented by a one-parameter exponential distribution. The threshold or base level is taken as a percentile of the observed mean discharges of a given month. This base level, which varies from month to month, can be viewed as the limit of an acceptable deficit (or energetic failure) associated to a given empirical probability of being in deficit. The second step of the approach is to estimate the value of the parameter for each distribution using the maximum likelihood method. Expressions for the estimator of a given percentile, , as well as its variance are deduced. Finally, the presented models are applied to observed data.     

Spherical harmonic synthesis and least squares computations in satellite gravity gradiometry

The computational requirements in the simulations of geopotential estimation from satellite gravity gradiometry are discussed. Fast algorithms for spherical harmonic synthesis and least squares accumulation on a vectorizing supercomputers are presented. Using these methods, in a test case estimation of 2595 coefficients of a degree and order 50 gravity field, sustained program execution speeds of 275 Mflops (87 % peak machine speed) on a single processor of a CRAY Y-MP were achieved, with spherical harmonics computation accounting for less than 1 % of total cost. From the results, it appears that brute-force estimation of a degree and order 180 field would require 537 Million Words of memory and 85 hours of CPU time, assuming mission duration of 1 month, and execution speed of 1 Gflops. Both memory size and execution speed requirements are within the capabilities of modern multi-processor supercomputers.