Hydrobiogeochemical Interactions along “Flow Tubes” Controls: Watershed Scale Contaminant Flow and Transformation at the Oak Ridge IFRC

Geochemical,spectrographic, microbiological and hydrogeologic studies at the ORIFRC site indicate that groundwater transport in structured media may behave as a system of parallel flow tubes. These tubes are preferred flowpaths that enable contaminant transport parallel to bedding planes (strike) over distances of 1000s of meters. A significant flux of groundwater is focused within an interval defined by the interface between the competent bedrock and overlying highly-weathered saprolite, commonly referred to as the"transition zone." Characteristics of this transition zone are dense fractures and the relative absence of weathering products (e.g. clays)results in a significantly higher permeability compared to both the overlying clay-saprolite and underlying bedrock. Several stratabound low seismic velocity zones located below the transition zone were identified during geophysics studies and were also determined to be fractured high permeability preferred contaminant transport pathways during subsequent drilling activities. XANES analysis of precipitates collected from these deeper flow zones indicate 95% or more of the U deposited is U(VI). Linear combination fitting of the EXAFS data shows that precipitates are ~51±5% U(VI)-carbonate-like phase (e.g., liebigite) and ~49±5% U(VI) associated with an iron oxide phase; inclusion of a third component in the fit suggests that up to 15% of the U(VI) may be associated with a phosphate phase or OH- phase (e.g.,schoepite). Although precipitates with similar U(VI)-carbonate and/or phosphate associations were identified in the transition zone pathways,there were also U(VI) complexes adsorbed to mineral surfaces that would tend to be more readily mobilized. Groundwater in the different flow tubes has been determined to consist of different water quality types that vary with the solid phase encountered (e.g., clays, carbonates, clastics) as contaminants migrate along the flow paths. This lateral and vertical variability in geochemistry, particularly pH, has a significant impact on microbiological community composition and activity. Ribosomal RNA gene analyses coupled with physiological and genomic analyses suggest that bacteria from the genus Rhodanobacter(a diverse population of denitrifiers that are moderately acid tolerant) have a high relative abundance in the acidic source zone at the ORIFRC site.Watershed-scale analysis across different flow paths/tubes revealed strong negative correlation between pH and the absolute and relative abundance of Rhodanobacter. Recent studies also confirmed that the ORIFRC site hosts a diverse fungal community, with significant differences observed between acidic (pH <5) and circumneutral (>5) wells. The lack of nitrous oxide reduction capability in fungi, and the detection of denitrification potential in slurry microcosms suggest that fungi may have aheretofore under appreciated role in biogeochemical transformations, with implications forsite remediation and greenhouse gas emissions. Further research is needed to determine if these organisms can influence U(VI) mobility either directly through immobilization or indirectly through the depletion of nitrate.In conclusion, additional studies are required to quantify the processes (e.g., solid phase reactions, recharge, diffusion, microbial interactions) that are occurring along the groundwater flow tubes identified at the ORIFRC so predictive models can be parameterized and used to assess long-term contaminant fate and transport and remedial options.

Anisotropic seismic inversion using a multigrid Monte Carlo approach

We propose a new approach for the inversion of anisotropic P-wave data based on Monte Carlo methods combined with a multigrid approach. Simulated annealing facilitates objective minimization of the functional characterizing the misfit between observed and predicted traveltimes, as controlled by the Thomsen anisotropy parameters (epsilon, delta). Cycling between finer and coarser grids enhances the computational efficiency of the inversion process, thus accelerating the convergence of the solution while acting as a regularization technique of the inverse problem. Multigrid perturbation samples the probability density function without the requirements for the user to adjust tuning parameters. This increases the probability that the preferred global, rather than a poor local, minimum is attained. Undertaking multigrid refinement and Monte Carlo search in parallel produces more robust convergence than does the initially more intuitive approach of completing them sequentially. We demonstrate the usefulness of the new multigrid Monte Carlo (MGMC) scheme by applying it to (a) synthetic, noise-contaminated data reflecting an isotropic subsurface of constant slowness, horizontally layered geologic media and discrete subsurface anomalies; and (b) a crosshole seismic data set acquired by previous authors at the Reskajeage test site in Cornwall, UK. Inverted distributions of slowness (s) and the Thomson anisotropy parameters (epsilon, delta) compare favourably with those obtained previously using a popular matrix-based method. Reconstruction of the Thomsen epsilon parameter is particularly robust compared to that of slowness and the Thomsen delta parameter, even in the face of complex subsurface anomalies. The Thomsen epsilon and delta parameters have enhanced sensitivities to bulk-fabric and fracture-based anisotropies in the TI medium at Reskajeage. Because reconstruction of slowness (s) is intimately linked to that epsilon and delta in the MGMC scheme, inverted images of phase velocity reflect the integrated effects of these two modes of anisotropy. The new MGMC technique thus promises to facilitate rapid inversion of crosshole P-wave data for seismic slownesses and the Thomsen anisotropy parameters, with minimal user input in the inversion process.

