The distribution of nitrate N/N in marine sediments and the impact of benthic nitrogen loss on the isotopic composition of oceanic nitrate

We report ¹⁵N/¹⁴N ratios of porewater nitrate in sediments from the Bering Sea basin, where microbial nitrate reduction has been identified as a significant sink for fixed nitrogen (N). Strong ¹⁵N enrichment in porewater nitrate is observed as one goes deeper in the sediments and nitrate concentration decreases (δ¹⁵N generally reaches 25-35‰). Analysis of profiles with a one-dimensional diffusion-reaction model yields organism-scale isotope effects for dissimilatory nitrate reduction (ε_cell) of 11‰ to 30‰, in the same range as measured in previous studies of cultures and the marine and lacustrine water column. Estimates of ε_cell, while uncertain, show a negative correlation with bottom water [O₂]; we propose that this relates to the at the depth of denitrification. The N isotope effect at the scale of nitrate sediment-water exchange (ε_app) is ∼0‰ in two unreactive deep sites and is typically <3‰ at more reactive sites at various depths. ε_app is much lower than ε_cell because nitrate consumption is nearly complete at the sediment depth of denitrification, minimizing the escape of ¹⁵N-enriched nitrate from the sediments. In reactive sediments, this is due to rapid denitrification, while in less reactive sediments, it is due to greater diffusive distances for nitrate to the depth of denitrification. The data suggest that low bottom water [O₂] tends to yield more complete expression of ε_cell at the sediment-water scale, due to higher at the depth of denitrification. While porewater ammonium-N isotopes were not measured, our porewater model suggests that, in sediments with high organic matter supply and/or low-[O₂] bottom waters, the efflux and subsequent oxidation of ammonium enriched in ¹⁵N by incomplete nitrification can significantly enhance the total net isotope effect of sedimentary N loss (ε_sed, equivalent to ε_app but including ammonium fluxes). Model analysis of representative sedimentary environments suggests a global mean ε_sed of ∼4‰ (∼2‰ if restricted to seafloor below 1 km depth).     

Contributions and unrealized potential contributions of cosmogenic-nuclide exposure dating to glacier chronology, 1990–2010

This paper reviews the application of cosmogenic-nuclide exposure dating to glacier chronology. Exposure dating of glacial landforms has made an outsize impact on this field because the technique filled an obvious need that had already been recognized by glacial geologists. By now, hundreds of studies have used cosmogenic-nuclide exposure dating to date glacial deposits, and in fact it is rare to find a study of glacial geology or glacier chronology, or any paleoclimate synthesis that makes use of such studies, that does not involve exposure dating. These developments have resulted in major contributions to glacier chronology and paleoclimate, in particular i) reconstructing Antarctic ice sheet change, ii) establishing the chronology of late Pleistocene and Holocene glacier change in mountain regions where it was previously unknown; iii) establishing the broad chronological outlines of mountain glaciations prior to the Last Glacial Maximum; and iv) gaining insight into subglacial erosional processes through the observation that many glaciated surfaces preserve cosmogenic-nuclide inventories from long past ice-free periods as well as the present one. An important potential future contribution will be the application of the large data set of exposure-dated glacier chronologies to better understand global and regional climate dynamics during Lateglacial and Holocene millennial-scale climate changes. However, this contribution cannot be realized without significant progress in two areas: i) understanding and accounting for geologic processes that cause apparent exposure ages on glacial landforms to differ from the true age of the landform, and ii) minimizing systematic uncertainties in exposure ages that stem from cosmogenic-nuclide production-rate estimates and scaling schemes. At present there exists an enormous data set of exposure ages on glacial deposits, but these data cannot be used to their full potential in paleoclimate syntheses due to an inadequate understanding of geologic scatter and production-rate uncertainties. The intent of this paper is to highlight this situation and suggest some strategies for realizing this potential.     

Exploring the limits of identifying sub-pixel thermal features using ASTER TIR data

