Continental-scale patterns in soil geochemistry and mineralogy: Results from two transects across the United States and Canada

In 2004, the US Geological Survey (USGS) and the Geological Survey of Canada (GSC) initiated a pilot study that involved collection of more than 1500 soil samples from 221 sites along two continental transects across Canada and the United States. The pilot study was designed to test and refine protocols for a soil geochemical survey of North America. The two transects crossed a wide array of soil parent materials, soil ages, climatic conditions, landforms, land covers and land uses. Sample sites were selected randomly at approximately 40-km intervals from a population defined as all soils of the continent. At each site, soils representing 0 to 5 cm depth, and the O, A, and C horizons, if present, were collected and analyzed for their near-total content of over 40 major and trace elements. Soils from 0–5 cm depth were also collected for analysis of organic compounds. Results from the transects confirm that soil samples collected at a 40-km spacing reveal coherent, continental- to subcontinental-scale geochemical and mineralogical patterns that can be correlated to aspects of underlying soil parent material, soil age and climate influence. The geochemical data also demonstrate that at the continental-scale the dominance of any of these major factors that control soil geochemistry can change across the landscape. Along both transects, soil mineralogy and geochemistry change abruptly with changes in soil parent materials. However, the chemical influence of a soil’s parent material can be obscured by changing climatic conditions. For the transects, increasing precipitation from west to east and increasing temperature from north to south affect both soil mineralogy and geochemistry because of climate effects on soil weathering and leaching, and plant productivity. Regional anomalous metal concentrations can be linked to natural variations in soil parent materials, such as high Ni and Cr in soils developed on ultramafic rocks in California or high P in soils formed on weathered Ordovician limestones in central Kentucky. On local scales, anomalous metal concentrations recognized in soil profiles, such as high P in soils from animal confinement sites, are consistent with local anthropogenic disturbances. At a larger scale, the distribution of Hg across the west to east transect demonstrates that it can be difficult to distinguish between natural or anthropogenic contributions and that many factors can contribute to an element’s spatial distribution.Only three samples in a subset of seventy-three 0–5 cm depth soil samples from the north to south transect had organochlorine pesticides values above the method detection limit, apparently related to historic usage of the pesticides DDT and dieldrin.     

Isotopic constraints on ice age fluids in active geothermal systems: Reykjanes, Iceland

The Reykjanes geothermal system is located on the landward extension of the Mid-Atlantic Ridge in southwest Iceland, and provides an on-land proxy to high-temperature hydrothermal systems of oceanic spreading centers. Previous studies of elemental composition and salinity have shown that Reykjanes geothermal fluids are likely hydrothermally modified seawater. However, δD values of these fluids are as low as −23‰, which is indicative of a meteoric water component. Here we constrain the origin of Reykjanes hydrothermal solutions by analysis of hydrogen and oxygen isotope compositions of hydrothermal epidote from geothermal drillholes at depths between 1 and 3 km. δD_EPIDOTE values from wells RN-8, -9, -10 and -17 collectively range from −60 to −78‰, and δ¹⁸O_EPIDOTE in these wells are between −3.0 and 2.3‰. The δD values of epidote generally increase along a NE trend through the geothermal field, whereas δ¹⁸O values generally decrease, suggesting a southwest to northeast migration of the geothermal upflow zone with time that is consistent with present-day temperatures and observed hydrothermal mineral zones. For comparative analysis, the meteoric-water dominated Nesjavellir and Krafla geothermal systems, which have a δD_FLUID of ∼ −79‰ and −89‰, respectively, show δD_EPIDOTE values of −115‰ and −125‰. In contrast, δD_EPIDOTE from the mixed meteoric-seawater Svartsengi geothermal system is −68‰; comparable to δD_EPIDOTE from well RN-10 at Reykjanes.Stable isotope compositions of geothermal fluids in isotopic equilibrium with the epidotes at Reykjanes are computed using published temperature dependent hydrogen and oxygen isotope fractionation curves for epidote-water, measured isotope composition of the epidotes and temperatures approximated from the boiling point curve with depth. Calculated δD and δ¹⁸O of geothermal fluids are less than 0‰, suggesting that fluids of meteoric or glacial origin are a significant component of the geothermal solutions. Additionally, δD_FLUID values in equilibrium with geothermal epidote are lower than those of modern-day fluids, whereas calculated δ¹⁸O_FLUID values are within range of the observed fluid isotope composition. We propose that modern δD_EPIDOTE and δD_FLUID values are the result of diffusional exchange between hydrous alteration minerals that precipitated from glacially-derived fluids early in the evolution of the Reykjanes system and modern seawater-derived geothermal fluids. A simplified model of isotope exchange in the Reykjanes geothermal system, in which the average starting δD_ROCK value is −125‰ and the water to rock mass ratio is 0.25, predicts a δD_FLUID composition within 1‰ of average measured values. This model resolves the discrepancy between fluid salinity and isotope composition of Reykjanes geothermal fluids, explains the observed disequilibrium between modern fluids and hydrothermal epidote, and suggests that rock-fluid interaction is the dominant control over the evolution of fluid isotope composition in the hydrothermal system.     

