<Emphasis Type="Italic">Deepwater Horizon</Emphasis> Oil Spill Impacts on Salt Marsh Fiddler Crabs (<Emphasis Type="Italic">Uca</Emphasis> spp.)

The Deepwater Horizon oil spill was the largest marine oil spill in US waters to date and one of the largest worldwide. Impacts of this spill on salt marsh vegetation have been well documented, although impacts on marsh macroinvertebrates have received less attention. To examine impacts of the oil spill on an important marsh invertebrate and ecosystem engineer, we conducted a meta-analysis on fiddler crabs (Uca spp.) using published sources and newly available Natural Resources Damage Assessment (NRDA) and Gulf of Mexico Research Initiative (GoMRI) data. Fiddler crabs influence marsh ecosystem structure and function through their burrowing and feeding activities and are key prey for a number of marsh and estuarine predators. We tested the hypothesis that the spill affected fiddler crab burrow density (crab abundance), burrow diameter (crab size), and crab species composition. Averaged across multiple studies, sites, and years, our synthesis revealed a negative effect of oiling on all three metrics. Burrow densities were reduced by 39 % in oiled sites, with impacts and incomplete recovery observed over 2010–2014. Burrow diameters were reduced from 2010 to 2011, but appeared to have recovered by 2012. Fiddler crab species composition was altered through at least 2013 and only returned to reference conditions where marsh vegetation recovered, via restoration planting in one case. Given the spatial and temporal extent of data analyzed, this synthesis provides compelling evidence that the Deepwater Horizon spill suppressed populations of fiddler crabs in oiled marshes, likely affecting other ecosystem attributes, including marsh productivity, marsh soil characteristics, and associated predators.     

Present-day oxidative subsidence of organic soils and mitigation in the Sacramento-San Joaquin Delta,California, USA

Subsidence of organic soils in the Sacramento-San Joaquin Delta threatens sustainability of the California (USA) water supply system and agriculture. Land-surface elevation data were collected to assess present-day subsidence rates and evaluate rice as a land use for subsidence mitigation. To depict Delta-wide present-day rates of subsidence, the previously developed SUBCALC model was refined and calibrated using recent data for CO₂ emissions and land-surface elevation changes measured at extensometers. Land-surface elevation change data were evaluated relative to indirect estimates of subsidence and accretion using carbon and nitrogen flux data for rice cultivation. Extensometer and leveling data demonstrate seasonal variations in land-surface elevations associated with groundwater-level fluctuations and inelastic subsidence rates of 0.5–0.8 cm yr^–1. Calibration of the SUBCALC model indicated accuracy of ±0.10 cm yr^–1 where depth to groundwater, soil organic matter content and temperature are known. Regional estimates of subsidence range from <0.3 to >1.8 cm yr^–1. The primary uncertainty is the distribution of soil organic matter content which results in spatial averaging in the mapping of subsidence rates. Analysis of leveling and extensometer data in rice fields resulted in an estimated accretion rate of 0.02–0.8 cm yr^–1. These values generally agreed with indirect estimates based on carbon fluxes and nitrogen mineralization, thus preliminarily demonstrating that rice will stop or greatly reduce subsidence. Areas below elevations of –2 m are candidate areas for implementation of mitigation measures such as rice because there is active subsidence occurring at rates greater than 0.4 cm yr^–1.     

Identification of Uranyl Minerals Using Oxygen K‐Edge X‐Ray Absorption Spectroscopy

Although most of the world's uranium exists as pitchblende or uraninite, this mineral can be weathered to a great variety of secondary uranium minerals, most containing the uranyl cation. Anthropogenic uranium compounds can also react in the environment, leading to spatial–chemical alterations that could be useful for nuclear forensics analyses. Soft X‐ray absorption spectroscopy (XAS) has the advantages of being non‐destructive, element‐specific and sensitive to electronic and physical structure. The soft X‐ray probe can also be focused to a spot size on the order of tens of nanometres, providing chemical information with high spatial resolution. However, before XAS can be applied at high spatial resolution, it is necessary to find spectroscopic signatures for a variety of uranium compounds in the soft X‐ray spectral region. To that end, we collected the near edge X‐ray absorption fine structure (NEXAFS) spectra of a variety of common uranyl‐bearing minerals, including uranyl carbonates, oxyhydroxides, phosphates and silicates. We find that uranyl compounds can be distinguished by class (carbonate, oxyhydroxide, phosphate or silicate) based on their oxygen K‐edge absorption spectra. This work establishes a database of reference spectra for future spatially resolved analyses. We proceed to show scanning X‐ray transmission microscopy (STXM) data from a schoepite particle in the presence of an unknown contaminant.     

