Approaches to the simulation of unconfined flow and perched groundwater flow in MODFLOW

Various approaches have been proposed to manage the nonlinearities associated with the unconfined flow equation and to simulate perched groundwater conditions using the MODFLOW family of codes. The approaches comprise a variety of numerical techniques to prevent dry cells from becoming inactive and to achieve a stable solution focused on formulations of the unconfined, partially-saturated, groundwater flow equation. Keeping dry cells active avoids a discontinuous head solution which in turn improves the effectiveness of parameter estimation software that relies on continuous derivatives. Most approaches implement an upstream weighting of intercell conductance and Newton-Raphson linearization to obtain robust convergence. In this study, several published approaches were implemented in a stepwise manner into MODFLOW for comparative analysis. First, a comparative analysis of the methods is presented using synthetic examples that create convergence issues or difficulty in handling perched conditions with the more common dry-cell simulation capabilities of MODFLOW. Next, a field-scale three-dimensional simulation is presented to examine the stability and performance of the discussed approaches in larger, practical, simulation settings.     

Compositional gradients surrounding spherulites in obsidian and their relationship to spherulite growth and lava cooling

Spherical masses of crystal fibers (spherulites) crystalize from rhyolitic melt/glass mainly in response to significant undercooling while lava cools. Spherulite growth should induce compositional gradients in the surrounding glass from expulsion of incompatible constituents and diffusion of those constituents away from the spherulite. Finite-difference numerical modeling of one-dimensional diffusion, in which diffusivities are allowed to vary with temperature, is used to investigate how compositional gradients reflect spherulite growth and lava cooling. Overall, three forms of gradients are identified. Elements that diffuse quickly are expelled from the spherulite but then migrate away too quickly to become enriched at the boundary of the spherulite. Elements that diffuse slowly are trapped within the growing spherulite. Between those endmembers are elements that are not trapped, yet diffuse slow enough that they become enriched at the contact. Their slow diffusion away then elevates their concentrations in the surrounding glass. How enriched those elements are at the spherulite-matrix interface and how far their enrichments extend outwards into the glass reflect how spherulites grow and thermal conditions during growth. Concentrations of H₂O, Rb, F, Li, Cl, Na, K, Sr, Cs, Ba, and Be were measured in and around spherulites in obsidian from a 4.7?±?1?km³ rhyolite lava dome erupted from Tequila volcano, Mexico. Measurable concentration gradients are found for H₂O, Rb, and F. Attributes of those gradients and the behaviors of the other elements are in accord with their experimentally constrained diffusivities. Spherulites appear to have grown following radial, rather than volumetric, growth. The observed gradients (and lack of others) are more consistent with growth mainly below the glass transition, which would necessitate the dome cooling at ca. 10^?5 to 10^?7?°C?s^?1. Such slow cooling is consistent with the relatively large volume of the dome.     

Simulating bubble number density of rhyolitic pumices from Plinian eruptions: constraints from fast decompression experiments

Decompression experiments of a crystal-free rhyolitic liquid with ≈ 6.6 wt. % H₂O were carried out at a pressure range from 250 MPa to 30–75 MPa in order to characterize effects of magma ascent rate and temperature on bubble nucleation kinetics, especially on the bubble number density (BND, the number of bubbles produced per unit volume of liquid). A first series of experiments at 800°C and fast decompression rates (10–90 MPa/s) produced huge BNDs (≈ 2 × 10¹⁴ m⁻³ at 10 MPa/s ; ≈ 2 × 10¹⁵ m⁻³ at 90 MPa/s), comparable to those in natural silicic pumices from Plinian eruptions (10¹⁵–10¹⁶ m⁻³). A second series of experiments at 700°C and 1 MPa/s produced BNDs (≈ 9×10¹² m⁻³) close to those observed at 800°C and 1 MPa/s (≈ 6 × 10¹² m⁻³), showing that temperature has an insignificant effect on BNDs at a given decompression rate. Our study strengthens the theory that the BNDs are good markers of the decompression rate of magmas in volcanic conduits, irrespective of temperature. Huge number densities of small bubbles in natural silicic pumices from Plinian eruptions imply that a major nucleation event occurs just below the fragmentation level, at which the decompression rate of ascending magmas is a maximum (≥ 1 MPa/s).     