Multiscale ensemble filtering for reservoir engineering applications

Reservoir management requires periodic updates of the simulation models using the production data available over time. Traditionally, validation of reservoir models with production data is done using a history matching process. Uncertainties in the data, as well as in the model, lead to a nonunique history matching inverse problem. It has been shown that the ensemble Kalman filter (EnKF) is an adequate method for predicting the dynamics of the reservoir. The EnKF is a sequential Monte-Carlo approach that uses an ensemble of reservoir models. For realistic, large-scale applications, the ensemble size needs to be kept small due to computational inefficiency. Consequently, the error space is not well covered (poor cross-correlation matrix approximations) and the updated parameter field becomes scattered and loses important geological features (for example, the contact between high- and low-permeability values). The prior geological knowledge present in the initial time is not found anymore in the final updated parameter. We propose a new approach to overcome some of the EnKF limitations. This paper shows the specifications and results of the ensemble multiscale filter (EnMSF) for automatic history matching. EnMSF replaces, at each update time, the prior sample covariance with a multiscale tree. The global dependence is preserved via the parent–child relation in the tree (nodes at the adjacent scales). After constructing the tree, the Kalman update is performed. The properties of the EnMSF are presented here with a 2D, two-phase (oil and water) small twin experiment, and the results are compared to the EnKF. The advantages of using EnMSF are localization in space and scale, adaptability to prior information, and efficiency in case many measurements are available. These advantages make the EnMSF a practical tool for many data assimilation problems.     

Responses of <Emphasis Type="Italic">Avicennia germinans</Emphasis> (Black Mangrove) and the Soil Microbial Community to Nitrogen Addition in a Hypersaline Wetland

The responses of dwarf black mangroves (Avicennia germinans) and components of the microbial community to the addition of nitrate over a 2-year period were examined. The field study was conducted in the Indian River Lagoon (Florida) in a mangrove-dominated impoundment that was established for purposes of mosquito control. The responses of mangroves to the regular addition of nitrate were insignificant or relatively minor compared to responses of the components of the microbial community. Denitrification rates, measured in the field and laboratory, increased significantly in fertilized plots and nitrous oxide emission rates were almost six times higher in fertilized plots. Nitrogen fixation was significantly lower in fertilized plots. Results suggest that mangrove systems in the N-limited Indian River Lagoon are likely to be long-term sinks for any increases in N loading.     

Landscape structure and use,climate, and water movement in the Mekong River basin

In this article the relative roles of precipitation and soil moisture in influencing runoff variability in the Mekong River basin are addressed. The factors controlling runoff generation are analysed in a calibrated macro‐scale hydrologic model, and it is demonstrated that, in addition to rainfall, simulated soil moisture plays a decisive role in establishing the timing and amount of generated runoff. Soil moisture is a variable with a long memory for antecedent hydrologic fluxes that is influenced by soil hydrologic parameters, topography, and land cover type. The influence of land cover on soil moisture implies significant hydrologic consequences for large‐scale deforestation and expansion of agricultural land. Copyright © 2007 John Wiley & Sons, Ltd.     

