Delineating runoff processes and critical runoff source areas in a pasture hillslope of the Ozark Highlands

The identification of runoff contributing areas would provide the ideal focal points for water quality monitoring and Best Management Practice (BMP) implementation. The objective of this study was to use a field‐scale approach to delineate critical runoff source areas and to determine the runoff mechanisms in a pasture hillslope of the Ozark Highlands in the USA. Three adjacent hillslope plots located at the Savoy Experimental Watershed, north‐west Arkansas, were bermed to isolate runoff. Each plot was equipped with paired subsurface saturation and surface runoff sensors, shallow groundwater wells, H‐flumes and rain gauges to quantify runoff mechanisms and rainfall characteristics at continuous 5‐minute intervals. The spatial extent of runoff source areas was determined by incorporating sensor data into a geographic information‐based system and performing geostatistical computations (inverse distance weighting method). Results indicate that both infiltration excess runoff and saturation excess runoff mechanisms occur to varying extents (0–58% for infiltration excess and 0–26% for saturation excess) across the plots. Rainfall events that occurred 1–5 January 2005 are used to illustrate the spatial and temporal dynamics of the critical runoff source areas. The methodology presented can serve as a framework upon which critical runoff source areas can be identified and managed for water quality protection in other watersheds. Copyright © 2008 John Wiley & Sons, Ltd.     

Hydrothermal alteration in the Reykjanes geothermal system: Insights from Iceland deep drilling program well RN-17

The Reykjanes geothermal system is a seawater-recharged hydrothermal system that appears to be analogous to seafloor hydrothermal systems in terms of host rock type and low water/rock alteration. The similarities make the Reykjanes system a useful proxy for seafloor vents. At some time during the Pleistocene, the system was dominated by meteoric water recharge, and fluid composition at Reykjanes has evolved through time as a result of changing proportions of meteoric water influx as well as differing pressure and temperature conditions. The purpose of this study is to characterize secondary mineralization, degree of metasomatic alteration, and bulk composition of cuttings from well RN-17 from the Reykjanes geothermal system. The basaltic host rock includes hyaloclastite, breccia, tuff, extrusive basalt, diabase, as well as a marine sedimentary sequence. The progressive hydrothermal alteration sequence observed with increasing depth results from reaction of geothermal fluids with the basaltic host rock. An assemblage of greenschist facies alteration minerals, including actinolite, prehnite, epidote and garnet, occurs at depths as shallow as 350 m; these minerals are commonly found in Icelandic geothermal systems at temperatures above 250 °C (Bird and Spieler, 2004). This requires hydrostatic pressures that exceed the present-day depth to boiling point curve, and therefore must record alteration at higher fluid pressures, perhaps as a result of Pleistocene glaciation. Major, minor, and trace element profiles of the cuttings indicate transitional MORB to OIB composition with limited metasomatic shifts in easily mobilized elements. Changes in MgO, K₂O and loss on ignition indicate that metasomatism is strongly correlated with protolith properties. The textures of alteration minerals reveal alteration style to be strongly dependent on protolith as well. Hyaloclastites are intensely altered with calc-silicate alteration assemblages comprising calcic hydrothermal plagioclase, grandite garnet, prehnite, epidote, hydrothermal clinopyroxene, and titanite. In contrast, crystalline basalts and intrusive rocks display a range in alteration intensity from essentially unaltered to pervasive and nearly complete albitization of igneous feldspar and uralitization of clinopyroxene. Hydrothermal anorthite (An92–An98) occurs in veins in the most altered basalt cuttings and is significantly more calcic than igneous feldspar (An48–An79). Amphibole compositions change from actinolite to hornblende at depth. Hydrothermal clinopyroxene, which occurs in veins, has greater variation in Fe content and is systematically more calcic than igneous pyroxene and also lacks uralitic textures. Solid solutions of prehnite, epidote, and garnet indicate evolving equilibria with respect to aluminum and ferric iron.     

