In situ estimation of the effective chemical diffusion coefficient of a rock matrix in a fractured aquifer

An in situ method of estimating the effective diffusion coefficient for a chemical constituent that diffuses into the primary porosity of a rock is developed by abruptly changing the concentration of the dissolved constituent in a borehole in contact with the rock matrix and monitoring the time-varying concentration. The experiment was conducted in a borehole completed in mudstone on the campus of the University of the Free State in Bloemfontein, South Africa. Numerous tracer tests were conducted at this site, which left a residual concentration of sodium chloride in boreholes that diffused into the rock matrix over a period of years. Fresh water was introduced into a borehole in contact with the mudstone, and the time-varying increase of chloride was observed by monitoring the electrical conductivity (EC) at various depths in the borehole. Estimates of the effective diffusion coefficient were obtained by interpreting measurements of EC over 34 d. The effective diffusion coefficient at a depth of 36 m was approximately 7.8×10^?6 m²/d, but was sensitive to the assumed matrix porosity. The formation factor and mass flux for the mudstone were also estimated from the experiment.     

Radar observations of asteroid 1580 Betulia

Radar observations of the asteroid 1580 Betulia, made at a wavelength of 12.6 cm, show a mean radar cross section of 2.2 ± 0.8 km² and a total spectral bandwidth of 26.5 ± 1.5 Hz. Combining our bandwidth measurements with the optically determined rotation period sets a lower limit to the asteroid's radius of 2.9 ± 0.2 km.     

Galilean satellites: 1976 radar results

Radar observations of the Galilean satellites, made in late 1976 using the 12.6-cm radar system of the Arecibo Observatory, have yielded mean geometric albedos of 0.04 ± , 0.69 ± 0.17, 0.37 ± 0.09, and 0.15 ± 0.04, for Io, Europa, Ganymede, and Callisto, respectively. The albedo for Io is about 40% smaller than that obtained approximately a year earlier, while the albedos for the outer three satellites average about 70% larger than the values previously reported for late 1975, raising the possibility of temporal variation. Very little dependence on orbital phase is noted; however, some regional scattering inhomogeneities are seen on the outer three satellites. For Europa, Ganymede, and Callisto, the ratios of the echo received in one mode of circular polarization to that received in the other were: 1.61 ± 0.20 1.48 ± 0.27, and 1.24 ± 0.19, respectively, with the dominant component having the same sence of circularity as that transmitted. This behavior has not previously been encountered in radar studies of solar system objects, whereas the corresponding observations with linear polarization are “normal.” Radii determined from the 1976 radar data for Europa and Ganymede are: 1530 ± 30 and 2670 ± 50 km, in fair agreement with the results from the 1975 radar observations and the best recent optical determinations. Doppler shifts of the radar echoes, useful for the improvement of the orbits of Jupiter and some of the Galilean satellites, are given for 12 nights in 1976 and 10 nights in 1975.     

Relative effects of physical and biological processes on nutrient and phytoplankton dynamics in a shallow estuary after a storm event

The effects of advection, dispersion, and biological processes on nitrogen and phytoplankton dynamics after a storm event in December 2002 are investigated in an estuary located on the northern New South Wales coast, Australia. Salinity observations for 16 d after the storm are used to estimate hydrodynamic transports for a one-dimensional box model. A biological model with nitrogen limited phytoplankton growth, mussel grazing, and a phytoplankton mortality term is forced by the calculated transports. The model captured important aspects of the temporal and spatial dynamics of the bloom. A quantitative analysis of hydrodynamic and biological processes shows that increased phytoplankton biomass due to elevated nitrogen loads after the storm was not primarily regulated by advection or dispersion in spite of an increase in river flow from <1 to 928×10³ m³ d⁻¹. Of the dissolved nitrogen that entered the surface layer of the estuary in the 16 d following the storm event, the model estimated that 28% was lost through exchange with the ocean or bottom layers, while 15% was removed by the grazing of just one mussel species,Xenostrobus securis, on phytoplankton, and 50% was lost through other biological phytoplankton loss processes.X. securis grazing remained an important loss process even when the estimated biological parameters in the model were varied by factors of ± 2. The intertidal mangrove pneumatophore habitat ofX. securis allows filtering of the upper water column from the lateral boundaries when the water column is vertically stratified, exerting top-down control on phytoplankton biomass.     

Cautions and Suggestions for Geochemical Sampling in Fractured Rock

