A sea bottom multichannel hydrophone array

The PUMA (Pull-Up Multichannel Array) is a sea-bottom instrument for remotely recording data from a 12 channel hydrophone array. Its purpose is to achieve (i) denser data coverage, leading to (ii) improved velocity analysis and (iii) multichannel processing of wide angle seismic data collected on the continental shelf. The instrument consists of a 1.2 km array terminating with a pressure case in which 8 FM cassette recorders, a power supply, microprocessor controller and internal clock are housed. It can be pre-programmed to switch on during shot windows for a total of four hours recording time.The PUMA was successfully used in an experiment west of Lewis, Outer Hebrides, U.K. in August–September 1984. We show an example of PUMA data from this experiment. Indications are that the instrument will provide improved constraints on seismic velocities in the lower continental crust and uppermost mantle.     

Phytoplankton dynamics within Gulf Stream intrusions on the southeastern United States continental shelf during summer 1981

During July and August 1981 subsurface intrusion of upwelled nutrient-rich Gulf Stream water was the dominant process affecting temporal and spatial changes in phytoplankton biomass and productivity of the southeastern United States continental shelf between 29 and 32°N latitude. Intruded waters in the study area covered as much as 10¹ km including virtually all of the middle and outer shelf and approximately 50% of the inner shelf area.Within 2 weeks following a large intrusion event in late July, middle shelf primary production and Chl a reached 3 to 4 gC m⁻ d⁻¹ and 75 mg m⁻, respectively. At the peak of the bloom 80% of the water column primary production occurred below the surface mixed-layer, and new primary production (i.e., NO₃-supported) exceeded 90% of the total. Chl a-normalized photosynthetic rates were very high as evidenced by high mean assimilation number (15.5 mg C mg Chl a⁻¹ h⁻¹), high mean α (14 mg C mg Chl a⁻¹ Ein⁻¹ m), and no photoinhibition. As a result of the high photosynthetic rates, mean light-utilization index (Ψ) was 2 to 3 times higher than reported for temperature sub-arctic and arctic waters.The results imply a seasonal (June to August) middle shelf production of 150 g C m⁻¹, about 15% higher than previous estimates of annual production on the middle shelf. Intrusions of the scale we observed in 1981 may not occur every summer. However, when such events do occur, they are by far the most important processes controlling summer phytoplankton dynamics of the middle and outer shelf and of the inner shelf in the southern half of the study area.     

W-shaped occultation signatures: Inference of entwined particle orbits in charged planetary ringlets

Observed W-shaped occultation signatures of certain narrow ringlets in the ring systems of Saturn and Uranus imply a concentration of material near their inner and outer radial edges. A model is proposed where edge bunching is a natural consequence of particles in entwined elliptical orbits, with the same particles alternately defining both edges. While such orbits cross over in radius, collisions would not occur if they have small inclinations, the same fixed argument of periapse ω, and other parameters whereby the particles would “fly in formation” along compressed helical paths relative to the core of the ringlet, which is taken to be a circle in the equatorial plane. For this model to match the observed ring thickness and ringlet widths, orbit inclinations i must be much smaller than their eccentricities e, which themselves would be very small compared to unity. Thus, the meridional cross section of the resultant torus would be a very thin ellipse of thickness proportional to i∣cos ω∣, tilted slightly from the equatorial plane by (i/e)∣sin ω∣ radians. However, gravitational perturbations due to the oblateness of the planet would cause a secular change in ω so that this cross section would collapse periodically to a tilted line, and collisions would then occur. If this collapse could be prevented, the torus could remain in a continuous state of nearly zero viscosity. Stabilization against collapse appears possible due to several remarkable characteristics that are added to the model when the particles are electrically charged. First, because of inherent features of the torus structure, a weak electric force could counter the key effect of the vastly larger oblateness force. Second, because the electric perturbation also affects i, there is a large region in ω,i space where stability against cross-sectional collapse is automatic. For this region, the thickness of the elliptical cross section would expand and contract in concert with the way that the major axis of the ellipse rocks back and forth relative to the equatorial plane. The period of these “rocking and breathing” changes would be from 1 to 3 weeks for a torus in the C ring of Saturn, for example. The electric effects could change considerably without driving the parameters of the torus from the stable domain where cross-sectional collapse does not occur. While specialized and in several important ways still incomplete, the proposed model could account for the W-shaped patterns and explain how very dense ringlets might endure without energy loss due to collisions. It also appears to be capable of explaining the observed sorting of particles by size within a ringlet. Several characteristics of the model suggest definitive tests of its applicability, including its prediction that a nonsymmetrical W-shaped occultation signature could be reversed a half orbit away, and that grazing solar illumination of tilted ringlets might cast shadows that change with time in a prescribed way.     