QSOs as probes of the early universe

The spectrum of the highest redshift QSO 2000–330 (z=3.78) contains four heavy-element absorption systems withz _abs>3.0. Interesting features include velocity structure atz _abs=3.552 which suggests a cluster origin and a purely low ionization system atz _abs=3.1881 typical of a galactic disk sightline.Paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 September–6 October, 1984.     

Holocene rockfalls in the southern Negev Desert,Israel and their relation to Dead Sea fault earthquakes

Rockfall ages in tectonically active regions provide information regarding frequency and magnitude of earthquakes. In the hyper-arid environment of the Dead Sea fault (DSF), southern Israel, rockfalls are most probably triggered by earthquakes. We dated rockfalls along the western margin of the DSF using terrestrial cosmogenic nuclides (TCN). At each rockfall site, samples were collected from simultaneously exposed conjugate boulders and cliff surfaces. Such conjugate samples initially had identical pre-fall (“inherited”) TCN concentrations. After boulder detachment, these surfaces were dosed by different production rates due to differences in post-fall shielding and geometry. However, in our study area, pre-rockfall inheritance and post-rockfall production rates of TCN cannot be evaluated. Therefore, we developed a numerical approach and demonstrated a way to overcome the above-mentioned problems. This approach can be applied in other settings where rockfalls cannot be dated by simple exposure dating. Results suggest rockfall ages between 3.6 ± 0.8 and 4.7 ± 0.7 ka. OSL ages of sediment accumulated behind the boulders range between 0.6 ± 0.1 and 3.4 ± 1.4 ka and support the TCN results. Our ages agree with dated earthquakes determined in paleoseismic studies along the entire length of the DSF and support the observation of intensive earthquake activity around 4–5 ka.     

The precision of radial-velocity determinations of solar-type stars by cross-correlation

The theoretical limit to the precision (defined as the inverse of the square of the random error) of a radial-velocity determination of a solar-type star obtained from the technique of digital cross-correlation is constructed as a function of simple parameters of the recording of a stellar spectrum.The relationship is used to establish conditions for maximizing the precision for a given object, exposure time, and observing instrument. Both fibre-fed and Cassegrain-mounted spectrographs are considered.     

Analysis of a deformable fracture in permeable material

The response of deformable fractures to changes in fluid pressure controls phenomena ranging from the flow of fluids near wells to the propagation of hydraulic fractures. We developed an analysis designed to simulate fluid flows in the vicinity of asperity‐supported fractures at rest, or fully open fractures that might be propagating. Transitions between at‐rest and propagating fractures can also be simulated. This is accomplished by defining contact aperture as the aperture when asperities on a closing fracture first make contact. Locations on a fracture where the aperture is less than the contact aperture are loaded by both fluid pressure and effective stress, whereas locations where the aperture exceeds the contact aperture are loaded only by fluid pressure. Fluid pressure and effective stress on the fracture are determined as functions of time by solving equations of continuity in the fracture and matrix, and by matching the global displacements of the fracture walls to the local deformation of asperities. The resulting analysis is implemented in a numerical code that can simulate well tests or hydraulic fracturing operations. Aperture changes during hydraulic well tests can be measured in the field, and the results predicted using this analysis are similar to field observations. The hydraulic fracturing process can be simulated from the inflation of a pre‐existing crack, to the propagation of a fracture, and the closure of the fracture to rest on asperities or proppant. Two‐dimensional, multi‐phase fluid flow in the matrix is included to provide details that are obscured by simplifications of the leakoff process (Carter‐type assumptions) used in many hydraulic fracture models. Execution times are relatively short, so it is practical to implement this code with parameter estimation algorithms to facilitate interpretation of field data. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical simulations of granular dynamics II: Particle dynamics in a shaken granular material

