Mapping of road‐salt‐contaminated groundwater discharge and estimation of chloride load to a small stream in southern New Hampshire,USA

Concentrations of chloride in excess of State of New Hampshire water‐quality standards (230 mg/l) have been measured in watersheds adjacent to an interstate highway (I‐93) in southern New Hampshire. A proposed widening plan for I‐93 has raised concerns over further increases in chloride. As part of this effort, road‐salt‐contaminated groundwater discharge was mapped with terrain electrical conductivity (EC) electromagnetic (EM) methods in the fall of 2006 to identify potential sources of chloride during base‐flow conditions to a small stream, Policy Brook. Three different EM meters were used to measure different depths below the streambed (ranging from 0 to 3 m). Results from the three meters showed similar patterns and identified several reaches where high EC groundwater may have been discharging. Based on the delineation of high (up to 350 mmhos/m) apparent terrain EC, seven‐streambed piezometers were installed to sample shallow groundwater. Locations with high specific conductance in shallow groundwater (up to 2630 mmhos/m) generally matched locations with high streambed (shallow subsurface) terrain EC. A regression equation was used to convert the terrain EC of the streambed to an equivalent chloride concentration in shallow groundwater unique for this site. Utilizing the regression equation and estimates of one‐dimensional Darcian flow through the streambed, a maximum potential groundwater chloride load was estimated at 188 Mg of chloride per year. Changes in chloride concentration in stream water during streamflow recessions showed a linear response that indicates the dominant process affecting chloride is advective flow of chloride‐enriched groundwater discharge. Published in 2010 by John Wiley & Sons, Ltd.     

Directional chemical variations in diamonds showing octahedral following cuboid growth

A progression from cuboid to octahedral growth has been observed in 16 natural diamonds from Yakutian kimberlites. X-ray and cathodoluminescence topography have revealed that the change in morphology of diamonds with cloudy cuboid cores may occur without mixed-habit growth but via generation of numerous octahedral apices on cuboid surfaces and subsequent gradual transformation into regular octahedral morphology. Nitrogen aggregation in both cuboid and octahedral domains of such diamonds suggests that they have had a long residence time under mantle conditions. Micro-inclusions in the cuboid domains of the diamonds testify to the nucleation and growth of cuboid cores from a hydrous-carbonatitic (oxidized) fluid. The transition from cuboid hummocky growth rich in inclusions to octahedral growth without inclusions may be linked to decreasing supersaturation in the parent fluid. Measurements of δ¹³C and N_ppm by ion microprobe show that the chemical variations observed between inner cuboid domains and outer octahedral zones commonly have a systematic character and as such they are probably not due to purely kinetic effects. The peripheral octahedral zones are always enriched in ¹³C in comparison with inner cuboid ones, and the total nitrogen content decreases with the change from cuboid to octahedral growth. The octahedral outer zones show a gradual progressive increase in δ¹³C, with an overall change of up to 5‰ from the cuboid core (δ¹³C usually between −8 and −6‰) to the diamond margin (δ¹³C usually between −4 and −2‰). Decreases in δ¹³C of this magnitude with a gradual increase in ¹³C may be attributed to the Rayleigh fractionation operating on a single parent fluid of close to normal mantle δ¹³C composition with diamond precipitating by the reduction of carbonatitic fluid in a closed system. However, one sample shows a variation of δ¹³C of approximately −17 to −6‰ and therefore suggests a possible change of fluid source composition from one containing subducted crustal organic carbon to one with common mantle carbon. An erratum to this article can be found at     

The impacts of the 2004 hurricanes on hydrology, water quality, and seagrass in the Central Indian River Lagoon, Florida

Between August 14 and September 26, 2004, four tropical weather systems (Charley, Frances, Ivan, and Jeanne) affected the central Indian River Lagoon (IRL). The central IRL received a prodigious amount of rainfall for the 2 mo, between 72 and 83 cm, which is a once-in-50-yr rainfall event. High stream discharges were generated that, combined with wind-suspended sediments, significantly reduced salinities and water transparency. In September, salinities among central IRL segments dropped from 30 psu or more to ≤15 psu, color increased from a low of 10 pcu to ≥100 pcu, and turbidity increased from ≤3 NTU up to 14 NTU. Evidence of the hurricanes' physical effects on seagrasses (burial, no scour) was limited to just one of the more than 25 sites inspected. Within 2 to 3 mo following the hurricane period, most parameters related to water transparency returned to or showed improvement over their prehurricane (February–July 2004) levels. Unseasonably low salinities (<20 psu) and moderately high color (>20 pcu) were observed through spring 2005, largely attributable to a relatively long residence time and a wetter-than-average spring season in 2005. By the end of the study period (July 2006), the central IRL generally showed a continuation of two opposite seagrass trends—an increase in depthlimit coverage but a decline in coverage density—that began before 2004. Also, within a limited reach of the central IRL, there was a temporary shift in species composition in summer 2005 (Ruppia maritima increased asHalodule wrightü decreased). It is likely that the persistently low salinities (not color) in 2004–2005 affected the species composition and coverage density. This study reveals that seagrasses are resilient to the acute effects of hurricanes and underscores the need to reduce chronic, an thropogenic effects on seagrasses.     

