Complex zoning between super-calcic pigeonite and augite from the Graveyard Point sill, Oregon: a record of the interplay between bulk and interstitial liquid fractionation

The 150 m thick late Miocene Graveyard Point sill (GPS) is situated at the Idaho-Oregon border near the southwestern edge of the western Snake River Plain. It records from bottom to top continuous fractional crystallization of a tholeiitic parent magma (lower chilled border, FeO/(FeO+MgO) = 0.59, Ni = 90 ppm) towards granophyres (late pods and dikes, FeO/(FeO+MgO) = 0.98, 78 wt% SiO₂ 3.5 wt% K₂O, <4 ppm Ni) showing a typical trend of Fe and P enrichment. Fractionating minerals are olivine (Fo₇₉-Fo₂), augite (X _Fe = 0.18−0.95), feldspars (An₈₀Or₁-An₁Or₆₂), Fe-Ti oxides (Ti-rich magnetite and ilmenite), apatite and in two samples super-calcic pigeonite (Wo_18–28 Fs_41–54). The granophyres may bear some quartz. Compositionally zoned minerals record a large interval of the fractionation process in every single sample, but this interval changes with stratigraphic height. In super-calcic pigeonite-bearing samples, olivine is scarce or lacking and because super-calcic pigeonite occurs as characteristic overgrowths on augite, its formation is interpreted to be related to the schematic reaction: augite + olivine (component in melt) + SiO₂ (in melt) = pigeonite, that defines the cotectic between augite and pigeonite in olivine-saturated basaltic systems. Line measurements with the electron microprobe reveal that the transition from augite to super-calcic pigeonite is continuous. However, some crystals show an abrupt “reversal” towards augite after super-calcic pigeonite growth. Two processes compete with each other in the GPS: fractional crystallization of the bulk liquid (the bulk melt separates from solids and interstitial liquids in the solidification front) and fractional crystallization of interstitial melt in the solidification front itself. Interplay between those two processes is proposed to account for the observed variations in mineral chemistry and mineral textures. Received: 25 November 1998 / Accepted: 14 June 1999     

QUANTIFYING AND MITIGATING FLOOD RISK IN RAPIDLY AGGRADING FLUVIAL SYSTEMS

ILOWERGUADALUPERIVER,SANJOSECALIFORNIATheLowerGuadalupeRiverProjectAreadoesnothavetheabilitytoconveytheflowofthe100-yearfloodevent.LocatedbetweentheSanrranciscoBayandDowntownSanJose,theLowerGuadalupeRiverProjectareaissubjecttomanynaturalandman-inducedagentsandnaturalphysicalprocessesthataffectthehydraulicandsedimenttransportcharacteristicsoftheRiver.Themostsignificantagentsandprocessesincludethelocalgeology,regionalseismicity,tidalprocesses,pastsubsidenceintheregion,risings…     

Architecture of a tide-influenced river delta in the Frontier Formation of central Wyoming, USA

The Frewens sandstone is composed of two elongate tide-influenced sandstone bodies that are positioned directly above and slightly landward of a more wave-influenced lobate sandstone. The 20-km-long, 3-km-wide Frewens sandstone bodies coarsen upwards and fine away from their axes, have gradational bases and margins and have eroded tops abruptly overlain by marine shales. These sandstones are superbly exposed in large cliffs on the banks of the South Fork of the Powder River in central Wyoming, USA. The deposits change upwards from thinly interbedded sandstones and mudstones to metre-thick heterolithic cross-strata and, finally, to metres-thick sandstone-dominated cross-strata. There is abundant evidence for tidal modulation of depositional flows; however, palaeocurrents were strongly ebb-dominated and nearly parallel the trend of sandstone-body elongation. Detailed mapping of stratal geometry and facies across these exposures shows a complex internal architecture. Large-scale bedding units within sandstone bodies are defined by alternations in facies, bed thickness and the abundance of shales. Such bedsets are inclined (5°–15°) in walls oriented parallel to palaeoflow and gradually decrease in dip over hundreds of metres as they extend from the sandstone-dominated deposits higher in a sandstone body to muddier deposits lower in the body. Where viewed perpendicular to palaeoflow, bedsets are 100-metre-wide lenses that shingle off the sandstone-body axis towards its margins. The sandstone bodies are interpreted as sand ridge deposits formed on the shoreface of a tide-influenced river delta. Metres-thick cross-strata in the upper parts of sandstone bodies resemble deposits of bars (sandwaves) formed where tidal currents moved across shallows and the tops of tidal ridges. Heterolithic deposits lower in sandstone bodies record fluctuating currents caused by ebb and flood tides and varying river discharge. Erosion surfaces capping sandstone bodies record tidal ravinement. The tidal ridges were abandoned following transgression and covered with marine mud as waters deepened.     