Understanding the characteristics of volcanic thermal emissions and how they change with time is important for forecasting and monitoring volcanic activity and potential hazards. Satellite instruments view volcanic thermal features across the globe at various temporal and spatial resolutions. Thermal features that may be a precursor to a major eruption, or indicative of important changes in an on-going eruption can be subtle, making them challenging to reliably identify with satellite instruments. The goal of this study was to explore the limits of the types and magnitudes of thermal anomalies that could be detected using satellite thermal infrared (TIR) data. Specifically, the characterization of sub-pixel thermal features with a wide range of temperatures is considered using ASTER multispectral TIR data. First, theoretical calculations were made to define a “thermal mixing detection threshold” for ASTER, which quantifies the limits of ASTER's ability to resolve sub-pixel thermal mixing over a range of hot target temperatures and % pixel areas. Then, ASTER TIR data were used to model sub-pixel thermal features at the Yellowstone National Park geothermal area (hot spring pools with temperatures from 40 to 90 °C) and at Mount Erebus Volcano, Antarctica (an active lava lake with temperatures from 200 to 800 °C). Finally, various sources of uncertainty in sub-pixel thermal calculations were quantified for these empirical measurements, including pixel resampling, atmospheric correction, and background temperature and emissivity assumptions.     

A Model-Averaging Method for Assessing Groundwater Conceptual Model Uncertainty

This study evaluates alternative groundwater models with different recharge and geologic components at the northern Yucca Flat area of the Death Valley Regional Flow System (DVRFS), USA. Recharge over the DVRFS has been estimated using five methods, and five geological interpretations are available at the northern Yucca Flat area. Combining the recharge and geological components together with additional modeling components that represent other hydrogeological conditions yields a total of 25 groundwater flow models. As all the models are plausible given available data and information, evaluating model uncertainty becomes inevitable. On the other hand, hydraulic parameters (e.g., hydraulic conductivity) are uncertain in each model, giving rise to parametric uncertainty. Propagation of the uncertainty in the models and model parameters through groundwater modeling causes predictive uncertainty in model predictions (e.g., hydraulic head and flow). Parametric uncertainty within each model is assessed using Monte Carlo simulation, and model uncertainty is evaluated using the model averaging method. Two model-averaging techniques (on the basis of information criteria and GLUE) are discussed. This study shows that contribution of model uncertainty to predictive uncertainty is significantly larger than that of parametric uncertainty. For the recharge and geological components, uncertainty in the geological interpretations has more significant effect on model predictions than uncertainty in the recharge estimates. In addition, weighted residuals vary more for the different geological models than for different recharge models. Most of the calibrated observations are not important for discriminating between the alternative models, because their weighted residuals vary only slightly from one model to another.     

Evaluation of Bias Associated with Capture Maps Derived from Nonlinear Groundwater Flow Models

The impact of groundwater withdrawal on surface water is a concern of water users and water managers, particularly in the arid western United States. Capture maps are useful tools to spatially assess the impact of groundwater pumping on water sources (e.g., streamflow depletion) and are being used more frequently for conjunctive management of surface water and groundwater. Capture maps have been derived using linear groundwater flow models and rely on the principle of superposition to demonstrate the effects of pumping in various locations on resources of interest. However, nonlinear models are often necessary to simulate head‐dependent boundary conditions and unconfined aquifers. Capture maps developed using nonlinear models with the principle of superposition may over‐ or underestimate capture magnitude and spatial extent. This paper presents new methods for generating capture difference maps, which assess spatial effects of model nonlinearity on capture fraction sensitivity to pumping rate, and for calculating the bias associated with capture maps. The sensitivity of capture map bias to selected parameters related to model design and conceptualization for the arid western United States is explored. This study finds that the simulation of stream continuity, pumping rates, stream incision, well proximity to capture sources, aquifer hydraulic conductivity, and groundwater evapotranspiration extinction depth substantially affect capture map bias. Capture difference maps demonstrate that regions with large capture fraction differences are indicative of greater potential capture map bias. Understanding both spatial and temporal bias in capture maps derived from nonlinear groundwater flow models improves their utility and defensibility as conjunctive‐use management tools.     

