A method for evaluating horizontal well pumping tests

Predicting the future performance of horizontal wells under varying pumping conditions requires estimates of basic aquifer parameters, notably transmissivity and storativity. For vertical wells, there are well-established methods for estimating these parameters, typically based on either the recovery from induced head changes in a well or from the head response in observation wells to pumping in a test well. Comparable aquifer parameter estimation methods for horizontal wells have not been presented in the ground water literature. Formation parameter estimation methods based on measurements of pressure in horizontal wells have been presented in the petroleum industry literature, but these methods have limited applicability for ground water evaluation and are based on pressure measurements in only the horizontal well borehole, rather than in observation wells. This paper presents a simple and versatile method by which pumping test procedures developed for vertical wells can be applied to horizontal well pumping tests. The method presented here uses the principle of superposition to represent the horizontal well as a series of partially penetrating vertical wells. This concept is used to estimate a distance from an observation well at which a vertical well that has the same total pumping rate as the horizontal well will produce the same drawdown as the horizontal well. This equivalent distance may then be associated with an observation well for use in pumping test algorithms and type curves developed for vertical wells. The method is shown to produce good results for confined aquifers and unconfined aquifers in the absence of delayed yield response. For unconfined aquifers, the presence of delayed yield response increases the method error.     

Impacts of Climate-Related Drivers on the Benthic Nutrient Filter in a Shallow Photic Estuary

In shallow photic systems, the benthic filter, including microphytobenthos and denitrifiers, is important in preventing or reducing release of remineralized NH₄ ⁺ to the water column. Its effectiveness can be impacted by climate-related drivers, including temperature and storminess, which by increasing wind and freshwater delivery can resuspend sediment, reduce salinity and deliver nutrients, total suspended solids, and chromophoric dissolved organic matter (CDOM) to coastal systems. Increases in temperature and freshwater delivery may initiate a cascade of responses affecting benthic metabolism with impacts on sediment properties, which in turn regulate nitrogen cycling processes that either sequester (via microphytobenthos), remove (via denitrification), or increase sediment nitrogen (via remineralization, nitrogen fixation, and dissimilatory nitrate reduction to ammonium). We conducted a seasonal study at shallow stations to assess the effects of freshwater inflow, temperature, wind, light, and CDOM on sediment properties, benthic metabolism, nitrogen cycling processes, and the effectiveness of the benthic filter. We also conducted a depth study to constrain seasonally varying parameters such as temperature to better assess the effects of light availability and water depth on benthic processes. Based on relationships observed between climatic drivers and response variables, we predict a reduction in the effectiveness of the benthic filter over the long term with feedbacks that will increase effluxes of N to the water column with the potential to contribute to system eutrophication. This may push shallow systems past a tipping point where trophic status moves from net autotrophy toward net heterotrophy, with new baselines characterized by degraded water quality.     

Europa's atmosphere, gas tori, and magnetospheric implications

A two-dimensional kinetic model calculation for the water group species (H₂O, H₂, O₂, OH, O, H) in Europa's atmosphere is undertaken to determine its basic compositional structure, gas escape rates, and velocity distribution information to initialize neutral cloud model calculations for the most important gas tori. The dominant atmospheric species is O₂ at low altitudes and H₂ at higher altitudes with average day-night column densities of 4.5×¹⁰¹⁴ and 7.7×¹⁰¹³ cm⁻², respectively. H₂ forms the most important gas torus with an escape rate of ∼2×¹⁰²⁷ s⁻¹ followed by O with an escape rate of ∼5×¹⁰²⁶ s⁻¹, created primarily as exothermic O products from O₂ dissociation by magnetospheric electrons. The circumplanetary distributions of H₂ and O are highly peaked about the satellite location and asymmetrically distributed near Europa's orbit about Jupiter, have substantial forward clouds extending radially inward to Io's orbit, and have spatially integrated cloud populations of 4.2×¹⁰³³ molecules for H₂ and 4.0×¹⁰³² atoms for O that are larger than their corresponding populations in Europa's local atmosphere by a factor of ∼200 and ∼1000, respectively. The cloud population for H₂ is a factor of ∼3 times larger than that for the combined cloud population of Io's O and S neutral clouds and provides the dominant neutral population beyond the so-called ramp region at 7.4-7.8 RJ in the plasma torus. The calculated brightness of Europa's O cloud on the sky plane is very dim at the sub-Rayleigh level. The H₂ and O tori provide a new source of europagenic molecular and atomic pickup ions for the thermal plasma and introduce a neutral barrier in which new plasma sinks are created for the cooler iogenic plasma as it is transported radially outward and in which new sinks are created to alter the population and pitch angle distribution of the energetic plasma as it is transported radially inward. The europagenic instantaneous pickup ion rates are peaked at Europa's orbit, dominate the iogenic pickup ion rates beyond the ramp region, and introduce new secondary plasma source peaks in the solution of the plasma transport problem. The H₂ torus is identified as the unknown Europa gas torus that creates both the observed loss of energetic H⁺ ions at Europa's orbit and the corresponding measured ENA production rate for H.     

