A Comparative Analysis of Weights of Evidence,Evidential Belief Functions,and Fuzzy Logic for Mineral Potential Mapping Using Incomplete Data at the Scale of Investigation

Large amounts of digital data must be analyzed and integrated to generate mineral potential maps, which can be used for exploration targeting. The quality of the mineral potential maps is dependent on the quality of the data used as inputs, with higher quality inputs producing higher quality outputs. In mineral exploration, particularly in regions with little to no exploration history, datasets are often incomplete at the scale of investigation with data missing due to incomplete mapping or the unavailability of data over certain areas. It is not always clear that datasets are incomplete, and this study examines how mineral potential mapping results may differ in this context. Different methods of mineral potential mapping provide different ways of dealing with analyzing and integrating incomplete data. This study examines the weights of evidence (WofE), evidential belief function and fuzzy logic methods of mineral potential mapping using incomplete data from the Carajás mineral province, Brazil to target for orogenic gold mineralization. Results demonstrate that WofE is the best one able to predict the location of known mineralization within the study area when either complete or unacknowledged incomplete data are used. It is suggested that this is due to the use of Bayes’ rule, which can account for “missing data.” The results indicate the effectiveness of WofE for mineral potential mapping with incomplete data.     

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.     

Modeling the Population Effects of Hypoxia on Atlantic Croaker (<Emphasis Type="Italic">Micropogonias undulatus</Emphasis>) in the Northwestern Gulf of Mexico: Part 1—Model Description and Idealized Hypoxia

We developed a spatially explicit, individual-based model to analyze how hypoxia effects on reproduction, growth, and mortality of Atlantic croaker in the northwestern Gulf of Mexico lead to population-level responses. The model follows the hourly growth, mortality, reproduction, and movement of individuals on a 300 × 800 spatial grid of 1-km² cells for 140 years. Chlorophyll-a concentration, water temperature, and dissolved oxygen (DO) were specified daily for each grid cell and repeated for each year of the simulation. A bioenergetics model was used to represent growth, mortality was assumed stage- and age-dependent, and the movement behavior of juveniles and adults was modeled based on temperature and avoidance of low DO. Hypoxia effects were imposed using exposure effect submodels that converted time-varying exposures to low DO to reduced hourly growth, increased hourly mortality, and reduced annual fecundity. Results showed that 100 years of either mild or intermediate hypoxia produced small reductions in population abundance, while repeated severe hypoxia caused a 19% reduction in long-term population abundance. Relatively few individuals were exposed to low DO each hour, but many individuals experienced some exposure. The response was dominated by a 5% average reduction in annual fecundity of individuals. Under conditions of random sequences of mild, intermediate, and severe hypoxia years occurring in proportion to their historical frequency, the model predicted a 10% decrease in the long-term population abundance of croaker. A companion paper substitutes hourly DO values from a three-dimensional water quality model for the idealized hypoxia and results in a more realistic population reduction of about 25%.     

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.     

Composition,Age, and Origin of Cretaceous Granitic Magmatism on the Eastern Chukchi Peninsula

New geochronological and isotopic geochemical data are given, which make it possible to recognize two types of granitic rocks on the eastern Chukchi Peninsula. Early Cretaceous Tkachen and Dolina granitic plutons with zircon ages (U–Pb SIMS) of 119–122 and 131–136 Ma are related to the first type. They cut through Devonian–Lower Carboniferous basement rocks and are overlain by the Aptian–Albian Etelkuyum Formation. Basal units of the latter contain fragments of granitic rocks. Late Cretaceous Provideniya and Rumilet granitic plutons, which contain zircons with ages of 94 and 85 Ma (U–Pb SIMS), respectively, belong to the second type. They cut through volcanic–sedimentary rocks of the Etelkuyum and Leurvaam formations pertaining to the Okhotsk–Chukotka Volcanic Belt. In petrographic and geochemical features, the Early Cretaceous granitic rocks of the Tkachen Pluton are commensurable with I-type granites, while Late Cretaceous granite of the Rumilet Pluton is comparable to A₂-type granite. The Sr–Nd isotopic data provide evidence that from the Early Cretaceous Tkachen and Dolina plutons to the Late Cretaceous Provideniya and Rumilet plutons, the degree of crustal assimilation of suprasubduction mantle-derived melts increases up to partial melting of heterogeneous continental crust enriched in rubidium. An unconformity and various degrees of secondary alteration of volcanic–sedimentary rocks have been established in the Okhotsk–Chukotka Volcanic Belt, and this was apparently caused by transition of the tectonic setting from suprasubduction to a transform margin with local extension.     

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.     

Built environment and agricultural value at risk from Hurricane Irma flooding in Florida (USA)

Natural Hazards - Flooding is the most costly type of natural disaster, as well as the most frequent. To provide risk-based flood insurance, providers must be able to accurately determine an...     

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     

Using Mesoproterozoic sedimentary geochemistry to reconstruct basin tectonic geography and link organic carbon productivity to nutrient flux from a Northern Australian large igneous Province

The Beetaloo Sub-basin, northern Australia, is considered the main depocentre of the 1,000-km scale Mesoproterozoic Wilton package of the greater McArthur Basin – the host to one of the oldest hydrocarbon global resources. The ca. 1.40–1.31 Ga upper Roper Group and the latest Mesoproterozoic to early Neoproterozoic unnamed group of the Beetaloo Sub-basin, together, record ca. 500 million years of depositional history within the North Australia Craton. Whole-rock shale Sm–Nd and Pb isotope data from these sediments reveal sedimentary provenance and their evolution from ca. 1.35 to 0.85 Ga. Furthermore, these data, together with shale major/trace elements data from this study and pyrolysis data from previous publications, are used to develop a dynamic tectonic geography model that links the organic carbon production and burial to an enhanced weathering of nutrients from a large igneous province. The ca. 1.35–1.31 Ga Kyalla Formation of the upper Roper Group is composed of isotopically evolved sedimentary detritus that passes up, into more isotopically juvenile Pb values towards the top of the formation. The increase in juvenile compositions coincides with elevated total organic carbon (TOC) contents of these sediments. The coherently enriched juvenile compositions and TOC the upper portions of the Kyalla Formation are interpreted to reflect an increase in nutrient supply associated with the weathering of basaltic sources (e.g. phosphorous). Possible, relatively juvenile, basaltic sources include the Wankanki Supersuite in the western Musgraves and the Derim Derim–Galiwinku large igneous province (LIP). The transition into juvenile, basaltic sources directly before a supersequence-bounding unconformity is here interpreted to reflect uplift and erosion of the Derim Derim–Galiwinku LIP, rather than an influx of southern Musgrave sources. A new baddeleyite crystallisation age of 1,312.9 ± 0.7 Ma provides both a tight constraint on the age of this LIP, along with its associated magmatic uplift, as well as providing a minimum age constraint for Roper Group deposition. The unconformably overlying lower and upper Jamison sandstones are at least 300 million years younger than the Kyalla Formation and were sourced from the Musgrave Province. An up-section increase in isotopically juvenile compositions seen in these rocks is interpreted to document the progressive exhumation of the western Musgrave Province. The overlying Hayfield mudstone received detritus from both the Musgrave and Arunta regions, and its isotopic geochemistry reveals affinities with other early Neoproterozoic basins (e.g. Amadeus, Victoria and Officer basins), indicating the potential for inter-basin correlations.