A tool for site prioritisation of marine protected areas under data poor conditions

A challenge of managing data poor fisheries lies in overcoming uncertainties associated with a lack of information on biological and socio-economic conditions. This paper focuses on site prioritisation for marine protected areas, a process that usually imposes moderate to heavy demands on data, time, local expertise, and funds. A fuzzy logic decision support tool for zoning marine areas that is suitable for use in data poor conditions is developed. This tool, the Protected Area Suitability Index (PASI), assesses the suitability of sites for protection based on fishers’ preferences for that site and the site’s conservation value. Only eight input attributes are required to run the PASI, which operates on a series of heuristic rules to estimate a site suitability score that ranges from 0 to 10, where 10 indicates that a site is highly suitable for being protected from fishing. Sensitivity analysis reveals that the PASI is relatively robust, and produces reliable results even as the system is degraded by the loss of data. Eighteen actively fished sites and 11 sites that are closed to fishing are assessed to evaluate how well the PASI is able to distinguish differences between sites. The estimated scores are significantly different (p<0.05), whereby non-fished sites are scored as being more suitable for protection than fished sites. The PASI can be used as a decision support tool to facilitate systematic marine spatial management under data poor conditions, especially in the task of identifying suitable sites for protection.     

Interactive Visualization to Advance Earthquake Simulation

The geological sciences are challenged to manage and interpret increasing volumes of data as observations and simulations increase in size and complexity. For example, simulations of earthquake-related processes typically generate complex, time-varying data sets in two or more dimensions. To facilitate interpretation and analysis of these data sets, evaluate the underlying models, and to drive future calculations, we have developed methods of interactive visualization with a special focus on using immersive virtual reality (VR) environments to interact with models of Earth’s surface and interior. Virtual mapping tools allow virtual “field studies” in inaccessible regions. Interactive tools allow us to manipulate shapes in order to construct models of geological features for geodynamic models, while feature extraction tools support quantitative measurement of structures that emerge from numerical simulation or field observations, thereby enabling us to improve our interpretation of the dynamical processes that drive earthquakes. VR has traditionally been used primarily as a presentation tool, albeit with active navigation through data. Reaping the full intellectual benefits of immersive VR as a tool for scientific analysis requires building on the method’s strengths, that is, using both 3D perception and interaction with observed or simulated data. This approach also takes advantage of the specialized skills of geological scientists who are trained to interpret, the often limited, geological and geophysical data available from field observations.     

Particle transport and distribution on the Mars Science Laboratory mission: Effects of triboelectric charging

We report on the nature of fine particle (<150 μm) transport under simulated martian conditions, in order to better understand the Mars Science Laboratory’s (MSL) sample acquisition, processing and handling subsystem (SA/SPaH). We find that triboelectric charging due to particle movement may have to be controlled in order for successful transport of fines that are created within the drill, processed through the Collection and Handling for In situ Martian Rock Analysis (CHIMRA) sample handing system, and delivered to the Sample Analysis at Mars (SAM) and Chemistry and Mineralogy (CheMin) instruments. These fines will be transferred from the surface material to the portioner, a 3 mm diameter, 8 mm deep distribution center where they will drop ∼2 cm to the instrument inlet funnels. In our experiments, movement of different material including terrestrial analogs and martian soil simulants (Mars Mojave Simulant - MMS) resulted in 1-7 nanocoulombs of charge to build up for several different experimental configurations. When this charging phenomenon occurs, several different results are observed including particle clumping, adherence of material on conductive surfaces, or electrostatic repulsion, which causes like-charged particles to move away from each other. This electrostatic repulsion can sort samples based upon differing size fractions, while adhesion causes particles of different sizes to bind into clods. Identifying these electrostatic effects can help us understand potential bias in the analytical instruments and to define the best operational protocols to collect samples on the surface of Mars.     

Implications of Warm Rain in Shallow Cumulus and Congestus Clouds for Large-Scale Circulations

Space-borne observations reveal that 20–40% of marine convective clouds below the freezing level produce rain. In this paper we speculate what the prevalence of warm rain might imply for convection and large-scale circulations over tropical oceans. We present results using a two-column radiative–convective model of hydrostatic, nonlinear flow on a non-rotating sphere, with parameterized convection and radiation, and review ongoing efforts in high-resolution modeling and observations of warm rain. The model experiments investigate the response of convection and circulation to sea surface temperature (SST) gradients between the columns and to changes in a parameter that controls the conversion of cloud condensate to rain. Convection over the cold ocean collapses to a shallow mode with tops near 850 hPa, but a congestus mode with tops near 600 hPa can develop at small SST differences when warm rain formation is more efficient. Here, interactive radiation and the response of the circulation are crucial: along with congestus a deeper moist layer develops, which leads to less low-level radiative cooling, a smaller buoyancy gradient between the columns, and therefore a weaker circulation and less subsidence over the cold ocean. The congestus mode is accompanied with more surface precipitation in the subsiding column and less surface precipitation in the deep convecting column. For the shallow mode over colder oceans, circulations also weaken with more efficient warm rain formation, but only marginally. Here, more warm rain reduces convective tops and the boundary layer depth—similar to Large-Eddy Simulation (LES) studies—which reduces the integrated buoyancy gradient. Elucidating the impact of warm rain can benefit from large-domain high-resolution simulations and observations. Parameterizations of warm rain may be constrained through collocated cloud and rain profiling from ground, and concurrent changes in convection and rain in subsiding and convecting branches of circulations may be revealed from a collocation of space-borne sensors, including the Global Precipitation Measurement (GPM) and upcoming Aeolus missions.     