Detection of soil compaction using seismic surface waves

Seismic geophysical methods have rarely been used in precision agriculture, predominantly due to the perception that they are slow and results require a complex evaluation. This paper explores the possibility of using a recently developed surface wave seismic geophysical approach, the multichannel analysis of surface waves (MASW) method, for assessment of agricultural compaction. This approach has the advantage of being non-intrusive, rapid and is able to produce 2D ground models with a relatively high density of spatial sampling points. The method, which was tested on a research site in Oakpark, Ireland, detected a significant difference in shear wave velocity between a heavily compacted headland and an uncompacted location. The results from this approach compared favourably with those obtained
from measurements of bulk density and penetrometer resistance and demonstrate that the MASW approach can distinguish between the extreme states of heavily compacted and uncompacted soil.

Vibration based seismic damage identification procedure for a steel concentrically braced frame structure

An approach for seismic damage identification of a single-storey steel concentrically braced frame (CBF) structure is presented through filtering and double integration of a recorded acceleration signal. A band-pass filter removes noise from the acceleration signal followed by baseline correction being used to reduce the drift in velocity and displacement during numerical integration. The pre-processing achieves reliable numerical integration that predicts the displacement response accurately when compared to the measured lateral in-plane displacement of the CBF structure. The lateral displacement of the CBF structure is used to infer buckling and yielding of bracing members through seismic tests. The level of interstorey drift of the CBF during a seismic excitation allows the yield and buckling of the bracing members to be identified and indirectly detects damage based on exceedance of calculated displacement limits. The calculated buckling and yielding displacement threshold limits used to identify damage are demonstrated to accurately identify initial buckling and yielding in the bracing members.

Seismic damage detection for a steel braced frame structure

A new approach for global detection of seismic damage in a single-storey steel concentrically braced frame (CBF) structure is presented. The filtered lateral in-plane acceleration response of the CBF structure is integrated twice to provide the lateral in-plane displacement which is used to infer buckling and yielding damage. The level of interstorey drift of the CBF during a seismic excitation allows the yield and buckling of the bracing members to be identified and indirectly detects damage based on exceedance of calculated lateral in-plane displacement limits. A band-pass filter removes noise from the acceleration signal followed by baseline correction being used to reduce the drift in velocity and displacement during numerical integration. This pre-processing results in reliable numerical integration of the frame acceleration that predicts the displacement response accurately when compared to the measured lateral displacement of the CBF structure. Importantly, the structural damage is not assumed through removal of bracing members, rather damage is induced through actual seismic loading. The buckling and yielding displacement threshold limits used to identify damage are demonstrated to accurately identify the initiation of buckling and yielding.

Difference equation approach to two-thermocouple sensor characterisation in constant velocity flow environments

Thermocouples are one of the most popular devices for temperature measurement due to their robustness, ease of manufacture and installation, and low cost. However, when used in certain harsh environments, for example, in combustion systems and engine exhausts, large wire diameters are required, and consequently the measurement bandwidth is reduced. This article discusses a software compensation technique to address the loss of high frequency fluctuations based on measurements from two thermocouples. In particular, a difference equation sDEd approach is proposed and compared with existing methods both in simulation and on experimental test rig data with constant flow velocity. It is found that the DE algorithm, combined with the use of generalized total least squares for parameter identification, provides better performance in terms of time constant estimation without any a priori assumption on the time constant ratios of the thermocouples.