Surfaces of planets and small bodies of our Solar System are often covered by a layer of granular material that can range from a fine regolith to a gravel-like structure of varying depths. Therefore, the dynamics of granular materials are involved in many events occurring during planetary and small-body evolution thus contributing to their geological properties.We demonstrate that the new adaptation of the parallel N-body hard-sphere code pkdgrav has the capability to model accurately the key features of the collective motion of bidisperse granular materials in a dense regime as a result of shaking. As a stringent test of the numerical code we investigate the complex collective ordering and motion of granular material by direct comparison with laboratory experiments. We demonstrate that, as experimentally observed, the scale of the collective motion increases with increasing small-particle additive concentration.We then extend our investigations to assess how self-gravity and external gravity affect collective motion. In our reduced-gravity simulations both the gravitational conditions and the frequency of the vibrations roughly match the conditions on asteroids subjected to seismic shaking, though real regolith is likely to be much more heterogeneous and less ordered than in our idealised simulations. We also show that collective motion can occur in a granular material under a wide range of inter-particle gravity conditions and in the absence of an external gravitational field. These investigations demonstrate the great interest of being able to simulate conditions that are to relevant planetary science yet unreachable by Earth-based laboratory experiments.     

Static and dynamic conceptual model of a complexly fractured crystalline rock aquifer

A conceptual model of anisotropic and dynamic permeability is developed from hydrogeologic and hydromechanical characterization of a foliated, complexly fractured, crystalline rock aquifer at Gates Pond, Berlin, Massachusetts. Methods of investigation include aquifer‐pumping tests, long‐term hydrologic monitoring, fracture characterization, downhole heat‐pulse flow meter measurements, in situ extensometer testing, and earth tide analysis. A static conceptual model is developed from observations of depth‐dependent and anisotropic permeability that effectively compartmentalizes the aquifer as a function of foliation intensity. Superimposed on the static model is dynamic permeability as a function of hydraulic head in which transient bulk aquifer transmissivity is proportional to changes in hydraulic head due to hydromechanical coupling. The dynamic permeability concept is built on observations that fracture aperture changes as a function of hydraulic head, as measured during in situ extensometer testing of individual fractures, and observed changes in bulk aquifer transmissivity as determined from earth tides during seasonal changes in hydraulic head, with higher transmissivity during periods of high hydraulic head, and lower transmissivity during periods of relatively lower hydraulic head. A final conceptual model is presented that captures both the static and dynamic properties of the aquifer. The workflow presented here demonstrates development of a conceptual framework for building numerical models of complexly fractured, foliated, crystalline rock aquifers that includes both a static model to describe the spatial distribution of permeability as a function of fracture type and foliation intensity and a dynamic model that describes how hydromechanical coupling impacts permeability magnitude as a function of hydraulic head fluctuation. This model captures important geologic controls on permeability magnitude, anisotropy, and transience and therefor offers potentially more reliable history matching and forecasts of different water management strategies, such as resource evaluation, well placement, permeability prediction, and evaluating remediation strategies.     

Investigating Bermuda's pollution history through stable isotope analyses of modern and museum-held gorgonian corals

For centuries, Bermuda has been challenged with wastewater management for the protection of human and environmental health. By quantifying the δ¹⁵N of the common sea fan Gorgonia ventalina sampled from 30 sites throughout Bermuda we show that sewage-derived nitrogen is detectable on nearshore coral reefs and declines across the lagoon to the outer rim. We also sampled gorgonians from two museum collections representing a 50y time-series (1958–2008). These samples revealed an increase in δ¹⁵N of > 4.0‰ until the mid-1970s, after which δ¹⁵N values slowly declined by ~ 2.0‰. A δ¹⁵N chronology from a gorgonian skeleton exhibited a similar decline over the last 30–40 years of approximately 0.6‰. We conclude that policies have been effective in reducing sewage impacts to Bermudian reefs. However, significant sources of sewage pollution persist and are likely have a strong impact on harbor and nearshore coral communities and human health.