Apatite as a monitor of late-stage magmatic processes at Volcán Irazú, Costa Rica

Apatite phenocrysts from the 1963 and 1723 eruptions of Irazú volcano (Costa Rica) record a volatile evolution history that confirms previous melt inclusion studies, and provides additional information concerning the relative and absolute timing of subvolcanic magmatic events. Measurements of H, Cl, and F by secondary ion mass spectrometry reveal multiple populations of apatite in both 1723 and 1963 magmas. Assuming nominal apatite/melt partition coefficients allows us to compare the pattern of melt inclusions and apatites in ternary space, demonstrating the fidelity of the record preserved in apatite, and revealing a complex history of magma mixing with at least two components. The preservation of heterogeneous populations of apatite and of internally heterogeneous crystals requires short timescales (days to years) for these magmatic processes to occur.     

Sedimentary and faunal events revealed by a revised correlation of post‐glacial Hirnantian (Late Ordovician) strata in the Welsh Basin,UK

The discovery of a previously unrecognized unconformity and of new faunas in the type Llandovery area underpins a revised correlation of Hirnantian strata in mid Wales. This has revealed the sedimentary and faunal events which affected the Lower Palaeozoic Welsh Basin during the global rise in sea level that followed the end‐Ordovician glacial maximum and has allowed their interpretation in the context of local and global influences. In peri‐basinal shelfal settings the onset of post‐glacial deepening is recorded by an unfossiliferous, transgressive shoreface sequence (Cwm Clyd Sandstone and Garth House formations) which rests unconformably on Rawtheyan rocks, deformed during an episode of pre‐Hirnantian tectonism. In the deep water facies of the basin centre, this same sequence boundary is now recognized as the contact between fine‐grained, re‐sedimented mudstones and an underlying regressive sequence of turbidite sandstones and conglomerates; it is at a level lower than previously cited and calls into question the established lithostratigraphy. In younger Hirnantian strata, graptolites associated with the newly recognized Ystradwalter Member (Chwefri Formation) demonstrate that this distal shelf unit correlates with the persculptus graptolite‐bearing Mottled Mudstone Member of the basinal succession. Together these members record an important macrofaunal recolonization of the Welsh Basin and mark a key event in the post‐glacial transgression. Further deepening saw the establishment of a stratified water column and the imposition of anoxic bottom water conditions across the basin floor. These post‐glacial Hirnantian events are consistent with the re‐establishment of connections between a silled Welsh Basin and the open Iapetus Ocean. However, a comparison with other areas suggests that each event records a separate deepening episode within a pulsed glacio‐eustatic transgression, while also reflecting changes in post‐glacial climate and patterns of oceanic circulation and associated biotic flux. British Geological Survey © NERC 2009. All rights reserved.     

The distribution of carbon monoxide in the lower atmosphere of Venus

We have obtained spatially resolved near-infrared spectroscopy of the Venus nightside on 15 nights over three observing seasons. We use the depth of the CO absorption band at 2.3 μm to map the two-dimensional distribution of CO across both hemispheres. Radiative transfer models are used to relate the measured CO band depth to the volume mixing ratio of CO. The results confirm previous investigations in showing a general trend of increased CO abundances at around 60° latitude north and south as compared with the equatorial regions. Observations taken over a few nights generally show very similar CO distributions, but significant changes are apparent over longer periods. In past studies it has been assumed that the CO latitudinal variation occurs near 35 km altitude, at which K-band sensitivity to CO is greatest. By modeling the detailed spectrum of the excess CO at high latitudes we show that it occurs at altitudes around 45 km, much higher than has previously been assumed, and that there cannot be significant contribution from levels of 36 km or lower. We suggest that this is most likely due to downwelling of CO-rich gas from the upper atmosphere at these latitudes, with the CO being removed by around 40 km through chemical processes such as the reaction with SO₃.     

Extending the Capture Map Concept to Estimate Discrete and Risk-Based Streamflow Depletion Potential

A popular and contemporary use of numerical groundwater models is to estimate the discrete relation between groundwater extraction and surface-water/groundwater exchange. Previously, the concept of a “capture map” has been put forward as a means to effectively summarize this relation for decision-making consumption. While capture maps have enjoyed success in the environmental simulation industry, they are deterministic, ignoring uncertainty in the underlying model. Furthermore, capture maps are not typically calculated in a manner that facilitates analysis of varying combinations of extraction locations and/or reaches. That is, they are typically constructed with focus on a single reach or group of reaches. The former of these limitations is important for conveying risk to decision makers and stakeholders, while the latter is important for decision-making support related to surface-water management, where future foci may include reaches that were not the focus of the original capture analysis. Herein, we use the concept of a response matrix to generalize the theory of the capture-map approach to estimate spatially discrete streamflow depletion potential. We also use first-order, second-moment uncertainty estimation techniques with the concept of “risk shifting” to place capture maps and streamflow depletion potential in a stochastic, risk-based framework. Our approach is demonstrated for an integrated groundwater/surface-water model of the lower San Antonio River, Texas, USA.     

Constraining basin thermal history and petroleum generation using palaeoclimate data in the Piceance Basin,Colorado