High Precision Sr‐Nd‐Hf‐Pb Isotopic Compositions of USGS Reference Material BCR‐2

The Lamont‐Doherty Earth Observatory radiogenic isotope group has been systematically measuring Sr‐Nd‐Pb‐Hf isotopes of USGS reference material BCR‐2 (Columbia River Basalt 2), as a chemical processing and instrumental quality control monitor for isotopic measurements. BCR‐2 is now a widely used geochemical inter‐laboratory reference material (RM), with its predecessor BCR‐1 no longer available. Recognising that precise and accurate data on RMs is important for ensuring analytical quality and for comparing data between different laboratories, we present a compilation of multiple digestions and analyses made on BCR‐2 during the first author's dissertation research. The best estimates of Sr, Nd and Hf isotope ratios and measurement reproducibilities, after filtering at the 2s level for outliers, were ⁸⁷Sr/⁸⁶Sr = 0.705000 ± 11 (2s, 16 ppm, n = 21, sixteen digestions, one outlier), ¹⁴³Nd/¹⁴⁴Nd = 0.512637 ± 13 (2s, 25 ppm, n = 27, thirteen digestions, one outlier) and ¹⁷⁶Hf/¹⁷⁷Hf = 0.282866 ± 11 (2s, 39 ppm, n = 25, thirteen digestions, no outliers). Mean Nd and Hf values were within error of those reported by Weis et al. (2006, 2007) in their studies of RMs; mean Sr values were just outside the 2s uncertainty range of both laboratories. Moreover, a survey of published Sr‐Nd‐Hf data shows that our results fall within the range of reported values, but with a smaller variability. Our Pb isotope results on acid leached BCR‐2 aliquots (n = 26, twelve digestions, two outliers) were ²⁰⁶Pb/²⁰⁴Pb = 18.8029 ± 10 (2s, 55 ppm), ²⁰⁷Pb/²⁰⁴Pb = 15.6239 ± 8 (2s, 52 ppm), ²⁰⁸Pb/²⁰⁴Pb = 38.8287 ± 25 (2s, 63 ppm). We confirm that unleached BCR‐2 powder is contaminated with Pb, and that sufficient leaching prior to digestion is required to achieve accurate values for the uncontaminated Pb isotopic compositions.     

Adjoint analysis of Buckley-Leverett and two-phase flow equations

This paper analyzes the adjoint equations and boundary conditions for porous media flow models, specifically the Buckley-Leverett equation, and the compressible two-phase flow equations in mass conservation form. An adjoint analysis of a general scalar hyperbolic conservation law whose primal solutions include a shock jump is initially presented, and the results are later specialized to the Buckley-Leverett equation. The non-convexity of the Buckley-Leverett flux function results in adjoint characteristics that are parallel to the shock front upstream of the shock and emerge from the shock front downstream of the shock. Thus, in contrast to the behavior of Burgers’ equation where the adjoint is continuous at a shock, the Buckley-Leverett adjoint, in general, contains a discontinuous jump across the shock. Discrete adjoint solutions from space-time discontinuous Galerkin finite element approximations of the Buckley-Leverett equation are shown to be consistent with the derived closed-form analytical solutions. Furthermore, a general result relating the adjoint equations for different (though equivalent) primal equations is used to relate the two-phase flow adjoints to the Buckley-Leverett adjoint. Adjoint solutions from space-time discontinuous Galerkin finite element approximations of the two-phase flow equations are observed to obey this relationship.     

Upstream—Downstream Shifts in Peak Recruitment of the Native Olympia Oyster in San Francisco Bay During Wet and Dry Years