A Generic Closed-Loop Contaminant Treatment System Package for MODFLOW 6

We present a contaminant treatment system (CTS) package for MODFLOW 6 that facilitates the simulation of pump-and-treat systems for groundwater remediation. Using the “nonintrusive” MODFLOW 6 application programming interface (API) capability, the CTS package can balance flows between extraction and injection wells within the outer flow solution loop and applies blended concentration/mass treatment efficiency within the outer transport solution loop. The former can be important when the requested extraction rate cannot be satisfied by the current simulated groundwater system conditions, while the latter can be important for simulating incomplete/imperfect treatment schemes. Furthermore, the CTS package allows users to temporally vary all aspects of a simulated CTS system, including the configuration and location of injection and extraction wells, and the CTS efficiency. This flexibility combined with the API-based implementation provide a generic and general CTS package that can be applied across the wide range of MODFLOW 6 simulation options and that evolves in step with MODFLOW 6 code modifications and advancements without needing to update the CTS package itself.     

Well tie for broadband seismic data

The seismic industry is increasingly acquiring broadband data in order to reap the benefits of extra low‐ and high‐frequency contents. At the low end, as the sharp low‐cut decay gets closer to zero frequency, it becomes harder for a well tie to estimate the low‐frequency response correctly. The fundamental difficulty is that well logs are too short to allow accurate estimation of the long‐period content of the data. Three distinctive techniques, namely parametric constant phase, frequency‐domain least squares with multi‐tapering, and Bayesian time domain with broadband priors, are introduced in this paper to provide a robust solution to the wavelet estimation problem for broadband seismic data. Each of these techniques has a different mathematical foundation that would enable one to explore a wide range of solutions that could be used on a case‐by‐case basis depending on the problem at hand. A case study from the North West Shelf Australia is used to analyse the performance of the proposed techniques. Cross‐validation is proposed as a robust quality control measure for evaluating well‐tie applications. It is observed that when the seismic data are carefully processed, then the constant phase approach would likely offer a good solution. The frequency‐domain method does not assume a constant phase. This flexibility makes it prone to over‐fitting when the phase is approximately constant. Broadband priors for the time‐domain least‐squares method are found to perform well in defining low‐frequency side lobes to the wavelet.     

Towards Improved Remedial Outcomes in Categorical Aquifers with an Iterative Ensemble Smoother

Categorical parameter distributions consisting of geologic facies with distinct properties, for example, high-permeability channels embedded in a low-permeability matrix, are common at contaminated sites. At these sites, low-permeability facies store solute mass, acting as secondary sources to higher-permeability facies, sustaining concentrations for decades while increasing risk and cleanup costs. Parameter estimation is difficult in such systems because the discontinuities in the parameter space hinder the inverse problem. This paper presents a novel approach based on Traveling Pilot Points (TRIPS) and an iterative ensemble smoother (IES) to solve the categorical inverse problem. Groundwater flow and solute transport in a hypothetical aquifer with a categorical parameter distribution are simulated using MODFLOW 6. Heads and concentrations are recorded at multiple monitoring locations. IES is used to generate posterior ensembles assuming a TRIPS prior and an approximate multi-Gaussian prior. The ensembles are used to predict solute concentrations and mass into the future. The evaluation also includes an assessment of how the number of measurements and the choice of the geological prior determine the characteristics of the posterior ensemble and the resulting predictions. The results indicate that IES was able to efficiently sample the posterior distribution and showed that even with an approximate geological prior, a high degree of parameterization and history matching could lead to parameter ensembles that can be useful for making certain types of predictions (heads, concentrations). However, the approximate geological prior was insufficient for predicting mass. The analysis demonstrates how decision-makers can quantify uncertainty and make informed decisions with an ensemble-based approach.     

Ecohydrologic feedbacks controlling sizes of cypress wetlands in a patterned karst landscape