Process oriented use of geostatistics to analyse creeping para‐glacial features

Complex para‐glacial systems may show signs of destabilization leading to frequent and potentially hazardous debris‐flows. Understanding creeping permafrost displacement over a period of time is therefore crucial for hazard management and risks assessment. This paper presents our methodology for estimating creeping permafrost displacement based upon data derived from various survey methods and demonstrates its relevance on the Glacier Bonnard system in southern Switzerland. Geostatistical processing allowed estimation of the displacement intensity over the area of interest, as well as assessment of the interpolation quality. Although the local measurement network needs to be refined, the results largely improve the understanding of the Bonnard para‐glacial system, though highlighting the need to locally refine the measurement network. In the present case, the destabilized front advances at ~1 m/a and the upstream creeping part at ~0.4 m/a. Variance analysis also provide objective thresholds that could be used to distinguish underlying physical processes. Copyright © 2015 John Wiley & Sons, Ltd.     

A finite-volume ELLAM for three-dimensional solute-transport modeling

A three-dimensional finite-volume ELLAM method has been developed, tested, and successfully implemented as part of the U.S. Geological Survey (USGS) MODFLOW-2000 ground water modeling package. It is included as a solver option for the Ground Water Transport process. The FVELLAM uses space-time finite volumes oriented along the streamlines of the flow field to solve an integral form of the solute-transport equation, thus combining local and global mass conservation with the advantages of Eulerian-Lagrangian characteristic methods. The USGS FVELLAM code simulates solute transport in flowing ground water for a single dissolved solute constituent and represents the processes of advective transport, hydrodynamic dispersion, mixing from fluid sources, retardation, and decay. Implicit time discretization of the dispersive and source/sink terms is combined with a Lagrangian treatment of advection, in which forward tracking moves mass to the new time level, distributing mass among destination cells using approximate indicator functions. This allows the use of large transport time increments (large Courant numbers) with accurate results, even for advection-dominated systems (large Peclet numbers). Four test cases, including comparisons with analytical solutions and benchmarking against other numerical codes, are presented that indicate that the FVELLAM can usually yield excellent results, even if relatively few transport time steps are used, although the quality of the results is problem-dependent.     

Arecibo radar observations of Mars surface characteristics in the northern hemisphere

Mars surface characteristics at and near the Viking Chryse and Tritonis Lacus landing areas were determined by radio scatter using the new 12.6 cm radar at the Arecibo Observatory during 1975–1976. Interpretation of each power spectrum suggests rms surface tilts of 4° at the final A1WNW (47.9°W, 22.5°N) site, 5° near the original A1 site, and 6° between the two. At the back-up site (A2) surface roughness estimates were about 4°. Striking changes in surface texture have been found near the eastern bases of Tharsis Montes and Albor Tholus, each volcanic feature marking the western boundary of very smooth surface units. The roughness sensed at 1 to 100 m scales by radar appears to be relatively independent of the surface units defined at large scale lengths by photogeologists. Radar properties thus provide an additional means by which planetary surfaces may be characterized.     

Saturn's rings: Particle size distributions for thin layer models

