A modelling study of hyporheic exchange pattern and the sequence,size, and spacing of stream bedforms in mountain stream networks,Oregon, USA

Studies of hyporheic exchange flows have identified physical features of channels that control exchange flow at the channel unit scale, namely slope breaks in the longitudinal profile of streams that generate subsurface head distributions. We recently completed a field study that suggested channel unit spacing in stream longitudinal profiles can be used to predict the spacing between zones of upwelling (flux of hyporheic water into the stream) and downwelling (flux of stream water into the hyporheic zone) in the beds of mountain streams. Here, we use two‐dimensional groundwater flow and particle tracking models to simulate vertical and longitudinal hyporheic exchange along the longitudinal axis of stream flow in second‐, third‐, and fourth‐order mountain stream reaches. Modelling allowed us to (1) represent visually the effect that the shape of the longitudinal profile has on the flow net beneath streambeds; (2) isolate channel unit sequence and spacing as individual factors controlling the depth that stream water penetrates the hyporheic zone and the length of upwelling and downwelling zones; (3) evaluate the degree to which the effects of regular patterns in bedform size and sequence are masked by irregularities in real streams. We simulated hyporheic exchange in two sets of idealized stream reaches and one set of observed stream reaches. Idealized profiles were constructed using regression equations relating channel form to basin area. The size and length of channel units (step size, pool length, etc.) increased with increasing stream order. Simulations of hyporheic exchange flows in these reaches suggested that upwelling lengths increased (from 2·7 m to 7·6 m), and downwelling lengths increased (from 2·9 m to 6·0 m) with increase in stream order from second to fourth order. Step spacing in the idealized reaches increased from 5·3 m to 13·7 m as stream size increased from second to fourth order. Simulated downwelling lengths increased from 4·3 m in second‐order streams to 9·7 m in fourth‐order streams with a POOL–RIFFLE–STEP channel unit sequence, and increased from 2·5 m to 6·1 m from second‐ to fourth‐order streams with a POOL–STEP–RIFFLE channel unit sequence. Upwelling lengths also increased with stream order in these idealized channels. Our results suggest that channel unit spacing, size, and sequence are all important in determining hyporheic exchange patterns of upwelling and downwelling. Though irregularities in the size and spacing of bedforms caused flow nets to be much more complex in surveyed stream reaches than in idealized stream reaches, similar trends emerged relating the average geomorphic wavelength to the average hyporheic wavelength in both surveyed and idealized reaches. This article replaces a previously published version (Hydrological Processes, 19 (17), 2915–2929 (2005) [ DOI:10.1002/hyp.5790 ]. See also retraction notice DOI:10.1002/hyp.6350 Copyright © 2006 John Wiley & Sons, Ltd.     

Separation of up-going and down-going wave fields of vertical cable data

The vertical cable method for acquiring and processing pre-stack 3-D marine seismic data is based on the technology developed by the US Navy for antisubmarine warfare. In order to achieve the maximum utilization of vertical cable field data, a new separation method of the up-going and down-going wave fields of the vertical cable data processing was developed in this paper, which is different from the separation of the down-going and up-going wave fields of normal VSP data processing. In tests with synthetic modeling data andactual field data, this newly developed method performs well and is also computationally simpler without pre-assumption conditions.     

The Gas Production Rate and Coma Structure of Comet C/1995 O1 (Hale–Bopp)

The University of Wisconsin–Madison and NASA–Goddard conducted acomprehensive multi-wavelength observing campaign of coma emissionsfrom comet Hale–Bopp, including OH 3080 Å, [O I] 6300 Å H₂O⁺ 6158 Å, H Balmer-α 6563 Å, NH₂ 6330 Å, [C I] 9850 ÅCN 3879 Å, C₂ 5141 Å, C₃ 4062 Å,C I 1657 Å, and the UV and optical continua. In thiswork, we concentrate on the results of the H₂O daughter studies.Our wide-field OH 3080 Å measured flux agrees with other, similarobservations and the expected value calculated from published waterproduction rates using standard H₂O and OH photochemistry.However, the total [O I] 6300 Å flux determined spectroscopically overa similar field-of-view was a factor of 3-4 higher than expected.Narrow-band [O I] images show this excess came from beyond theH₂O scale length, suggesting either a previously unknown source of[O I] or an error in the standard OH + ν→ O(¹ D) + H branching ratio. The Hale–Bopp OH and[O I] distributions, both of which were imaged tocometocentric distances >1 × 10⁶ km, were more spatiallyextended than those of comet Halley (after correcting for brightnessdifferences), suggesting a higher bulk outflow velocity. Evidence ofthe driving mechanism for this outflow is found in the Hα lineprofile, which was narrower than in comet Halley (though likelybecause of opacity effects, not as narrow as predicted by Monte-Carlomodels). This is consistent with greater collisional coupling betweenthe suprathermal H photodissociation products and Hale–Bopp's densecoma. Presumably because of mass loading of the solar wind by ionsand ions by the neutrals, the measured acceleration of H₂O⁺ downthe ion tail was much smaller than in comet Halley. Tailwardextensions in the azimuthal distributions of OH 3080 Å,[O I], and [C I] , as well as a Doppler asymmetry in the[O I] line profile, suggest ion-neutral coupling. While thetailward extension in the OH can be explained by increased neutralacceleration, the [O I] 6300 Å and [C I] 9850 Å emissions show 13%and >200% excesses in this direction (respectively), suggesting anon-negligible contribution from dissociative recombination of CO⁺and/or electron collisional excitation. Thus, models including theeffects of photo- and collisional chemistry are necessary for the fullinterpretation of these data.     

The Tiny-Scale Atomic Structure: Gas Cloudlets or Scintillation Phenomenon?

We present preliminary results from the recent multi-epoch HI absorption measurements toward several pulsars using the Arecibo telescope. We do not find significant variations in optical depth profiles over periods of 4 months and 9 years. The upper limits on the optical depth variations in directions of B0823+26 and B1133+16 are similar on scales of 10–20, 350 and 500 AU. The large number of non detections of the tiny scale atomic structure suggests that the AU-sized structure is not ubiquitous in the interstellar medium.     

Tectonic histories between Alba Patera and Syria Planum, Mars

Syria Planum and Alba Patera are two of the most prominent features of magmatic-driven activity identified for the Tharsis region and perhaps for all of Mars. In this study, we have performed a Geographic Information System-based comparative investigation of their tectonic histories using published geologic map information and Mars Orbiter Laser Altimetry (MOLA) data. Our primary objective is to assess their evolutional histories by focusing on their extent of deformation in space and time through stratigraphic, paleotectonic, topographic, and geomorphologic analyses. Though there are similarities among the two prominent features, there are several distinct differences, including timing deformational extent, and tectonic intensity of formation. Whereas Alba Patera displays a major pulse of activity during the Late Hesperian/Early Amazonian, Syria Planum is a long-lived center that displays a more uniform distribution of simple graben densities ranging from the Noachian to the Amazonian, many of which occur at greater distances away from the primary center of activity. The histories of the two features presented here are representative of the complex, long-lived evolutional history of Tharsis.     

Quantifying the Effects of Commercial Clam Aquaculture on C and N Cycling: an Integrated Ecosystem Approach

Increased interest in using bivalve cultivation to mitigate eutrophication requires a comprehensive understanding of the net carbon (C) and nitrogen (N) budgets associated with cultivation on an ecosystem scale. This study quantified C and N processes related to clam (Mercenaria mercenaria) aquaculture in a shallow coastal environment (Cherrystone Inlet, VA) where the industry has rapidly increased. Clam physiological rates were compared with basin-wide ecosystem fluxes including primary production, benthic nutrient regeneration, and respiration. Although clam beds occupy only 3 % of the ecosystem’s surface area, clams filtered 7–44 % of the system’s volume daily, consumed an annual average of 103 % of the phytoplankton production, creating a large flux of particulate C and N to the sediments. Annually, N regenerated and C respired by clam and microbial metabolism in clam beds were ～3- and ～1.5-fold higher, respectively, than N and C removed through harvest. Due to the short water residence time, the low watershed load, and the close vicinity of clam beds to the mouth of Cherrystone Inlet, cultivated clams are likely subsidized by phytoplankton from the Chesapeake Bay. Consequently, much of the N released by mineralization associated with clam cultivation is “new” N as it would not be present in the system without bivalve facilitation. Macroalgae that are fueled by the enhanced N regeneration from clams represents a eutrophying process resulting from aquaculture. This synthesis demonstrates the importance of considering impacts of bivalve aquaculture in an ecosystem context especially relative to the potential of bivalves to remove nutrients and enhance C sinks.     

Doubling of coastal erosion under rising sea level by mid-century in Hawaii

This study discusses the scaling properties of the spatial distribution of the December 26, 2004, Sumatra aftershocks. We estimate the spatial correlation dimension D ₂ of the epicentral distribution of aftershocks recorded by a local network operated by Geological Survey of India. We estimate the value of D ₂ for five blocks in the source area by using generalized correlation integral approach. We assess its bias due to finite data points, scaling range, effects of location errors, and boundary effects theoretically and apply it to real data sets. The correlation dimension was computed both for real as well as synthetic data sets that include randomly generated point sets obtained using uniform distributions and mimicking the number of events and outlines of the effective areas filled with epicenters. On comparing the results from the real data and random point sets from simulations, we found the lower limit of bias in D ₂ estimates from limited data sets to be 0.26. Thus, the spatial variation in correlation dimensions among different blocks using local data sets cannot be directly compared unless the influence of bias in the real aftershock data set is taken into account. They cannot also be used to infer the geometry of the faults. We also discuss the results in order to add constraints on the use of synthetic data and of different approaches for uncertainty analysis on spatial variation of D ₂. A difference in D ₂ values, rather than their absolute values, among small blocks is of interest to local data sets, which are correlated with their seismic b values. Taking into account the possible errors and biases, the average D ₂ values vary from 1.05 to 1.57 in the Andaman–Nicobar region. The relative change in D ₂ values can be interpreted in terms of clustering and diffuse seismic activity associated with the low and high D ₂ values, respectively. Overall, a relatively high D ₂ and low b value is consistent with high-magnitude, diffuse activity in space in the source region of the 2004 Sumatra earthquake.     

Celebrating the Career of Dr. John R. Garratt: Long-Term Proponent of Boundary-Layer Meteorology and International Man of Mystery

Boundary-Layer Meteorology -