Many landforms on Earth are profoundly influenced by biota. In particular, biota play a significant role in creating karst biogeomorphology, through biogenic CO₂ accelerating calcite weathering. In this study, we explore the ecohydrologic feedback mechanisms that have created isolated depressional wetlands on exposed limestone bedrock in South Florida – Big Cypress National Preserve –as a case study for karst biogeomorphic processes giving rise to regularly patterned landscapes. Specifically, we are interested in: (1) whether cypress depressions on the landscape have reached (or will reach) equilibrium size; (2) if so, what feedback mechanisms stabilize the size of depressions; and (3) what distal interactions among depressions give rise to the even distribution of depressions in the landscape. We hypothesize three feedback mechanisms controlling the evolution of depressions and build a numerical model to evaluate the relative importance of each mechanism. We show that a soil cover feedback (i.e. a smaller fraction of CO₂ reaches the bedrock surface for weathering as soil cover thickens) is the major feedback stabilizing depressions, followed by a biomass feedback (i.e. inhibited biomass growth with deepening standing water and extended inundation period as depressions expand in volume). Strong local positive feedback between the volume of depressions and rate of volume expansion and distal negative feedback between depressions competing for water likely lead to the regular patterning at the landscape scale. The individual depressions, however, are not yet in steady state but would be in ~0.2–0.4 million years. This represents the first study to demonstrate the decoupling of landscape-scale self-organization and the self-organization of its constituent agents. © 2018 John Wiley & Sons, Ltd.     

Effectiveness of Monitored Natural Attenuation (MNA) as a Groundwater Remedy for Arsenic in Phosphatic Wastes

A study was conducted to evaluate monitored natural attenuation (MNA) as a remedy for arsenic in groundwater at a former phosphate mining and manufacturing facility. The mineralogy, speciation, and lability of arsenic in phosphatic wastes present in soils were characterized using sequential extraction procedures, leaching experiments, batch adsorption tests, and microchemical speciation analysis. A PHREEQC-based reactive transport model was also parameterized using these laboratory results, and it was used to evaluate the importance of identified attenuation mechanisms on arsenic concentrations along a vertical flow path from a shallow, alluvial aquifer to the underlying Floridan aquifer. Arsenic was found to occur in several chemical forms in phosphatic wastes, including unstable sulfide minerals, adsorbed surface complexes, and relatively insoluble phosphate and oxide minerals. Most arsenic was associated with stable minerals. The reactive transport model predicted that historical leaching of solid-phase waste materials in soils would not have generated enough arsenic to explain the concentrations observed in downgradient groundwater; instead, the source of arsenic to groundwater was likely acidic and saline process water that infiltrated though unlined ponds and ditches during historical manufacturing operations. A key factor affecting the long-term effectiveness of natural attenuation of arsenic in groundwater is the occurrence and stability of iron oxyhydroxides in aquifer sediments. According to laboratory and reactive transport model results, sufficient levels were found to be present at the site to effectively limit arsenic migration at concentrations exceeding drinking water standards in the future in the Floridan aquifer. This study presents the geochemical evaluations that are needed to satisfy EPA guidelines on determining whether or not MNA is an acceptable remedy for a site. It specifically details the characterization and modeling that were used to demonstrate effectiveness at a site where MNA was ultimately selected as the remedy for arsenic in groundwater.     

Recent shifts in coastline change and shoreline stabilization linked to storm climate change

Since cuspate coastlines are especially sensitive to changes in wave climate, they serve as potential indicators of initial responses to changing wave conditions. Previous work demonstrates that Cape Hatteras and Cape Lookout, North Carolina, which are largely unaffected by shoreline stabilization efforts, have become increasingly asymmetric over the past 30 years, consistent with model predictions for coastline response to increases in Atlantic Ocean summer wave heights and resulting changes in the distribution of wave‐approach angles. Historic and recent shoreline change observations for Cape Fear, North Carolina, and model simulations of coastline response to an increasingly asymmetric wave climate in the presence of beach nourishment, produce comparable differences in shoreline change rates in response to changes in wave climate. Results suggest that the effect of beach nourishment is to compensate for – and therefore to mask – natural responses to wave climate change that might otherwise be discernible in patterns of shoreline change alone. Therefore, this case study suggests that the effects of wave climate change on human‐modified coastlines may be detectable in the spatial and temporal patterns of shoreline stabilization activities. Similar analyses of cuspate features in areas where the change in wave climate is less pronounced (i.e. Fishing Point, Maryland/Virginia) and where local geology appears to exert control on coastline shape (i.e. Cape Canaveral, Florida), suggest that changes in shoreline configuration that may be arising from shifting wave climate are currently limited to sandy wave‐dominated coastlines where the change in wave climate has been most pronounced. However, if hurricane‐generated wave heights continue to increase, large‐scale shifts in patterns of erosion and accretion will likely extend beyond sensitive cuspate features as the larger‐scale coastline shape comes into equilibrium with changing wave conditions. Copyright © 2014 John Wiley & Sons, Ltd.