Comparison of an algebraic multigrid algorithm to two iterative solvers used for modeling ground water flow and transport

Numerical solution of large-scale ground water flow and transport problems is often constrained by the convergence behavior of the iterative solvers used to solve the resulting systems of equations. We demonstrate the ability of an algebraic multigrid algorithm (AMG) to efficiently solve the large, sparse systems of equations that result from computational models of ground water flow and transport in large and complex domains. Unlike geometric multigrid methods, this algorithm is applicable to problems in complex flow geometries, such as those encountered in pore-scale modeling of two-phase flow and transport. We integrated AMG into MODFLOW 2000 to compare two- and three-dimensional flow simulations using AMG to simulations using PCG2, a preconditioned conjugate gradient solver that uses the modified incomplete Cholesky preconditioner and is included with MODFLOW 2000. CPU times required for convergence with AMG were up to 140 times faster than those for PCG2. The cost of this increased speed was up to a nine-fold increase in required random access memory (RAM) for the three-dimensional problems and up to a four-fold increase in required RAM for the two-dimensional problems. We also compared two-dimensional numerical simulations of steady-state transport using AMG and the generalized minimum residual method with an incomplete LU-decomposition preconditioner. For these transport simulations, AMG yielded increased speeds of up to 17 times with only a 20% increase in required RAM. The ability of AMG to solve flow and transport problems in large, complex flow systems and its ready availability make it an ideal solver for use in both field-scale and pore-scale modeling.     

Surface coatings on quartz grains in bentonites and their relevance to human health

The cytotoxicity of quartz in the human lung is recognized to be dependent on both the inherent properties of the silica dust and external factors related to the history of the dust and including the presence of surface contamination. In this study, the physical and chemical surface properties of quartz grains in commercial bentonite deposits from the western (South Dakota) and southern (Alabama) USA were investigated. Measured quartz contents of bentonites range from 1.9 to 8.5 wt% with the <10 μm size fraction comprising 6–45% of this total. Trace element contents (Fe–Ti–Al) of quartz grains from any given bentonite are similar, indicating a single origin for the quartz with little if any contamination from other sources. Surface coatings are pervasive on all quartz grains and resist removal by repeated vigorous washings and reaction with HCl. Textural attributes and XPS and EDS analyses of these coatings are consistent with most being montmorillonite and, less frequently, mixtures of montmorillonite and opaline silica. Opaline silica (opal-A and opal-CT) occurs in two texturally distinct generations: an early massive grain-coating event and as later lepispheres. Montmorillonite coating thicknesses range from <1 μm to more than 10 μm thick. Surfaces of plagioclase, K-feldspar, and biotite grains are conspicuously devoid of montmorillonite coatings, but may show sparse distributions of opal-CT lepispheres. HRTEM has not confirmed a topotactic relationship or atomic structural concordance between montmorillonite coatings and underlying quartz grains. Alternatively, a precursor volcanic glass phase that coats the quartz surfaces during volcanic eruption and/or preferential early precipitation of opaline silica on quartz may provide substrates for development of montmorillonite coatings. Estimations of montmorillonite biodurability under pulmonary pH conditions suggest possible prolonged sequestration of respired bentonite quartz grains from contact with lung materials and modified cytotoxic reactivity.     

A comparison of photolysis rate parameters estimated from measured and simulated actinic flux for wintertime conditions at Storm Peak Laboratory,Colorado

Photolysis rate parameters depend upon solar actinic flux and chemical species dependent quantum yields and cross sections. Spectrally resolved measurements of actinic flux should be preferred over flux derived from models for the analysis of field observations. Actinic flux can be difficult to derive from the irradiance measurements of flat-plate radiometers. It is also difficult to estimate from models due to uncertainties in the ozone column, aerosol concentrations and distributions, cloud cover, optical depth and surface albedo. A series of actinic flux measurements were performed at Storm Peak Laboratory (3,210 m above sea level), Colorado, United States with spectroradiometers during the wintertime (January 07–10, 2004). The site is relatively remote with a clean atmosphere and during the wintertime the ground is generally covered by fresh snow with a high albedo. The actinic flux measurements were used to estimate the photolysis rate parameters of ozone, nitrogen dioxide and formaldehyde. The measured actinic flux and the photolysis rate parameters derived from the flux were compared to calculations using the Tropospheric Ultraviolet-Visible Model (TUV), version 4.2 (Madronich and Flocke, 1998). The TUV modeled actinic flux, the measured flux and the photolysis rate parameters derived from them had similar temporal patterns. However there were significant differences in their magnitude due to uncertainties in the data available to initialize the TUV model and the calibration of the spectroradiometer.     

