Multidimensional upstream weighting for multiphase transport on general grids

The governing equations for multiphase flow in porous media have a mixed character, with both nearly elliptic and nearly hyperbolic variables. The flux for each phase can be decomposed into two parts: (1) a geometry- and rock-dependent term that resembles a single-phase flux; and (2) a mobility term representing fluid properties and rock–fluid interactions. The first term is commonly discretized by two- or multipoint flux approximations (TPFA and MPFA, respectively). The mobility is usually treated with single-point upstream weighting (SPU), also known as dimensional or donor cell upstream weighting. It is well known that when simulating processes with adverse mobility ratios, SPU suffers from grid orientation effects. An important example of this, which will be considered in this work, is the displacement of a heavy oil by water. For these adverse mobility ratio flows, the governing equations are unstable at the modeling scale, rendering a challenging numerical problem. These challenges must be addressed in order to avoid systematic biasing of simulation results. In this work, we present a framework for multidimensional upstream weighting for multiphase flow with buoyancy on general two-dimensional grids. The methodology is based on a dual grid, and the resulting transport methods are provably monotone. The multidimensional transport methods are coupled with MPFA methods to solve the pressure equation. Both explicit and fully implicit approaches are considered for time integration of the transport equations. The results show considerable reduction of grid orientation effects compared to SPU, and the explicit multidimensional approach allows larger time steps. For the implicit method, the total number of non-linear iterations is also reduced when multidimensional upstream weighting is used.     

Estimates of groundwater discharge to a coastal wetland using multiple techniques: Taylor Slough, Everglades National Park, USA

Quantifying water exchange between a coastal wetland and the underlying groundwater is important for closing water, energy and chemical budgets. The coastal wetlands of the Florida Everglades (USA) are at the forefront of a large hydrologic restoration project, and understanding of groundwater/surface-water interactions is needed to comprehend the effects of the project. Four independent techniques were used to identify water exchange at varying spatial and temporal scales in Taylor Slough, Everglades National Park. The techniques included a water-budget study and measurements of hydraulic head gradients, geochemical tracers, and temperature. During the 18-month study, the four methods converged as to the timing of groundwater discharge, typically between June and September, contemporaneous with the wet season and increasing surface-water levels. These results were unexpected, as groundwater discharge was predicted to be greatest when surface-water levels were low, typically during the dry season. Either a time lag of 1?C5?months in the response of groundwater discharge to low surface-water levels or precipitation-induced groundwater discharge may explain the results. Groundwater discharge was a significant contributor (27?%) to the surface water in Taylor Slough with greater rates of discharge observed towards the coastline in response to seawater intrusion.     

Implications of recent sea level rise science for low-elevation areas in coastal cities of the conterminous U.S.A.

Recently published work estimates that global sea level rise (SLR) approaching or exceeding 1 m by 2100 is plausible, thus significantly updating projections by the Fourth Assessment of the Intergovernmental Panel on Climate Change. Furthermore, global greenhouse gas (GHG) emissions over the 21st century will not only influence SLR in the next ??90 years, but will also commit Earth to several meters of additional SLR over subsequent centuries. In this context of worsening prospects for substantial SLR, we apply a new geospatial dataset to calculate low-elevation areas in coastal cities of the conterminous U.S.A. potentially impacted by SLR in this and following centuries. In total, 20 municipalities with populations greater than 300,000 and 160 municipalities with populations between 50,000 and 300,000 have land area with elevations at or below 6 m and connectivity to the sea, as based on the 1 arc-second National Elevation Dataset. On average, approximately 9% of the area in these coastal municipalities lies at or below 1 m. This figure rises to 36% when considering area at or below 6 m. Areal percentages of municipalities with elevations at or below 1?C6 m are greater than the national average along the Gulf and southern Atlantic coasts. In contrast to the national and international dimensions of and associated efforts to curb GHG emissions, our comparison of low-elevation areas in coastal cities of the conterminous U.S.A. clearly shows that SLR will potentially have very local, and disproportionate, impacts.     