A simple algorithm for generating streamflow networks for grid‐based,macroscale hydrological models

A simple algorithm for generating streamflow networks for macroscale hydrological models (MHMs) from digital elevation models (DEMs) is presented. Typically these hydrological models are grid based, with the simulated runoff produced within each cell routed through a stream network which connects the centers of cells in the direction of the major streams. Construction of such stream networks is a time consuming task, which has generally been done by hand with the aid of maps. Results indicate that the algorithm works satisfactorily in areas of both high and low relief, and for a wide range of model cell resolutions, although some manual adjustments may be necessary. Copyright © 1999 John Wiley & Sons, Ltd.     

LNAPL Recovery Endpoints: Lessons Learnt Through Modeling,Experiments, and Field Trials

It is important to estimate what light nonaqueous phase liquid (LNAPL) recovery can be practicably achieved from subsurface environments. Over the last decade, research to address this included a broad field program, laboratory measurements and experimentation, and modeling approaches. Here, we consolidate key findings from the research in the context of current literature and understanding, with a focus on a well-validated, multiphase multicomponent modeling approach to achieve estimates of reasonable endpoints for LNAPL recovery. Simple analytical models can provide approximate saturation distributions and estimates of LNAPL recoverability via transmissivity approximation, but are insufficient to predict LNAPL saturation- and composition-based recovery endpoints for various recovery technologies. This is because they cannot account for multiphase, multicomponent fate and transport and key processes such as hysteresis. Recent advances to improve estimates of the fraction of recoverable LNAPL and its transmissivity are summarized. These advances include further development and application of a well-validated model to characterize active LNAPL recovery endpoints. We present key factors that affect the determination of LNAPL recovery endpoints, and outline how recovery endpoints are affected by natural source zone depletion (NSZD—currently gaining acceptance as a LNAPL remediation option). Major factors include geo-physical characteristics of the formation, magnitude of an LNAPL release and partitioning properties of the key LNAPL constituents of concern. Based on the capabilities of the validated model, the paper also provides a basis to optimize LNAPL recovery efforts.     

K‐Ar age determination,whole‐rock and oxygen isotope geochemistry of the post‐collisional Bizmişen and Çaltı plutons,SW Erzincan,eastern Central Anatolia,Turkey

Post‐collisional granitoid plutons intrude obducted Neo‐Tethyan ophiolitic rocks in central and eastern Central Anatolia. The Bizmişen and Çaltı plutons and the ophiolitic rocks that they intrude are overlain by fossiliferous and flyschoidal sedimentary rocks of the early Miocene Kemah Formation. These sedimentary rocks were deposited in basins that developed at the same time as tectonic unroofing of the plutons along E–W and NW–SE trending faults in Oligo‐Miocene time. Mineral separates from the Bizmişen and Çaltı plutons yield K‐Ar ages ranging from 42 to 46 Ma, and from 40 to 49 Ma, respectively. Major, trace, and rare‐earth element geochemistry as well as mineralogical and textural evidence reveals that the Bizmişen pluton crystallized first, followed at shallower depth by the Çaltı pluton from a medium‐K calcalkaline, I‐type hybrid magma which was generated by magma mixing of coeval mafic and felsic magmas. Delta ¹⁸O values of both plutons fall in the field of I‐type granitoids, although those of the Çaltı pluton are consistently higher than those of the Bizmişen pluton. This is in agreement with field observations, petrographic and whole‐rock geochemical data, which indicate that the Bizmişen pluton represents relatively uncontaminated mantle material, whereas the Çaltı pluton has a significant crustal component. Structural data indicating the middle Eocene emplacement age and intrusion into already obducted ophiolitic rocks, suggest a post‐collisional extensional origin. However, the pure geochemical discrimination diagrams indicate an arc origin which can be inherited either from the source material or from an upper mantle material modified by an early subduction process during the evolution of the Neo‐Tethyan ocean. Copyright © 2005 John Wiley & Sons, Ltd.