Nature of the iogenic plasma source in Jupiter's magnetosphere: I. Circumplanetary distribution

Three-dimensional calculations are presented for the circumplanetary nature of the iogenic plasma source (pickup ions produced by electron and charge exchange processes in the plasma torus) created by O and S gases located above Io's exobase in its corona and escaping extended neutral clouds (designated as the “Outer Region”). These calculations are undertaken using neutral cloud models for O and S with realistic incomplete collisional cascade source velocity distributions and rates at Io's exobase and realistic spacetime loss processes in the plasma torus. The resulting spatial distributions for O and S about Jupiter are highly peaked at Io but extend at much lower density levels all about the planet, particularly within Io's orbit where they may play a role in the pitch angle scattering and energy loss of radially inward diffusing energetic electrons for the synchrotron radiation belts of Jupiter, in producing bite-outs in the energy distribution of energetic heavy ions near Io's orbit, and in providing a charge exchange source for energetic neutral atoms (ENAs) detected both near and far from Jupiter. For the iogenic plasma source created by these neutrals, two-dimensional distributions produced by integrating the three-dimensional information along the magnetic field lines are presented for the instantaneous values of the pickup ion rates, the total- and net-mass loading rates, the mass-per-unit-magnetic-flux source rate, the pickup conductivity, the pickup radial current, and the pickup ion power (or energy rate). On the circumplanetary spatial scale, the instantaneous iogenic plasma source is highly peaked about Io's position on its orbit around Jupiter. The degree of orbital asymmetry and its physical origin are discussed, and overall spatially integrated rates are presented. The spatially integrated net-mass loading rate is 154 kg s⁻¹ and the total (electron impact and charge exchange) mass loading rate is 275 kg s⁻¹. Rough minimum estimates are made for the spatially integrated total-mass loading rate created by the “Inner Region” (spatial region below Io's exobase) and are at least ∼1 to 2.5 times larger than that for the Outer Region. Implications of the iogenic plasma source created by the Outer Region and the Inner Region are discussed.     

Shear instability and gravity wave saturation in an asymmetrically stratified jet

Motivated by the mean current and stratification structure associated with the equatorial undercurrent (EUC), we examine the stability and wave propagation characteristics of a highly idealized model flow: the asymmetrically stratified jet. This is a parallel shear flow in which the depth-varying current has the sech² form of a Bickley jet. The stratification has a step function structure: the buoyancy frequency takes uniform values above and below the center of the jet, with the larger value occurring below. The spectrum contains three classes of unstable normal modes. Two are extensions of the sinuous and varicose modes of the unstratified Bickley jet; the third has not been described previously. The asymmetric stratification structure allows instabilities to radiate gravity wave energy from the upper flank of the jet to the lower flank, where it encounters a critical layer. From here, wave energy may be reflected, absorbed or transmitted. Absorption results in wave saturation and momentum transfer to the mean flow, in close analogy with the breaking of orographic gravity waves in the middle atmosphere. Transmission beyond the lower flank may partly account for wave signals observed in the deep equatorial oceans. All of these processes exert zonal forces on the jet that alter its speed and shape. The wave structures and associated fluxes developed by the idealized model are compared with observations of the EUC.     

Thermal expansion of hydrated six-coordinate silicon in thaumasite,Ca<Subscript>3</Subscript>Si(OH)<Subscript>6</Subscript>(CO<Subscript>3</Subscript>)(SO<Subscript>4</Subscript>)·12H<Subscript>2</Subscript>O