Careful assessment of basin thermal history is critical to modelling petroleum generation in sedimentary basins. In this paper, we propose a novel approach to constraining basin thermal history using palaeoclimate temperature reconstructions and study its impact on estimating source rock maturation and hydrocarbon generation in a terrestrial sedimentary basin. We compile mean annual temperature (MAT) estimates from macroflora assemblage data to capture past surface temperature variation for the Piceance Basin, a high‐elevation, intermontane, sedimentary basin in Colorado, USA. We use macroflora assemblage data to constrain the temporal evolution of the upper thermal boundary condition and to capture the temperature change with basin uplift. We compare these results with the case where the upper thermal boundary condition is based solely upon a simplified latitudinal temperature estimate with no elevation effect. For illustrative purposes, 2 one‐dimensional (1‐D) basin models are constructed using these two different upper thermal boundary condition scenarios and additional geological and geochemical input data in order to investigate the impact of the upper thermal boundary condition on petroleum source rock maturation and kerogen transformation processes. The basin model predictions indicate that the source rock maturation is very sensitive to the upper thermal boundary condition for terrestrial basins with variable elevation histories. The models show substantial differences in source rock maturation histories and kerogen transformation ratio over geologic time. Vitrinite reflectance decreases by 0.21%Ro, source rock transformation ratio decreases 10.5% and hydrocarbon mass generation decreases by 16% using the macroflora assemblage data. In addition, we find that by using the macroflora assemblage data, the modelled depth profiles of vitrinite reflectance better matches present‐day measurements. These differences demonstrate the importance of constraining thermal boundary conditions, which can be addressed by palaeotemperature reconstructions from palaeoclimate and palaeo‐elevation data for many terrestrial basins. Although the palaeotemperature reconstruction compiled for this study is region specific, the approach presented here is generally applicable for other terrestrial basin settings, particularly basins which have undergone substantial subaerial elevation change over time.     

A GPS-Supported Gravity Survey in the Amazon of Ecuador

The growing importance and effectiveness of the Global Positioning System (GPS) in positioning geophysical surveys is now well established. In Ecuador's Amazon basin the use of GPS was the determining factor in the successful completion of a gravity survey carried out in 1994. Over 500 gravity stations, spaced from 1 to 5 km apart and positioned by GPS, were collected in the jungle environment in the country's southeast, where the preexisting data, of variable and unknown quality, were mostly confined to areas along the major rivers. This station densification, which includes the reoccupation of selected points, contributed to increase the gravity coverage; it also allowed us to evaluate and correct 1100 preexisting data points in the area, so that the different surveys can be merged. The composite data set is used for a geologic analysis of the study area, where previously unknown structures are discovered. ? 1999 John Wiley & Sons, Inc.     

Growth and mortality of juvenile winter flounder in two New England estuaries

A length-based model for calculating growth and mortality of juvenile winter flounder (Pseudopleuronectes americanus) populations has been developed. This model is based on work by Sullivan et al. (1990) and incorporates the von Bertalanffy growth equation, including stochasticity in growth, and a mortality rate that decreases exponentially with size. The length-based model was fit to observed size-frequency distributions, and model likelihood profiles were generated to produce 95% confidence intervals about parameter estimates. We analyzed size-frequency distributions of 3 to 15 cm juvenile winter flounder, collected with a 1-m beam trawl, at monthly intervals from June to October during 1993 and 1994. Growth rates were higher at a contaminated site, New Haven Harbor, than at a clean site, the Connecticut River estuary, however, the parameter estimates had overlapping 95% confidence intervals. Mortality rates were similar at the two sites.     

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.     

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.     

Evaluating Subsurface Parameterization to Simulate Hyporheic Exchange: The Steinlach River Test Site

Hyporheic exchange is the interaction of river water and groundwater, and is difficult to predict. One of the largest contributions to predictive uncertainty for hyporheic exchange has been attributed to the representation of heterogeneous subsurface properties. Our study evaluates the trade-offs between intrinsic (irreducible) and epistemic (reducible) model errors when choosing between homogeneous and highly complex subsurface parameter structures. We modeled the Steinlach River Test Site in Southwest Germany using a fully coupled surface water-groundwater model to simulate hyporheic exchange and to assess the predictive errors and uncertainties of transit time distributions. A highly parameterized model was built, treated as a “virtual reality” and used as a reference. We found that if the parameter structure is too simple, it will be limited by intrinsic model errors. By increasing subsurface complexity through the addition of zones or heterogeneity, we can begin to exchange intrinsic for epistemic errors. Thus, the appropriate level of detail to represent the subsurface depends on the acceptable range of intrinsic structural errors for the given modeling objectives and the available site data. We found that a zonated model is capable of reproducing the transit time distributions of a more detailed model, but only if the geological structures are known. An interpolated heterogeneous parameter field (cf. pilot points) showed the best trade-offs between the two errors, indicating fitness for practical applications. Parameter fields generated by multiple-point geostatistics (MPS) produce transit time distributions with the largest uncertainties, however, these are reducible by additional hydrogeological data, particularly flux measurements.