Effect of a COL1A1 Sp1 Binding Site Polymorphism on Arterial Pulse Wave Velocity: An Index of Compliance.

Reduced arterial compliance precedes changes in blood pressure, which may be mediated through alterations in vessel wall matrix composition. We investigated the effect of the collagen type I-1 gene (COL1A1) +2046G>T polymorphism on arterial compliance in healthy individuals. We recruited 489 subjects (251 men and 238 women; mean age, 22.6±1.6 years). COL1A1 genotypes were determined using polymerase chain reaction and digestion by restriction enzyme Bal1. Arterial pulse wave velocities were measured in 3 segments, aortoiliac (PWVA), aortoradial (PWVB), and aorto-dorsalis-pedis (PWVF), as an index of compliance using a noninvasive optical method. Data were available for 455 subjects. The sample was in Hardy-Weinberg equilibrium with genotype distributions and allele frequencies that were not significantly different from those reported previously. The T allele frequency was 0.22 (95% confidence interval, 0.19 to 0.24). Two hundred eighty-three (62.2%) subjects were genotype GG, 148 (35.5%) subjects were genotype GT, and 24 (5.3%) subjects were genotype TT. A comparison of GG homozygotes with GT and TT individuals demonstrated a statistically significant association with arterial compliance: PWVF 4.92±0.03 versus 5.06±0.05 m/s (ANOVA, P=0.009), PWVB 4.20±0.03 versus 4.32±0.04 m/s (ANOVA, P=0.036), and PWVA 3.07±0.03 versus 3.15±0.03 m/s (ANOVA, P=0.045). The effects of genotype were independent of age, gender, smoking, mean arterial pressure, body mass index, family history of hypertension, and activity scores. We report an association between the COL1A1 gene polymorphism and arterial compliance. Alterations in arterial collagen type 1A deposition may play a role in the regulation of arterial compliance

CaII K interstellar observations towards early-type disc and halo stars - distances to intermediate- and high-velocity clouds

We compare existing high spectral resolution (R = lambda/Deltalambda similar to 40 000) Ca II Kobservations (lambda(air) = 3933.66 Angstrom) towards 88 mainly B-type stars, and new observations taken using the Intermediate dispersion Spectrograph and Imaging System (ISIS) on the William Herschel Telescope at R similar to 10 000 towards three stars taken from the Palomar-Green Survey, with 21-cm HI emission-line profiles, in order to search for optical absorption towards known intermediate- and high-velocity cloud complexes. Given certain assumptions, limits to the gas phase abundance of Ca II are estimated for the cloud components. We use the data to derive the following distances from the Galactic plane (z). (i) Tentative lower z-height limits of 2800 and 4100 pc towards complex C using lack of absorption in the spectra of HD341617 and PG 0855 + 294, respectively. (ii) A weak lower z-height of 1400 pc towards complex WA-WB using lack of absorption in EC 09470-1433 and a weak lower limit of 2470 pc using lack of absorption in EC 09452-1403. (iii) An upper z- height of 2470 pc towards a southern intermediate- velocity cloud (IVC) with v(LSR) = -55 km s(-1) using PG 2351 + 198. (iv) Detection of a possible IVC in Ca II absorption at v(LSR) = +52 km s(-1) using EC 20104-2944. No associated HI in emission is detected. At this position, normal Galactic rotation predicts velocities of up to similar to+ 25 km s(-1). The detection puts an upper z-height of 1860 pc to the cloud. (v) Tentative HI and Ca II K detections towards an IVC at similar to+70 km s(-1) in the direction of high-velocity cloud (HVC) complex WE, sightline EC 06387-8045, indicating that the IVC may be at a z-height lower than 1770 pc. (vi) Detection of Ca II K absorption in the spectrum of PG 0855 + 294 in the direction of IV20, indicating that this IVC has a z-height smaller than 4100 pc. (vii) A weak lower z-height of 4300 pc towards a small HVC with v(LSR) = +115 km s(-1) at l, b = 200degrees, + 52degrees, using lack of absorption in the Ca II K spectrum of PG 0955 + 291.