Factors Controlling the Regional Distribution of Vanadium in Groundwater

Although the ingestion of vanadium (V) in drinking water may have possible adverse health effects, there have been relatively few studies of V in groundwater. Given the importance of groundwater as a source of drinking water in many areas of the world, this study examines the potential sources and geochemical processes that control the distribution of V in groundwater on a regional scale. Potential sources of V to groundwater include dissolution of V rich rocks, and waste streams from industrial processes. Geochemical processes such as adsorption/desorption, precipitation/dissolution, and chemical transformations control V concentrations in groundwater. Based on thermodynamic data and laboratory studies, V concentrations are expected to be highest in samples collected from oxic and alkaline groundwater. However, the extent to which thermodynamic data and laboratory results apply to the actual distribution of V in groundwater is not well understood. More than 8400 groundwater samples collected in California were used in this study. Of these samples, high (≥50 µg/L) and moderate (25 to 49 µg/L) V concentrations were most frequently detected in regions where both source rock and favorable geochemical conditions occurred. The distribution of V concentrations in groundwater samples suggests that significant sources of V are mafic and andesitic rock. Anthropogenic activities do not appear to be a significant contributor of V to groundwater in this study. High V concentrations in groundwater samples analyzed in this study were almost always associated with oxic and alkaline groundwater conditions, which is consistent with predictions based on thermodynamic data.     

Iron in the southeastern Bering Sea: Elevated leachable particulate Fe in shelf bottom waters as an important source for surface waters

Surface transects and vertical profiles of total and leachable particulate Fe, Mn, Al and P, along with dissolved and soluble Fe were obtained during August 2003 in the southeastern Bering Sea. High concentrations of leachable particulate Fe were observed in the bottom waters over the Bering Sea shelf with an unusually high percentage of the suspended particulate Fe being leachable. Leachable particulate Fe averaged 81% of total particulate Fe, and existed at elevated concentrations that averaged 23 times greater than dissolved Fe in the isolated cool pool waters over the mid shelf where substantial influence of sedimentary denitrification was apparent. The elevated leachable particulate Fe is suggested to be a result of sedimentary Fe reduction in surficial sediments, diffusion of Fe(II) from the sediments to the bottom waters, and subsequent oxidation and precipitation of reduced Fe in the overlying bottom waters. Eddies and meanders of the Bering Slope Current can mix this Fe-rich water into the Green Belt at the outer shelf-break front. Elevated levels of leachable particulate Fe were observed in surface waters near the Pribilof Islands where enhanced vertical mixing exists. Storm events and/or cooling during fall/winter with the resultant destruction of the thermally stratified two-layer system can also mix the subsurface water into surface waters where the elevated leachable particulate Fe is a substantial source of biologically available Fe.     

Cartoblography: Experiments in Using and Organising the Spatial Context of Micro‐blogging

Blogs, micro‐blogs and online forums underpin a more interconnected world. People communicate ever more and are increasingly keen to explain and illustrate their lives; showing where they are and what they are doing. Desktop, online and mobile mapping landscapes have never been as rich or diverse yet this challenges cartography to adapt and remain relevant in the modern mapping world. We explore the spatial expression and potential value of micro‐blogging and Twitter as a social networking tool. Examples of “twitter maps” are reviewed that leverage the Twitter API and online map services to locate some component of the “tweet”. Scope, function and design are illustrated through development of two proof‐of‐concept map mashups that support collaborative real‐time mapping and the organisation and display of information for mass user events. Through the experiments in using and organising data in this way we demonstrate the value of “cartoblography”– a framework for mapping the spatial context of micro‐blogging.     

Insight into sediment transport processes on saline rangeland hillslopes using three‐dimensional soil microtopography changes