Slab dehydration and basalt petrogenesis in subduction systems involving very young oceanic lithosphere

Compositions of post-Miocene basalts erupted in the Garibaldi and Central America volcanic arcs exhibit significant correlations with the age of the subducted plate. In general, SiO₂, Al₂O₃, CaO, V, and (Sr/P)_N decrease and FeO, MgO, TiO₂ and Na₂O increase as the age of the subducted plate decreases. Variations in CaO/Al₂O₃, SiO₂, (Sr/P)_N, and Ba are compatible with lesser slab input, and hence less hydrous melting conditions in the mantle wedge in segments of the arcs overlying the youngest oceanic lithosphere. This interpretation is supported by comparison with peridotite melting experiments, which suggest higher melt pressures and temperatures in the mantle wedge above very young oceanic lithosphere. These observations point to a model in which dehydration of the downgoing slab occurs at shallow depths in subduction systems involving oceanic lithosphere younger than about 20 Ma. Because young oceanic lithosphere is relatively warm, little post-subduction heating is required to produce metamorphic reactions that release slab volatiles. Geodynamic models indicate most volatile-liberating reactions will occur within the seismogenic zone in oceanic lithosphere younger than 20 Ma, thus limiting the volatile flux beneath the arc and encouraging drier, higher temperature and higher pressure melting conditions in the mantle wedge in comparison to typical arc systems. Liberation of volatiles in the downgoing plate is strongly dependant on the shear stress on the fault, but is predicted to occur within the seismogenic zone for shear stresses greater than 33 MPa. Similarly, early loss of volatiles is predicted over a wide range of convergence rates, plate dips, and convergence angles. These results are shown to be robust for realistic ranges of slab dip, convergence angle, and shear stress, suggesting that volatile-poor melt generation is a characteristic of modern and ancient arc systems that involve subduction of young oceanic lithosphere.     

Large-Eddy Simulation of Flow and Pollutant Dispersion in High-Aspect-Ratio Urban Street Canyons with Wall Model

A large-eddy simulation (LES) with a one-equation subgrid-scale (SGS) model was developed to investigate the flow field and pollutant dispersion inside street canyons of high aspect ratio (AR). A 1/7th power-law wall model was implemented near rigid walls to mitigate the demanding near-wall resolution requirements in LES. This LES model had been extensively validated against experimental results for street canyons of AR = 1 and 2 before it was applied to the cases of AR = 3 and 5. A ground-level passive pollutant line source, located in the middle of the street, was used to simulate vehicular emissions. Three and five vertically aligned primary recirculations were developed in the street canyons of AR 3 and 5, respectively. The ground-level mean wind speed was less than 0.5% of the free stream value, which makes it difficult for the pollutant to be transported upward for removal. High pollutant concentration and variance were found near the buildings where the air flow is upwards. It was found that the velocity fluctuation, pollutant concentration and variance were all closely related to the interactions between the primary recirculations and/or the free surface layer. Several quantities, which are non-linear functions of AR, were introduced to quantify the air quality in street canyons of different configurations.     

Assessing Climate Change Implications for Water Resources Planning

Numerous recent studies have shown that existing water supply systems are sensitive to climate change. One apparent implication is that water resources planning methods should be modified accordingly. Few of these studies, however, have attempted to account for either the chain of uncertainty in projecting water resources system vulnerability to climate change, or the adaptability of system operation resulting from existing planning strategies. Major uncertainties in water resources climate change assessments lie in a) climate modeling skill; b) errors in regional downscaling of climate model predictions; and c) uncertainties in future water demands. A simulation study was designed to provide insight into some aspects of these uncertainties. Specifically, the question that is addressed is whether a different decision would be made in a reservoir reallocation decision if knowledge about future climate were incorporated (i.e., would planning based on climate change information be justified?). The case study is possible reallocation of flood storage to conservation (municipal water supply) on the Green River, WA. We conclude that, for the case study, reservoir reallocation decisions and system performance would not differ significantly if climate change information were incorporated in the planning process.