Understanding the conditions that drive variation in recruitment of key estuarine species can be important for effective conservation and management of their populations. The Olympia oyster (Ostrea lurida) is native to the Pacific coast of North America and has been a target of conservation efforts, though relatively little information on larval recruitment exists across much of its range. This study examined the recruitment of Olympia oysters at biweekly to monthly intervals at four sites in northern San Francisco Bay from 2010 to 2015 (except 2013). Mean monthly temperatures warmed at all sites during the study, while winter (January–April) mean monthly salinity decreased significantly during a wet year (2011), but otherwise remained high as a result of a drought. A recurring peak in oyster recruitment was identified in mid-estuary, in conditions corresponding to a salinity range of 25–30 and >16 °C at the time of settlement (April–November). Higher average salinities and temperatures were positively correlated with greater peak recruitment. Interannual variation in the timing of favorable conditions for recruitment at each site appears to explain geographic and temporal variation in recruitment onset. Higher winter/spring salinities and warmer temperatures at the time of recruitment corresponded with earlier recruitment onset within individual sites. Across all sites, higher winter/spring salinities were also correlated with earlier onset and earlier peak recruitment. Lower winter salinities during 2011 also resulted in a downstream shift in the location of peak recruitment.     

Defining the Potential of Nanoscale Re‐Os Isotope Systematics Using Atom Probe Microscopy

Atom probe microscopy (APM) is a relatively new in situ tool for measuring isotope fractions from nanoscale volumes (< 0.01 μm³). We calculate the theoretical detectable difference of an isotope ratio measurement result from APM using counting statistics of a hypothetical data set to be ± 4δ or 0.4% (2s). However, challenges associated with APM measurements (e.g., peak ranging, hydride formation and isobaric interferences), result in larger uncertainties if not properly accounted for. We evaluate these factors for Re‐Os isotope ratio measurements by comparing APM and negative thermal ionisation mass spectrometry (N‐TIMS) measurement results of pure Os, pure Re, and two synthetic Re‐Os‐bearing alloys from Schwander et al. (2015, Meteoritics and Planetary Science, 50, 893) [the original metal alloy (HSE) and alloys produced by heating HSE within silicate liquid (SYN)]. From this, we propose a current best practice for APM Re‐Os isotope ratio measurements. Using this refined approach, mean APM and N‐TIMS ¹⁸⁷Os/¹⁸⁹Os measurement results agree within 0.05% and 2s (pure Os), 0.6–2% and 2s (SYN) and 5–10% (HSE). The good agreement of N‐TIMS and APM ¹⁸⁷Os/¹⁸⁹Os measurements confirms that APM can extract robust isotope ratios. Therefore, this approach permits nanoscale isotope measurements of Os‐bearing alloys using the Re‐Os geochronometer that could not be measured by conventional measurement principles.     

Pattern,style and timing of British–Irish Ice Sheet advance and retreat over the last 45 000 years: evidence from NW Scotland and the adjacent continental shelf

