Simulated impacts of artificial groundwater recharge and discharge on the source area and source volume of an Atlantic Coastal Plain stream, Delaware, USA

A numerical groundwater-flow model was used to characterize the source area and volume of Phillips Branch, a baseflow-dominated stream incising a highly permeable unconfined aquifer on the low relief Delmarva Peninsula, USA. Particle-tracking analyses indicate that the source area (5.51 km²) is ～20% smaller than the topographically defined watershed (6.85 km²), and recharge entering ～37% of the surface watershed does not discharge to Phillips Branch. Groundwater residence time within the source volume ranges from a few days to almost 100 years, with 95% of the volume “flushing” within 50 years. Artificial discharge from groundwater pumping alters the shape of the source area and reduces baseflow due to the interception of stream flow paths, but has limited impacts on the residence time of groundwater discharged as baseflow. In contrast, artificial recharge from land-based wastewater disposal substantially reduces the source area, lowers the range in residence time due to the elimination of older flow paths to the stream, and leads to increased discharge to adjacent surface-water bodies. This research suggests that, in this and similar hydrogeologic settings, the “watershed” approach to water-resource management may be limited, particularly where anthropogenic stresses alter the transport of soluble contaminants through highly permeable unconfined aquifers.     

3D model of fracture zones at Soultz-sous-Forêts based on geological data,image logs,induced microseismicity and vertical seismic profiles

This study presents a new deterministic 3D model of the fracture zones observed in the granitic reservoir of the Soultz European geothermal project. The major fracture zones encountered around 6 wells (4550, EPS1, GPK1, GPK2, GPK3 and GPK4) consist in 53 main structures that are located and characterized in terms of size and orientation: 39 fracture zones, 8 microseismic structures and 6 structures derived from vertical seismic profiles are represented in the 3D model using Discrete Fracture Network tools of the gOcad modelling platform (Paradigm?, Earth Decision?). This work illustrates the complexity of 3D fracture zone correlation and interpretation in crystalline rock masses characterized at meter scale (borehole) and at the reservoir scale (kilometer) thanks to geophysical imaging techniques.     

Modelling the mean and turbulent structure of the summertime Arctic cloudy boundary layer

This paper addresses the problem of modelling the summertime Arctic cloudy boundary layer. Specifically we consider the problem of multi-layered clouds in the boundary layer that includes the decoupling of the turbulence between upper and lower clouds. A high-resolution one-dimensional model with second-order turbulence closure and spectral radiative transfer is used to simulate a case study that was obtained during the 1980 Arctic Stratus Experiment. The effects of radiation, large-scale vertical motion and drizzle are investigated in sensitivity studies. Results of this study show that radiative transfer is important to the maintenance of the multiple cloud layers, and suggest that weak rising vertical motion is the most favorable situation to maintain two separate cloud layers.     

Field measurements of snow-drift threshold and mass fluxes,and related model simulations

Field measurements were carried out to calculate the threshold friction velocity for snow saltation, and mass fluxes during snow drift. The wind was measured in three components by an ultrasonic anemometer, and the mass fluxes were determined using an optical sensor (snow particle counter), acoustic sensors (Flowcapt) and mechanical traps. The threshold friction velocity was found to be correlated to the grain size (R²=0.75). The mass flux measurements were compared with numerical simulations of snow drift, and it was demonstrated that the maximum snow transport takes place at shear stress values of roughly two times the average shear stress over 20 min. By implementing a probability distribution for the shear stress the mass flux was simulated with only the mean measured value of the shear stress as input. This procedure enables the future use of the numerical model for operational applications.     

Assessing the Goodness-of-Fit of Statistical Distributions When Data Are Grouped