The Entropy Score and its Uses in Earthquake Forecasting

Suppose a forecasting scheme associates a probability p^* with some observed outcome. The entropy score given to this forecast is then –logp^*. This article provides a review of the background to this scoring method, its main properties, and its relationships to concepts such as likelihood, probability gain, and Molchans - diagram. It is shown that, in terms of this score, an intrinsic characterization can be given for the predictability of a given statistical forecasting model. Uses of the score are illustrated by applications to the stress release and ETAS models, electrical signals, and M8.     

Origin of silicic volcanic rocks in Central Costa Rica: a study of a chemically variable ash-flow sheet in the Tiribí Tuff

Chemical heterogeneities of pumice clasts in an ash-flow sheet can be used to determine processes that occur in the magma chamber because they represent samples of magma that were erupted at the same time. The dominant ash-flow sheet in the Tiribí Tuff contains pumice clasts that range in composition from 55.1 to 69.2 wt% SiO₂. It covers about 820 km² and has a volume of about 25 km³ dense-rock equivalent (DRE). Based on pumice clast compositions, the sheet can be divided into three distinct chemical groupings: a low-silica group (55.1-65.6 wt% SiO₂), a silicic group (66.2-69.2 wt% SiO₂), and a mingled group (58.6-67.7 wt% SiO₂; all compositions calculated 100% anhydrous). Major and trace element modeling indicates that the low-silica magma represents a mantle melt that has undergone fractional crystallization, creating a continuous range of silica content from 55.1-65.6 wt% SiO₂. Eu/Eu*, MREE, and HREE differences between the two groups are not consistent with crystal fractionation of the low-silica magma to produce the silicic magma. The low-silica group and the silicic group represent two distinct magmas, which did not evolve in the same magma chamber. We suggest that the silicic melts resulted from partial melting of relatively hot, evolved calc-alkaline rocks that were previously emplaced and ponded at the base of an over-thickened basaltic crust. The mingled group represents mingling of the two magmas shortly before eruption. Electronic supplementary material to this paper can be obtained by using the Springer LINK server located at http://dx.doi.org/10.1007/s00445-001-0188-8.     

Arctic aerosol and Arctic climate: Results from an energy budget model

An energy budget model is used to study the effect on Arctic climate of optically active aerosol in the Arctic atmosphere. The dependence of the change in surface temperature on the vertical distribution of the aerosol and on the radiative properties of the aerosol-free atmosphere, the Arctic surface, and the aerosol, itself, are calculated. An extensive sensitivity analysis is performed to assess the degree to which the results of the model are dependent upon the assumptions underlying it.List of Symbols Used I ₀ Solar flux at the top of the Arctic Atmosphere (Arctic here means 70° N latitude to the pole) - a _S Surface albedo of the Arctic (a _S ^c is the value of surface albedo at which the sign of the surface temperature perturbation changes) - Reflection coefficient of the aerosol-free Arctic atmosphere - Absorption coefficient of the aerosol-free Arctic atmosphere - Transmission coefficient of the aerosol-free Arctic atmosphere - RI ₀ Total flux of sunlight reflected from the Arctic - A _A I ₀ Total flux of sunlight absorbed in the Arctic atmosphere - A _S I ₀ Total flux of sunlight absorbed at the Arctic surface - A _aer I ₀ Total flux of sunlight absorbed in the Arctic aerosol - Q _A Net atmospheric flow of energy, per unit of Arctic surface area, north across 70° N latitude - Q _S Net oceanic flow of energy, per unit of Arctic surface area, north across 70° N latitude - E Convective plus latent heat fluxes from surface to atmosphere - F _A Net flow of energy to the Arctic atmosphere - F _S Net flow of energy to the Arctic surface - T _A An effective temperature of the Arctic atmosphere - T _S Surface temperature of the Arctic - w Single-scattering albedo of the aerosol - t Optical depth of the aerosol - g Fraction of incident radiation scattered forward by the aerosol - Reflection coefficient of the aerosol - Absorption coefficient of the aerosol - Transmission coefficient of the aerosol - p,q Number of atmospheric layers and the inverse of the fraction of incident IR absorbed in each layer in the energy budget model - F,G,H Measures of the amount of IR-active atmosphere above the surface, the aerosol, and the clouds