Thaumasite, Ca₃Si(OH)₆(CO₃)(SO₄)12H₂O, occurs as a low-temperature secondary alteration phase in mafic igneous and metamorphic rocks, and is recognized as a product and indicator of sulfate attack in Portland cement. It is also the only mineral known to contain silicon in six-coordination with hydroxyl (OH)^? that is stable at ambient P–T conditions. Thermal expansion of the various components of this unusual structure has been determined from single-crystal X-ray structure refinements of natural thaumasite at 130 and 298 K. No phase transitions were observed over this temperature range. Cell parameters at room temperature are: a= 11.0538(6) Å, c=10.4111(8) Å and V=1101.67(10) ų, and were measured at intervals of about 50 K between 130 and 298 K, resulting in mean axial and volumetric coefficients of thermal expansion (×10^?5K^?1); α _a =1.7(1), α _c =2.1(2), and α _V =5.6(2). Although the unit cell and ^VIIICaO₈ polyhedra show significant positive thermal expansion over this temperature range, the silicate octahedron, sulfate tetrahedron, and carbonate group show zero or negative thermal expansion, with α _V (^VISiO₆) = ?0.6 ± 1.1, α _V (^IVSO₄)=?5.8 ± 1.4, and α _R (C–O)= 0.0 ± 1.8 (×10^?5 K^?1). Most of the thermal expansion is accommodated by lengthening of the R(O...O) hydrogen bond distances by on average 5σ, although the hydrogen bonds involving hydroxyl sites on ^VISi expand twice as much as those on molecular water, causing the [Ca₃Si(OH)₆(H₂O)₁₂]⁴⁺ columns to expand in diameter more than they move apart over this temperature range. The average Si–OH bond length of the six-coordinated Si atom 〈R(^VISi–OH)〉 in thaumasite is 1.783(1) Å, being about 0.02 Å (?20σ) shorter than ^VISi–OH in the dense hydrous magnesium silicate, phase D, MgSi₂H₂O₆.     

On‐line identification of non‐linear hysteretic structural systems using a variable trace approach

In this paper, an adaptive on‐line parametric identification algorithm based on the variable trace approach is presented for the identification of non‐linear hysteretic structures. At each time step, this recursive least‐square‐based algorithm upgrades the diagonal elements of the adaptation gain matrix by comparing the values of estimated parameters between two consecutive time steps. Such an approach will enforce a smooth convergence of the parameter values, a fast tracking of the parameter changes and will remain adaptive as time progresses. The effectiveness and efficiency of the proposed algorithm is shown by considering the effects of excitation amplitude, of the measurement units, of larger sampling time interval and of measurement noise. The cases of exact‐, under‐, over‐parameterization of the structural model have been analysed. The proposed algorithm is also quite effective in identifying time‐varying structural parameters to simulate cumulative damage in structural systems. Copyright © 2001 John Wiley & Sons, Ltd.     

Point defects and non-stoichiometry in forsterite

Existing experimental evidence suggests that in pure Mg-forsterite with a fixed Mg/Si ratio, vacancies and interstitials on the Mg sublattice are the dominant intrinsic point defects, while O vacancies and conduction-band electrons are the primary oxygen-deficient defects. Orthopyroxene appears to dissolve readily in pure Mg-forsterite at high temperatures, the change in Mg/Si ratio being accommodated by Mg vacancies and Si interstitials; as a consequence, upper-mantle forsterite may be substantially non-stoichiometric. It is postulated that the amount of Fe that can be oxidized in the forsterite phase is directly related to the excess Si present, either incorporated in the structure, or available from another phase richer in SiO₂ such as pyroxene.     

The spatial distribution of long lived gas clouds emitted by satellites in the outer solar system

Simple models are presented for the spatial distribution of gases emitted by satellites in the outer solar system with emphasis on Io and Titan. The models, valid for long lived species in regions of space outside the gravitational zone of the parent satellite are applied to observed hydrogen and sodium clouds orbiting Jupiter and to an expected hydrogen cloud around Saturn.     

Nonlocal fluxes and Stokes drift effects in the K-profile parameterization

 The K-profile parameterization of upper-ocean mixing is tested and extended using observations and large eddy simulations of upper-ocean response to a westerly windburst. A nonlocal momentum flux term is added, and the amplitude of the nonlocal scalar flux is recalibrated. Parameterizations of Stokes drift effects are added following recent work by McWilliams and Sullivan (2001). These changes allow the parameterization to produce both realistic gradients of momentum and scalars in the nocturnal boundary layer and enhanced mixing during stable conditions. The revised parameterization is expected to produce improved representations of lateral advection and sea-surface temperature in large-scale models. Received: 31 August 2001 / Accepted: 15 December 2001