Morphology and variability of the Titan ringlet and Huygens ringlet edges

We present a forward modeling approach for determining, in part, the ring particle spatial distribution in the vicinity of sharp ring or ringlet edges. Synthetic edge occultation profiles are computed based on a two-parameter particle spatial distribution model. One parameter, h, characterizes the vertical extent of the ring and the other, δ, characterizes the radial scale over which the ring optical depth transitions from the background ring value to zero. We compare our synthetic occultation profiles to high resolution stellar occultation light curves observed by the Cassini Ultraviolet Imaging Spectrograph (UVIS) High Speed Photometer (HSP) for occultations by the Titan ringlet and Huygens ringlet edges.More than 100 stellar occultations of the Huygens ringlet and Titan ringlet edges were studied, comprising 343 independent occultation cuts of the edges of these two ringlets. In 237 of these profiles the measured light-curve was fit well with our two-parameter edge model. Of the remaining edge occultations, 69 contained structure that could only be fit with extremely large values of the ring-plane vertical thickness (h > 1 km) or by adopting a different model for the radial profile of the ring optical depth. An additional 37 could not be fit by our two-parameter model.Certain occultations at low ring-plane incidence angles as well as occultations nearly tangent to the ring edge allow the direct measurement of the radial scale over which the particle packing varies at the edge of the ringlet. In 24 occultations with these particular viewing geometries, we find a wide variation in the radial scale of the edge. We are able to constrain the vertical extent of the rings at the edge to less than ∼300 m in the 70% of the occultations with appropriate viewing geometry, however tighter constraints could not be placed on h due to the weaker sensitivity of the occultation profile to vertical thickness compared to its sensitivity to δ.Many occultations of a single edge could not be fit to a single value of δ, indicating large temporal or azimuthal variability, although the azimuthal variation in δ with respect to the longitudes of various moons in the system did not show any discernible pattern.     

Comparison of numerical methods for simulating strongly nonlinear and heterogeneous reactive transport problems—the MoMaS benchmark case

Although multicomponent reactive transport modeling is gaining wider application in various geoscience fields, it continues to present significant mathematical and computational challenges. There is a need to solve and compare the solutions to complex benchmark problems, using a variety of codes, because such intercomparisons can reveal promising numerical solution approaches and increase confidence in the application of reactive transport codes. In this contribution, the results and performance of five current reactive transport codes are compared for the 1D and 2D subproblems of the so-called easy test case of the MoMaS benchmark (Carrayrou et al., Comput Geosci, 2009, this issue). This benchmark presents a simple fictitious reactive transport problem that highlights the main numerical difficulties encountered in real reactive transport problems. As a group, the codes include iterative and noniterative operator splitting and global implicit solution approaches. The 1D easy advective and 1D easy diffusive scenarios were solved using all codes, and, in general, there was a good agreement, with solution discrepancies limited to regions with rapid concentration changes. Computational demands were typically consistent with what was expected for the various solution approaches. The differences between solutions given by the three codes solving the 2D problem are more important. The very high computing effort required by the 2D problem illustrates the importance of parallel computations. The most important outcome of the benchmark exercise is that all codes are able to generate comparable results for problems of significant complexity and computational difficulty.     

Upscaling facies models to preserve connectivity of designated facies

An upscaling algorithm has been developed that generates an irregular coarse grid that preserves flow connectivity by applying a rule-based upscaling algorithm to a fine-scale facies distribution. The algorithm is demonstrated using stochastically generated paleo-fluvial facies distributions. First, an irregular grid honoring the channel facies is created, followed by computation of effective anisotropic parameters for all coarse-grid cells. For the apparent layer-cake geometry of overbank deposits seen in outcrop, two local upscaling methods are compared: (1) the layered system approximation and (2) the mode. To assess upscaling performance, flow simulations for the original and upscaled grids are compared. The horizontal layered approximation (arithmetic mean) performs poorly, over-predicting lateral connectivity where even infrequent disconnection becomes important. Performance of the mode as an upscaling algorithm depends on the probability that a coarse-grid cell will be dominated by a single facies, and it performs surprisingly well because the upscaled grid-generation algorithm honors the channels, informing the upscaling process. Lastly, the irregular coarse grid was compared to a uniform coarse grid, showing superior performance with the irregular grid. The reduction in grid size achieved by irregular-grid generation will be a function of the geometrical complexity of the geologic objects to be honored.     