Predicting the future response of ice sheets to climate warming and rising global sea level is important but difficult. This is especially so when fast-flowing glaciers or ice streams, buffered by ice shelves, are grounded on beds below sea level. What happens when these ice shelves are removed? And how do the ice stream and the surrounding ice sheet respond to the abruptly altered boundary conditions? To address these questions and others we present new geological, geomorphological, geophysical and geochronological data from the ice-stream-dominated NW sector of the last British–Irish Ice Sheet (BIIS). The study area covers around 45 000 km² of NW Scotland and the surrounding continental shelf. Alongside seabed geomorphological mapping and Quaternary sediment analysis, we use a suite of over 100 new absolute ages (including cosmogenic-nuclide exposure ages, optically stimulated luminescence ages and radiocarbon dates) collected from onshore and offshore, to build a sector-wide ice-sheet reconstruction combining all available evidence with Bayesian chronosequence modelling. Using this information we present a detailed assessment of ice-sheet advance/retreat history, and the glaciological connections between different areas of the NW BIIS sector, at different times during the last glacial cycle. The results show a highly dynamic, partly marine, partly terrestrial, ice-sheet sector undergoing large size variations in response to sub-millennial-scale climatic (Dansgaard–Oeschger) cycles over the last 45 000 years. Superimposed on these trends we identify internally driven instabilities, operating at higher frequency, conditioned by local topographic factors, tidewater dynamics and glaciological feedbacks during deglaciation. Specifically, our new evidence indicates extensive marine-terminating ice-sheet glaciation of the NW BIIS sector during Greenland Stadials 12 to 9 – prior to the main ‘Late Weichselian’ ice-sheet glaciation. After a period of restricted glaciation, in Greenland Interstadials 8 to 6, we find good evidence for rapid renewed ice-sheet build-up in NW Scotland, with the Minch ice-stream terminus reaching the continental shelf edge in Greenland Stadial 5, perhaps only briefly. Deglaciation of the NW sector took place in numerous stages. Several grounding-zone wedges and moraines on the mid- and inner continental shelf attest to significant stabilizations of the ice-sheet grounding line, or ice margin, during overall retreat in Greenland Stadials 3 and 2, and to the development of ice shelves. NW Lewis was the first substantial present-day land area to deglaciate, in the first half of Greenland Stadial 3 at a time of globally reduced sea-level c. 26 kabp , followed by Cape Wrath at c. 24 kabp. The topographic confinement of the Minch straits probably promoted ice-shelf development in early Greenland Stadial 2, providing the ice stream with additional support and buffering it somewhat from external drivers. However, c. 20–19 kabp , as the grounding-line migrated into shoreward deepening water, coinciding with a marked change in marine geology and bed strength, the ice stream became unstable. We find that, once underway, grounding-line retreat proceeded in an uninterrupted fashion with the rapid loss of fronting ice shelves – first in the west, then the east troughs – before eventual glacier stabilization at fjord mouths in NW Scotland by ~17 kabp. Around the same time, ~19–17 kabp , ice-sheet lobes readvanced into the East Minch – possibly a glaciological response to the marine-instability-triggered loss of adjacent ice stream (and/or ice shelf) support in the Minch trough. An independent ice cap on Lewis also experienced margin oscillations during mid-Greenland Stadial 2, with an ice-accumulation centre in West Lewis existing into the latter part of Heinrich Stadial 1. Final ice-sheet deglaciation of NW mainland Scotland was punctuated by at least one other coherent readvance at c. 15.5 kabp , before significant ice-mass losses thereafter. At the glacial termination, c. 14.5 kabp , glaciers fed outwash sediment to now-abandoned coastal deltas in NW mainland Scotland around the time of global Meltwater Pulse 1A. Overall, this work on the BIIS NW sector reconstructs a highly dynamic ice-sheet oscillating in extent and volume for much of the last 45 000 years. Periods of expansive ice-sheet glaciation dominated by ice-streaming were interspersed with periods of much more restricted ice-cap or tidewater/fjordic glaciation. Finally, this work indicates that the role of ice streams in ice-sheet evolution is complex but mechanistically important throughout the lifetime of an ice sheet – with ice streams contributing to the regulation of ice-sheet health but also to the acceleration of ice-sheet demise via marine ice-sheet instabilities.     

Correction: Hydraulic tomography analysis of municipal-well operation data with geology-based groundwater models

Hydrogeology Journal - A Correction to this paper has been published: https://doi.org/10.1007/s10040-021-02351-x     

An occurrence of grandidierite,kornerupine, and tourmaline in southeastern Ontario,Canada

Grandidierite, kornerupine, and tourmaline occur in high-grade pelitic gneisses from southeastern Ontario, Canada. The kornerupine occurs in quartz-bearing layers associated with biotite, cordierite, garnet, ilmenite, K-feldspar, magnetite, quartz, and, less commonly, sillimanite. Grandidierite is found in quartz-poor, cordierite+sillimanite segregations in contact with biotite, cordierite, ilmenite, K-feldspar, magnetite, sillimanite, and, more rarely, garnet. Tourmaline is sporadically distributed in all compositional layers, but is not in contact with the other borosilicates. There is no textural evidence for a reactive relationship among the three borosilicates. Neither chemical or textural equilibrium has been achieved on the scale of a thin section.It is proposed that the granite, K-feldspar-rich leucosomes, and different borosilicate assemblages in adjacent compositional layers evolved along a path of decreasing pressure and increasing temperature. The P-T path intersected a series of dehydration and melting reactions. This P-T path indicates that uplift had occurred before cooling had started and before the maximum temperature was reached. Corona and symplectite textures developed at various times during uplift both before and after cooling had started.