In arid and semi‐arid rangeland environments, an accurate understanding of runoff generation and sediment transport processes is key to developing effective management actions and addressing ecosystem response to changes. Yet, many primary processes (namely sheet and splash and concentrated flow erosion, as well as deposition) are still poorly understood due to a historic lack of measurement techniques capable of parsing total soil loss into these primary processes. Current knowledge gaps can be addressed by combining traditional erosion and runoff measurement techniques with image‐based three‐dimensional (3D) soil surface reconstructions. In this study, data (hydrology, erosion and high‐resolution surface microtopography changes) from rainfall simulation experiments on 24 plots in saline rangelands communities of the Upper Colorado River Basin were used to improve understanding on various sediment transport processes. A series of surface change metrics were developed to quantify and characterize various erosion and transport processes (e.g. plot‐wide versus concentrated flow detachment and deposition) and were related to hydrology and biotic and abiotic land surface characteristics. In general, erosivity controlled detachment and transport processes while factors modulating surface roughness such as vegetation controlled deposition. The extent of the channel network was a positive function of slope, discharge and vegetation. Vegetation may deflect runoff in many flow paths but promoted deposition. From a management perspective, this study suggests that effective runoff soil and salt load reduction strategies should aim to promote deposition of transported sediments rather than reducing detachment which might not be feasible in these resource‐limited environments. Copyright © 2016 John Wiley & Sons, Ltd.     

F,Cl and S input via serpentinite in subduction zones: implications for the nature of the fluid released at depth

The abundances of F, Cl and S in arc magmas are systematically higher than in other mantle‐derived magmas, suggesting that these elements are added from the slab along with H₂O. We present ion probe microanalyses of F, Cl and S in serpentine minerals that represent the P–T evolution of the oceanic lithosphere, from its serpentinization at the ridge, to its dehydration at around 100 km depth during subduction. F, Cl and S are incorporated early into serpentine during its formation at mid‐ocean ridges, and serpentinized lithosphere then carries these elements to subduction zones. More than 50% of the F, Cl and S are removed from serpentine during the prograde metamorphic lizardite/antigorite transition. Due to the low solubility of F in water, and to the low amount of water released during this phase transition, the fluids mobilizing these elements must be dominated by SO_X rather than H₂O.     

Isotopic fractionation associated with biosynthesis of fatty acids by a marine bacterium under oxic and anoxic conditions

Shewanella putrefaciens (Strain MR-4), a gram negative facultative marine bacterium, was grown to stationary phase under both aerobic and anaerobic conditions using lactate as the sole carbon source. Aerobically-produced cells were slightly enriched in ¹³C (+1.5‰) relative to the lactate carbon source, whereas those from anaerobic growth were depleted in ¹³C (−2.2‰). The distribution of fatty acids produced under aerobic conditions was similar to that resulting from anaerobic growth, being dominated by C_16:1 ω7 and C_16:0 fatty acids with a lesser amount of the C_18:1 ω7 component. Low concentrations of saturated even numbered normal fatty acids in the C₁₄ to C₁₈ range, and iso-C_15:0 were synthesized under both conditions. Fatty acids from anaerobic cultures (average δ¹³C=−37.8‰) were considerably depleted in ¹³C relative to their aerobically-synthesized counterparts (−28.8‰). The distinct differences in isotopic composition of both whole cells and individual fatty acid components result from differences in assimilation pathways. Under aerobic conditions, the primary route of assimilation involves the pyruvate dehydrogenase enzyme complex which produces acetyl-CoA, the precursor to lipid synthesis. In contrast, under anaerobic conditions formate, and not acetate, is the central intermediate in carbon assimilation with the precursors to fatty acid synthesis being produced via the serine pathway. Anaerobically-produced bacterial fatty acids were depleted by up to 12‰ relative to the carbon source. Therefore, detection of isotopically depleted fatty acids in sediments may be falsely attributed to a terrestrial origin, when in fact they are the result of bacterial resynthesis.