Modeling statistical distributions of phenomena can be compromised by the choice of goodness-of-fit statistics. The Pearson chi-square test is the most commonly used test in the geosciences, but the lesser known empirical distribution function (EDF) statistics should be preferred in many test situations. Using a data set from geomorphology, the Anderson–Darling test for grouped exponential distributions is employed to illustrate ease of use and statistical advantages of this EDF test. Attention to the issues discussed will result in more informed statistic selection and increased rigor in the identification of distribution functions that describe random variables.     

Mapping Curvilinear Structures with Local Anisotropy Kriging

Interpolating geo-data with curvilinear structures using geostatistics is often disappointing. Channels, for example, become disconnected sets of lakes when interpolated from point data. In order to improve the interpolation of geological structures (e.g., curvilinear structures), we present a new form of kriging, local anisotropy kriging (LAK). Local anisotropy kriging combines a gradient algorithm from image analysis with kriging in an iterative way. After an initial standard kriging interpolation, the gradient algorithm determines the local anisotropy for each cell in the grid using a search area around the cell. Subsequently, kriging is carried out with the spatially varying anisotropy. The anisotropy calculation and subsequent kriging steps will then succeed until the result is satisfactory in the way of reproducing the curvilinear structures. Depending on the size of the search area more or less detail in the geological structures can be reproduced with LAK. Using test examples we show that LAK interpolates data with curvilinear structures more realistically than standard kriging. In a real world case, using bathymetric data of the Oosterschelde estuary, LAK also proves to be quantitatively superior to standard kriging. Absolute interpolation errors are decreased by 23%. Local anisotropy kriging only uses information from point data, which makes the method very objective, it only presents “what the data can tell.”     

Effects of tributyltin on activity and burrowing behavior of the fiddler crab,Uca pugilator

Fiddler crabs,Uca pugilator, collected from the field showed no avoidance to burrowing in 1 μg g^?1 tributyltin (TBT)-contaminated sand held in laboratory trays. Treatment ofU. pugilator with levels of tributyltin as low as 0.5 μg l^?1 for 1–3 wk resulted in an acceleration of the righting reflex in females, indicative of hyperactivity. Crabs of both sexes exhibited a reduction in burrowing activity, as measured by the number of burrows dug at 15 and 60 min after release into laboratory trays containing sand, and by the number of crabs within burrows at those times. The reduction in burrowing activity was not dose-dependent at concentrations of 0.5 to 50 μg l^?1, and was unchanged between one and three weeks of exposure. Since fiddler crabs are dependent on burrows for many aspects of their lives, the reduction in burrowing behavior, should it occur in nature, would have serious consequences for the species.     

Astronomical imaging with InSb arrays

Ten years ago, Forrest presented the first astronomical images with an SBRC 32×32 InSb array camera at the first NASA-Ames Infrared Detector Technology Workshop. Soon after, SBRC began development of 58×62 InSb arrays, both for ground-based astronomy and for SIRTF. By the time of the 1987 Hilo workshop Ground-based Astronomical Observations with Infrared Array Detectors astronomical results from cameras based on SBRC 32×32 and 58×62 InSb arrays, a CE linear InSb array, and a French 32×32 InSb CID array were presented. And at the Tucson 1990 meeting Astrophysics with Infrared Arrays, it was clear that this new technology was no longer the province of IR pundits, but provided a tool for all astronomers. At this meeting, the first astronomical observations with SBRC's new, gateless passivation 256×256 InSb arrays will be presented: they perform spectacularly!In this review, I can only broadly brush on the interesting science completed with InSb array cameras. Because of the broad wavelength coverage (1–5.5 m) of InSb, and the extremely high performance levels throughout the band, InSb cameras are used not only in the near IR, but also from 3–5.5 m, where unique science is achieved. For example, the point-like central engines of AGNs are delineated at L and M, and Br and 3.29 m dust emission images of galactic and extragalactic objects yield excitation conditions. Examples of imaging spectroscopy, high spatial resolution imaging, as well as deep, broad-band imaging with InSb cameras at this meeting illustrate the power of InSb array cameras.