Landscape variations in stream water SO4 and δSSO4 in a boreal stream network

Despite reduced anthropogenic deposition during the last decades, deposition sulphate may still play an important role in the biogeochemical cycles of S and many catchments may act as net sources of S that may remain for several decades. The aim of this study is to elucidate the temporal and spatial dynamics of both SO₄²⁻ and δ³⁴S_SO4 in stream water from catchments with varying percentage of wetland and forest coverage and to determine their relative importance for catchment losses of S. Stream water samples were collected from 15 subcatchments ranging in size from 3 to 6780 ha, in a boreal stream network, northern Sweden. In forested catchments (<2% wetland cover) S-SO₄²⁻ concentrations in stream water averaged 1.7 mg L⁻¹ whereas in wetland dominated catchments (>30% wetland cover) the concentrations averaged 0.3 mg L⁻¹. A significant negative relationship was observed between S-SO₄²⁻ and percentage wetland coverage (r² = 0.77, p < 0.001) and the annual export of stream water SO₄²⁻ and wetland coverage (r² = 0.76, p < 0.001). The percentage forest coverage was on the other hand positively related to stream water SO₄²⁻ concentrations and the annual export of stream water SO₄²⁻ (r² = 0.77 and r² = 0.79, respectively). The annual average δ³⁴S_SO4 value in wetland dominated streams was +7.6‰ and in streams of forested catchments +6.7‰. At spring flood the δ³⁴S_SO4 values decreased in all streams by 1‰ to 5‰. The δ³⁴S_SO4 values in all streams were higher than the δ³⁴S_SO4 value of +4.7‰ in precipitation (snow). The export of S ranged from 0.5 kg S ha⁻¹ yr⁻¹ (wetland headwater stream) to 3.8 kg S ha⁻¹ yr⁻¹ (forested headwater stream). With an average S deposition in open field of 1.3 kg S ha⁻¹ yr⁻¹ (2002-2006) the mass balance results in a net export of S from all catchments, except in catchments with >30% wetland. The high temporal and spatial resolution of this study demonstrates that the reducing environments of wetlands play a key role for the biogeochemistry of S in boreal landscapes and are net sinks of S. Forested areas, on the other hand were net sources of S.     

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.     

The yield and isotopic composition of radiolytic H2, a potential energy source for the deep subsurface biosphere

The production rate and isotopic composition of H₂ derived from radiolytic reactions in H₂O were measured to assess the importance of radiolytic H₂ in subsurface environments and to determine whether its isotopic signature can be used as a diagnostic tool. Saline and pure, aerobic and anaerobic water samples with pH values of 4, 7, and 10 were irradiated in sealed vials at room temperature with an artificial γ source, and the H₂ abundance in the headspace and its isotopic composition were measured. The H₂ concentrations were observed to increase linearly with dosage at a rate of 0.40 ± 0.04 molecules (100 eV)⁻¹ within the dosage range of 900 to 3500 Gray (Gy; Gy = 1 J Kg⁻¹) with no indication of a maximum limit on H₂ concentration. At ∼2000 Gy, the H₂ concentration varied only by 16% across the experimental range of pH, salinity, and O₂. Based upon this measured yield and H₂ yields for α and β particles, a radiolytic H₂ production rate of 10⁻⁹ to 10⁻⁴ nM s⁻¹ was estimated for the range of radioactive element concentrations and porosities typical of crustal rocks. The δD of H₂ was independent of the dosage, pH (except for pH 4), salinity, and O₂ and yielded an αD_H2O-H2 of 2.05 ± 0.07 (αD_H2O-H2 = (D/H)_H2O to (D/H)_H2), slightly less than predicted radiolytic models. Although this radiolytic fractionation value is significantly heavier than that of equilibrium isotopic exchange between H₂ and H₂O, the isotopic exchange rate between H₂ and H₂O will erase the heavy δD of radiolytic H₂ if the age of the groundwater is greater than ∼10³ to 10⁴ yr. The millimolar concentrations of H₂ observed in the groundwater of several Precambrian Shields are consistent with radiolysis of water that has resided in the subsurface for a few million years. These concentrations are well above those required to support H₂-utilizing microorganisms and to inhibit H₂-producing, fermentative microorganisms.     

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.     

Late Quaternary sea‐level change and evolution of Belfast Lough,Northern Ireland: new offshore evidence and implications for sea‐level reconstruction

The interplay of eustatic and isostatic factors causes complex relative sea‐level (RSL) histories, particularly in paraglacial settings. In this context the past record of RSL is important in understanding ice‐sheet history, earth rheology and resulting glacio‐isostatic adjustment. Field data to develop sea‐level reconstructions are often limited to shallow depths and uncertainty exists as to the veracity of modelled sea‐level curves. We use seismic stratigraphy, 39 vibrocores and 26 radiocarbon dates to investigate the deglacial history of Belfast Lough, Northern Ireland, and reconstruct past RSL. A typical sequence of till, glacimarine and Holocene sediments is preserved. Two sea‐level lowstands (both max. ?40 m) are recorded at c. 13.5 and 11.5k cal a bp . Each is followed by a rapid transgression and subsequent periods of RSL stability. The first transgression coincides temporally with a late stage of Meltwater Pulse 1a and the RSL stability occurred between c. 13.0 and c. 12.2k cal a bp (Younger Dryas). The second still/slowstand occurred between c. 10.3 and c. 11.5k cal a bp . Our data provide constraints on the direction and timing of RSL change during deglaciation. Application of the Depth of Closure concept adds an error term to sea‐level reconstructions based on seismic stratigraphic reconstructions.