Collecting water samples for geochemical analyses in open bedrock boreholes or in discrete intervals of boreholes intersected by multiple fractures is likely to yield ambiguous results for ground water chemistry because of the variability in the transmissivity, storativity, and hydraulic head of fractures intersecting the borehole. Interpreting chemical analyses of water samples collected in bedrock boreholes requires an understanding of the hydraulic conditions in the borehole under the ambient flow regime in the aquifer as well as during sampling. Pumping in open boreholes, regardless of the pumping rate and the location of the pump intake, first draws water from the borehole and then from fractures intersecting the borehole. The time at which the volumetric rate of water entering the borehole from fractures is approximately equal to the pumping rate can be identified by monitoring the logarithm of drawdown in the borehole as a function of the logarithm of time. Mixing of water entering the borehole from fractures with water in the borehole must be considered in estimating the time at which the pump discharge is representative of aquifer water. In boreholes intersected by multiple fractures, after the contribution from the borehole volume has diminished, the contribution of fractures to the pump discharge will be weighted according to their transmissivity, regardless of the location of the pump intake. This results in a flux-averaged concentration in the pump discharge that is biased by the chemical signature of those fractures with the highest transmissivity. Under conditions where the hydraulic head of fractures varies over the length of the borehole, open boreholes will be subject to ambient flow in the water column in the borehole. In some instances, the magnitude of the ambient flow may be similar to the designated pumping rate for collecting water samples for geochemical analyses. Under such conditions, the contributions to the pump discharge from individual fractures will be a function not only of the transmissivity of the fractures, but also of the distribution of hydraulic head in fractures intersecting the borehole. To reduce or eliminate the deleterious effects of conducting geochemical sampling in open boreholes, a straddle-packer apparatus that isolates a single fracture or a series of closely spaced fractures is recommended. It is also recommended that open boreholes be permanently outfitted with borehole packers or borehole liners in instances where maintaining the hydraulic and chemical stratification in the aquifer is of importance. In a field example, a comparison of results from sampling in an open borehole and in discrete intervals of the same borehole showed dramatic differences in the concentrations of chemical constituents in the water samples, even though chemical field parameters stabilized prior to both open borehole and discrete interval sampling.     

Effective elastic properties of randomly fractured soils: 3D numerical experiments

This paper is concerned with numerical tests of several rock physical relationships. The focus is on effective velocities and scattering attenuation in 3D fractured media. We apply the so‐called rotated staggered finite‐difference grid (RSG) technique for numerical experiments. Using this modified grid, it is possible to simulate the propagation of elastic waves in a 3D medium containing cracks, pores or free surfaces without applying explicit boundary conditions and without averaging the elastic moduli. We simulate the propagation of plane waves through a set of randomly cracked 3D media. In these numerical experiments we vary the number and the distribution of cracks. The synthetic results are compared with several (most popular) theories predicting the effective elastic properties of fractured materials. We find that, for randomly distributed and randomly orientated non‐intersecting thin penny‐shaped dry cracks, the numerical simulations of P‐ and S‐wave velocities are in good agreement with the predictions of the self‐consistent approximation. We observe similar results for fluid‐filled cracks. The standard Gassmann equation cannot be applied to our 3D fractured media, although we have very low porosity in our models. This is explained by the absence of a connected porosity. There is only a slight difference in effective velocities between the cases of intersecting and non‐intersecting cracks. This can be clearly demonstrated up to a crack density that is close to the connectivity percolation threshold. For crack densities beyond this threshold, we observe that the differential effective‐medium (DEM) theory gives the best fit with numerical results for intersecting cracks. Additionally, it is shown that the scattering attenuation coefficient (of the mean field) predicted by the classical Hudson approach is in excellent agreement with our numerical results.     

A revised evolutionary history of armadillos (Dasypus) in North America based on ancient mitochondrial DNA

The large, beautiful armadillo, Dasypus bellus, first appeared in North America about 2.5 million years ago, and was extinct across its southeastern US range by 11 thousand years ago (ka). Within the last 150 years, the much smaller nine‐banded armadillo, D. novemcinctus, has expanded rapidly out of Mexico and colonized much of the former range of the beautiful armadillo. The high degree of morphological similarity between these two species has led to speculation that they might be a single, highly adaptable species with phenotypical responses and geographical range fluctuations resulting from environmental changes. If this is correct, then the biology and tolerance limits for D. novemcinctus could be directly applied to the Pleistocene species, D. bellus. To investigate this, we isolated ancient mitochondrial DNA from late Pleistocene‐age specimens of Dasypus from Missouri and Florida. We identified two genetically distinct mitochondrial lineages, which most likely correspond to D. bellus (Missouri) and D. novemcinctus (Florida). Surprisingly, both lineages were isolated from large specimens that were identified previously as D. bellus. Our results suggest that D. novemcinctus, which is currently classified as an invasive species, was already present in central Florida around 10 ka, significantly earlier than previously believed.     

Very high CO2 exchange fluxes at the peak of the rainy season in a West African grazed semi-arid savanna ecosystem

Africa is a sink of carbon, but there are large gaps in our knowledge regarding the CO₂ exchange fluxes for many African ecosystems. Here, we analyse multi-annual eddy covariance data of CO₂ exchange fluxes for a grazed Sahelian semi-arid savanna ecosystem in Senegal, West Africa. The aim of the study is to investigate the high CO₂ exchange fluxes measured at the peak of the rainy season at the Dahra field site: gross primary productivity and ecosystem respiration peaked at values up to ?48 μmol CO₂ m^?2 s^?1 and 20 μmol CO₂ m^?2 s^?1, respectively. Possible explanations for such high fluxes include a combination of moderately dense herbaceous C4 ground vegetation, high soil nutrient availability and a grazing pressure increasing the fluxes. Even though the peak net CO₂ uptake was high, the annual budget of ?229 ± 7 ± 49 g C m^?2 y^?1 (±random errors ± systematic errors) is comparable to that of other semi-arid savanna sites due the short length of the rainy season. An inter-comparison between the open-path and a closed-path infrared sensor indicated no systematic errors related to the instrumentation. An uncertainty analysis of long-term NEE budgets indicated that corrections for air density fluctuations were the largest error source (11.3% out of 24.3% uncertainty). Soil organic carbon data indicated a substantial increase in the soil organic carbon pool for the uppermost .20 m. These findings have large implications for the perception of the carbon sink/source of Sahelian ecosystems and its response to climate change.