Thermal history and origin of the Tanzanian Craton from Pb isotope thermochronology of feldspars from lower crustal xenoliths

Common and radiogenic Pb isotopic compositions of plagioclase and anti-perthitic feldspars from granulite-facies lower crustal xenoliths from the Labait Volcano on the eastern margin of the Tanzanian Craton have been measured via laser ablation MC-ICP-MS. Common Pb in plagioclase and a single stage Pb evolution model indicate that the lower crust of the Tanzanian Craton was extracted from mantle having a ²³⁸U/²⁰⁴Pb of 8.1 ± 0.3 and a ²³²Th/²³⁸U of 4.3 ± 0.1 at 2.71 ± 0.09 Ga (all uncertainties are 2σ). Since 2.4 Ga, some orthoclase domains within anti-perthites have evolved with a maximum ²³⁸U/²⁰⁴Pb of 6 and ²³²Th/²³⁸U of 4.3. The spread in Pb isotopic composition in the anti-perthitic feldspars yields single crystal Pb–Pb isochrons of ～ 2.4 Ga, within uncertainty of U–Pb zircon ages from the same sample suite. The Pb isotopic heterogeneities imply that these granulites resided at temperatures < 600 °C in the lower crust of the Tanzanian Craton from ca. 2.4 Ga to the present. In concert with the chemistry of surface samples, mantle xenoliths, and lower crustal xenoliths, our data imply that the cratonic lithosphere in Tanzania formed ca. ～ 2.7 Ga, in a convergent margin setting, and has remained undisturbed since 2.7 Ga.     

Thermodynamic characterization of boltwoodite and uranophane: Enthalpy of formation and aqueous solubility

Boltwoodite and uranophane are uranyl silicates common in oxidized zones of uranium ore deposits. An understanding of processes that impact uranium transport in the environment, especially pertaining to the distribution of uranium between solid phases and aqueous solutions, ultimately requires determination of thermodynamic parameters for such crystalline materials. We measured formation enthalpies of synthetic boltwoodites, K(UO₂)(HSiO₄)·H₂O and Na(UO₂)(HSiO₄)·H₂O, and uranophane, Ca(UO₂)₂(HSiO₄)₂·5H₂O, by high temperature oxide melt solution calorimetry. We also studied the aqueous solubility of these phases from both saturated and undersaturated conditions at a variety of pH. The combined data permit the determination of standard enthalpies, entropies and Gibbs free energies of formation for each phase and analysis of its potential geological impact from a thermodynamic point of view.     

Determination of arsenic concentration and distribution in the Floridan Aquifer System

Here we report on the different sampling strategies for almost seven years of sampling rocks/sediments for the determination of As within the Intermediate Aquifer System (IAS) and upper Floridan Aquifer System (FAS), a very large and productive limestone aquifer spanning from Georgia into Florida. In the FAS, As contamination has become a recurring problem during aquifer storage and recovery (ASR), particularly in central and south Florida.To investigate these phenomena, samples from solid drill cores and rock cuttings were collected from the Hawthorn Group, Suwannee Limestone, Ocala Limestone and Avon Park Formation. Samples were taken along drill cores and rock cuttings (referred to as ‘interval’ samples) or from particular drill core sections and rock cuttings (referred to as ‘targeted’ samples) likely to contain elevated concentrations of As as indicated by the presence of pyrite, hydrous ferric oxide, organic matter, clay minerals, fracture surfaces, and high permeable (moldic) zones.Arsenic was present in all of the stratigraphic units at low concentrations, close to the global average for As in limestone of 2.6 mg/kg. The highest As concentration was 69 mg/kg. In all units, however, the average bulk As concentration in the targeted samples was substantially higher than that in the interval samples. Based on direct spot measurements by electron microprobe and indirect calculations, pyrite was identified as the main source of As in the FAS. Concentrations in pyrite ranged from less than 100 mg/kg to more than 11,000 mg/kg. Because pyrite is heterogeneously distributed, both vertically and horizontally in the sampled stratigraphic units, the same was observed for the distribution of As. However, As concentrations generally decreased with depth, i.e., highest As values in the Hawthorn Group and lowest As values in the Ocala Limestone and Avon Park Formation. Compared to pyrite, other trace minerals contained much less As.The average As concentrations of the two types of sample media (solid cores and rock cuttings) were quite similar. These results indicate that if simply the average bulk rock As concentration of a geologic unit is the desired outcome of an investigation, either interval or targeted sampling of rock cuttings, seems to be sufficient. This is particularly important when time and money are a factor. This approach could work equally well for any other trace element. Structural sedimentary information, such as fractures, etc., is likely lost, however, when sampling rock cuttings. Thus, if this information is required, solid core samples need to be collected by hollow core diamond drilling.     