The small physical thickness of Saturn's rings requires that radio occultation observations be interpreted using scattering models with limited amounts of multiple scatter. A new model in which the possible order of near-forward scatter is strictly limited allows for the small physical thickness, and can be used to relate Voyager 1 observations of 3.6-and 13-cm wavelength microwave scatter from Saturn's rings to the ring particle size distribution function n(a), for particles with radius 0.001 ≤ a ≤ 20 m. This limited-scatter model yields solutions for particle size distribution functions for eight regions in Saturn's rings, which exhibit approximately inverse-cubic power-law behavior, with large-size cutoffs in particle radius ranging from about 5 m in ring C to about 10 m in parts of ring A. The power-law index is about 3.1 in ring C, about 2.8 in the Cassini division, and increases systematically with radial location in ring A from 2.7 at 2.10R_s to slightly more than 3.0 at 2.24R_s. Corresponding mass densities are 32–43 kg/m² in ring C, 188 kg/m² in the Cassini division, and 244–344 kg/m² in ring A, under the assumption that the material density of the particles is 0.9 g/cm³. These values are a factor of 1 to 2 lower than first-order mass loading estimates derived from resonance phenomena. In view of the uncertainties in the measurements and in the linear density wave model, and the strong arguments for icy particles with specific gravity not greater than about 1, we interpret this discrepancy as being indicative of possible differences in the regions studied, or systematic errors in the interpretation of the scattering results, the density wave phenomena, or some combination of the above.     

A New Natural Secondary Reference Material for Garnet U-Pb Dating by TIMS and LA-ICP-MS

Gem quality andradite-rich garnet (IUC-1), obtained from the Miocene trachyte dome near Ankara city (Turkey), has been identified as a potential natural secondary reference material for U-Pb dating. In this study, U-Pb dating was performed in five different laboratories using isotope dilution TIMS and laser ablation ICP-MS to determine the homogeneity of euhedral garnet crystals. The U-Pb ID-TIMS data for IUC-1 yielded ²⁰⁷Pb/²³⁵U and ²⁰⁶Pb/²³⁸U ages of 20.9 ± 0.4 and 20.6 ± 0.8 Ma respectively, and these values are consistent with U-Pb LA-ICP-MS analyses, in which different garnet crystals yielded ages of 20.8 ± 0.1, 20.7 ± 0.1, 20.7 ± 0.2 and 20.2 ± 0.1 Ma. An andradite (IUC-2) from the Serçeören wollastonite skarn (Turkey) can be used as a secondary reference material provided detailed imaging of the crystals is undertaken. ID-TIMS data yielded ²⁰⁷Pb/²³⁵U and ²⁰⁶Pb/²³⁸U ages of 20.4 ± 0.4 and 20.9 ± 1.0 Ma respectively, and yielded U-Pb ages on different grains of 20.5 ± 0.1, 20.7 ± 1.0, 20.8 ± 1.7 and 20.9 ± 1.6 Ma. The assigned weighted mean age of IUC-1 (20.4 ± 0.5 Ma, 2s) is proposed as a 2023 reference value. IUC-1 garnet is expected to contribute significantly to rapidly developing garnet geochronology in the near future.     

Modeling effects of multinode wells on solute transport

Long-screen wells or long open boreholes with intraborehole flow potentially provide pathways for contaminants to move from one location to another in a ground water flow system. Such wells also can perturb a flow field so that the well will not provide water samples that are representative of ground water quality a short distance away from the well. A methodology is presented to accurately and efficiently simulate solute transport in ground water systems that include wells longer than the grid spacing used in a simulation model of the system and hence are connected to multiple nodes of the grid. The methods are implemented in a MODFLOW-compatible solute-transport model and use MODFLOW's Multi-Node Well Package but are generic and can be readily implemented in other solute-transport models. For nonpumping multinode wells (used to simulate open boreholes or observation wells, for example) and for low-rate pumping wells (in which the flow between the well and the ground water system is not unidirectional), a simple routing and local mixing model was developed to calculate nodal concentrations within the borehole. For high-rate pumping multinode wells (either withdrawal or injection, in which flow between the well and the ground water system is in the same direction at all well nodes), complete and instantaneous mixing in the wellbore of all inflows is assumed.     

The joint effect of rank and grain size on the flotation of Australian bituminous coals

Laboratory batch flotation tests were carried out to determine the effect of grain size on the recovery of four Australian bituminous coals of high and low rank. The flotation concentrates, taken over successive time intervals, and the tailings, were sized and the recovery-size curves plotted. A method was developed for determining the proportions of liberated and composite grains in the various flotation products. Flotation tests were carried out in the presence of methylisobutylcarbinol (MIBC) and MIBC + kerosene. The shapes of the recovery-size curves depended on the rank of the coal and whether or not kerosene was present. The variation of the recovery with size was explained in terms of the joint effect of rank and size on the floatability of composite grains.