HST/STIS observations of the high-velocity interstellar cloud HVC 291.2-41.2+80: a warm, mainly ionized high-velocity cloud

We present intermediate-resolution HST/STIS spectra of a high- velocity interstellar cloud ((LSR)-L-upsilon = + 80 kms(-1)) towards DI1388, a young star in the Magellanic Bridge located between the Small and Large Magellanic Clouds. The STIS data have a signal-to-noise ratio (S/N) of 20-45 and a spectral resolution of about 6.5 km s(-1) (FWHM), The high-velocity cloud absorption is observed in the lines of C II, O I, Si II, Si III, Si IV and S III. Limits can be placed on the amount of S II and Fe II absorption that is present. An analysis of the relative abundances derived from the observed species, particularly C II and O I, suggests that this high-velocity gas is warm (T-k similar to 10(3)-10(4) K) and predominantly ionized, This hypothesis is supported by the presence of absorption produced by highly ionized species, such as Si IV, This sightline also intercepts two other high-velocity clouds that produce weak absorption features at (LSR)-L-upsilon = + 113 and + 130kms(-1) in the STIS spectra.

Multiwavelength observations of the M15 intermediate-velocity cloud

We present Westerbork Synthesis Radio Telescope HI images, Lovell telescope multibeam H I wide-field mapping, William Herschel Telescope long-slit echelle Ca II observations, Wisconsin Halpha Mapper (WHAM) facility images, and IRAS ISSA 60- and 100-mum co-added images towards the intermediate- velocity cloud (IVC) at + 70 km s(-1), located in the general direction of the M15 globular cluster. When combined with previously published Arecibo data, the H I gas in the IVC is found to be clumpy, with a peak H I column density of similar to1.5 x 10(20) cm(-2), inferred volume density (assuming spherical symmetry) of similar to24 cm(-3)/D (kpc) and a maximum brightness temperature at a resolution of 81 x 14 arcsec(2) of 14 K. The major axis of this part of the IVC lies approximately parallel to the Galactic plane, as does the low- velocity H I gas and IRAS emission. The H I gas in the cloud is warm, with a minimum value of the full width at half-maximum velocity width of 5 km s(-1) corresponding to a kinetic temperature, in the absence of turbulence, of similar to540 K. From the H I data, there are indications of two-component velocity structure. Similarly, the Ca II spectra, of resolution 7 km s(-1), also show tentative evidence of velocity structure, perhaps indicative of cloudlets. Assuming that there are no unresolved narrow-velocity components, the mean values of log(10)[N(Ca II K) cm(2)] similar to 12.0 and Ca II/H I similar to2 5 x 10(-8) are typical of observations of high Galactic latitude clouds. This compares with a value of Ca II/H I>10(-6) for IVC absorption towards HD 203664, a halo star of distance 3 kpc, some 3.degrees1 from the main M15 IVC condensation. The main IVC condensation is detected by WHAM in Halpha with central local-standard-of-rest velocities of similar to60-70 km s(-1), and intensities uncorrected for Galactic extinction of up to 1.3 R, indicating that the gas is partially ionized. The FWHM values of the Halpha IVC component, at a resolution of 1degrees, exceed 30 km s(-1). This is some 10 km s(-1) larger than the corresponding H I value at a similar resolution, and indicates that the two components may not be mixed. However, the spatial and velocity coincidence of the Halpha and H I peaks in emission towards the main IVC component is qualitatively good. If the Halpha emission is caused solely by photoionization, the Lyman continuum flux towards the main IVC condensation is similar to2.7 x 10(6) photon cm(-2) s(-1). There is not a corresponding IVC Halpha detection towards the halo star HD 203664 at velocities exceeding similar to60 km s(- 1). Finally, both the 60- and 100-mum IRAS images show spatial coincidence, over a 0.675 x 0 625 deg(2) field, with both low- and intermediate-velocity H I gas (previously observed with the Arecibo telescope), indicating that the IVC may contain dust. Both the Halpha and tentative IRAS detections discriminate this IVC from high-velocity clouds, although the H I properties do not. When combined with the H I and optical results, these data point to a Galactic origin for at least parts of this IVC.