The domestic divide: Access to the Internet in New Zealand

We use unit records of the 2006 census to show that access to the Internet in the home varies geographically in New Zealand primarily as a result of demographic and socio‐economic differences among individuals. Of particular significance is the much lower rates of domestic access experienced by Māori and Pacific Island individuals even after controlling for differences in their age, gender, education, income, occupation and settlement type. While differences in Internet access by ethnicity has been noted before, it is the magnitude and persistence of this difference in New Zealand after controlling for correlated factors that renders this study unique. Our results have important implications in an education environment increasingly reliant on Web access, but they also raise questions about the extent of access to the Internet outside the home.     

An empirical model of subcritical bedform migration

The ability to predict bedform migration in rivers is critical for estimating bed material load, yet there is no relation for predicting bedform migration (downstream translation) that covers the full range of conditions under which subcritical bedforms develop. Here, the relation between bedform migration rates and transport stage is explored using a field and several flume data sets. Transport stage is defined as the non‐dimensional Shields stress divided by its value at the threshold for sediment entrainment. Statistically significant positive correlations between both ripple and dune migration rates and transport stage are found. Stratification of the data by the flow depth to grain‐size ratio improved the amount of variability in migration rates that was explained by transport stage to ca 70%. As transport stage increases for a given depth to grain‐size ratio, migration rates increase. For a given transport stage, the migration rate increases as the flow depth to grain‐size ratio gets smaller. In coarser sediment, bedforms move faster than in finer sediment at the same transport stage. Normalization of dune migration rates by the settling velocity of bed sediment partially collapses the data. Given the large amount of variability that arises from combining data sets from different sources, using different equipment, the partial collapse is remarkable and warrants further testing in the laboratory and field.     

It starts at home? Climate policies targeting household consumption and behavioral decisions are key to low-carbon futures

Through their consumption behavior, households are responsible for 72% of global greenhouse gas emissions. Thus, they are key actors in reaching the 1.5 °C goal under the Paris Agreement. However, the possible contribution and position of households in climate policies is neither well understood, nor do households receive sufficiently high priority in current climate policy strategies. This paper investigates how behavioral change can achieve a substantial reduction in greenhouse gas emissions in European high-income countries. It uses theoretical thinking and some core results from the HOPE research project, which investigated household preferences for reducing emissions in four European cities in France, Germany, Norway and Sweden. The paper makes five major points: First, car and plane mobility, meat and dairy consumption, as well as heating are the most dominant components of household footprints. Second, household living situations (demographics, size of home) greatly influence the household potential to reduce their footprint, even more than country or city location. Third, household decisions can be sequential and temporally dynamic, shifting through different phases such as childhood, adulthood, and illness. Fourth, short term voluntary efforts will not be sufficient by themselves to reach the drastic reductions needed to achieve the 1.5 °C goal; instead, households need a regulatory framework supporting their behavioral changes. Fifth, there is a mismatch between the roles and responsibilities conveyed by current climate policies and household perceptions of responsibility. We then conclude with further recommendations for research and policy.