The geoarchaeological and paleoenvironmental context of Paleoindian sites in western Middle Park,Colorado, USA

Geoarchaeological investigations in western Middle Park provide important information for understanding the soil‐stratigraphic context of Paleoindian components, as well as the latest Quaternary environmental change and landscape evolution in a Southern Rocky Mountain intermontane basin. Paleoindian components are associated with the oldest two of four latest Quaternary stratigraphic units (1–4) recognized in co‐alluvial mantles (combined slopewash and colluvium) in uplands and in alluvial valley fills. Limited data suggest accumulation of unit 1 as early as ∼12,500 ¹⁴C yr B.P. in alluvial valleys and by at least ∼11,000 ¹⁴C yr B.P. in uplands was followed by brief stability and soil formation. A relatively widespread disconformity marks earliest Holocene erosion and substantial removal of latest Pleistocene deposits in upland and alluvial settings followed by unit 2 deposition ∼10,000–9000 ¹⁴C yr B.P., perhaps signaling the abrupt onset of an intensified summer monsoon. In situ Paleoindian components in uplands are found in a moderately developed buried soil (the Kremmling soil) formed in units 1 and 2 in thin (≤1m) hillslope co‐alluvial mantles. The Kremmling soil reflects geomorphic stability in upland and alluvial settings ∼9000–4500 ¹⁴C yr BP, and represents a buried landscape with the potential to contain additional Paleoindian components, although elsewhere in western Middle Park Early Archaic components are documented in morphologically similar soils. Kremmling soil morphology, the relative abundance of charcoal in unit 2 relative to younger units, and charcoal morphology indicate the expansion of forest cover, including Pinus, and grass cover during the early and middle Holocene, suggesting conditions moister than present. © 2010 Wiley Periodicals, Inc.     

Extracting olivine (Fo–Fa) compositions from Raman spectral peak positions

The dominant feature of the olivine Raman spectrum is a doublet that occurs in the spectral region of 815–825 cm⁻¹ (DB1) and 838–857 cm⁻¹ (DB2). These features arise from coupled symmetric and asymmetric stretching vibrational modes of the constituent SiO₄ tetrahedra. The frequencies of both peaks show monotonic shifts following cation substitution between forsterite and fayalite. We present a calibration for extracting olivine Fo contents (Fo = Mg/(Mg + Fe) molar ratio; Fo_0–100) from the peak positions of this doublet, permitting estimates of chemical composition from Raman spectra (acquired in the laboratory or field) as well as providing information on crystal structure (distinction of polymorphs). Eight samples spanning the compositional range from forsterite to fayalite were used to develop the calibration equations for the DB1 and DB2 peaks individually and together. The data indicate that the DB1 peak is more reliable for calculating the compositions of Fe-rich olivine but that the DB2 peak is better for magnesian compositions. The two-peak calibration overcomes the limitations of the single-peak calibrations and is capable of calculating olivine compositions to within ±10 Fo units.     

What Can Groundwater Monitoring Tell Us About Gas Migration? A Numerical Modeling Study

Groundwater monitoring to measure a variety of indicator parameters including dissolved gas concentrations, total dissolved gas pressure (TDGP), and redox indicators is commonly used to evaluate the impacts of gas migration (GM) from energy development in shallow aquifer systems. However, these parameters can be challenging to interpret due to complex free-phase gas source architecture, multicomponent partitioning, and biogeochemical reactions. A series of numerical simulations using a gas flow model and a reactive transport model were conducted to delineate the anticipated evolution of indicator parameters following GM in an aquifer under a variety of physical and biogeochemical conditions. The simulations illustrate how multicomponent mass transfer processes and biogeochemical reactions create unexpected spatial and temporal variations in several analytes. The results indicate that care must be taken when interpreting measured indicator parameters including dissolved hydrocarbon concentrations and TDGP, as the presence of dissolved gases in background groundwater and biogeochemical processes can cause potentially misleading conclusions about the impact of GM. Based on the consideration of multicomponent gas partitioning in this study, it is suggested that dissolved background gases such as N₂ and Ar can provide valuable insights on the presence, longevity and fate of free-phase natural gas in aquifer systems. Overall, these results contribute to developing a better understanding of indicators for GM in groundwater, which will aid the planning of future monitoring networks and subsequent data interpretation.