Anisotropic dual-continuum representations for multiscale poroelastic materials: Development and numerical modelling

Dual-continuum (DC) models can be tractable alternatives to explicit approaches for the numerical modelling of multiscale materials with multiphysics behaviours. This work concerns the conceptual and numerical modelling of poroelastically coupled dual-scale materials such as naturally fractured rock. Apart from a few exceptions, previous poroelastic DC models have assumed isotropy of the constituents and the dual-material. Additionally, it is common to assume that only one continuum has intrinsic stiffness properties. Finally, little has been done into validating whether the DC paradigm can capture the global poroelastic behaviours of explicit numerical representations at the DC modelling scale. We address the aforementioned knowledge gaps in two steps. First, we utilise a homogenisation approach based on Levin's theorem to develop a previously derived anisotropic poroelastic constitutive model. Our development incorporates anisotropic intrinsic stiffness properties of both continua. This addition is in analogy to anisotropic fractured rock masses with stiff fractures. Second, we perform numerical modelling to test the DC model against fine-scale explicit equivalents. In doing, we present our hybrid numerical framework, as well as the conditions required for interpretation of the numerical results. The tests themselves progress from materials with isotropic to anisotropic mechanical and flow properties. The fine-scale simulations show that anisotropy can have noticeable effects on deformation and flow behaviour. However, our numerical experiments show that the DC approach can capture the global poroelastic behaviours of both isotropic and anisotropic fine-scale representations.     

Making data useful for modelers to understand complex Earth systems

The Arctic is changing rapidly with dramatic local and global effect. To understand that change requires understanding the Arctic as a system. Models of different processes and at various scales are necessary tools for analyzing and understanding the Earth system. Models are extremely diverse, yet they all require quality data. Through a series of case studies, augmented with ethnographic observation around the International Polar Year, this work examines how modelers assess, acquire, and prepare data for their models. By comparing specific case studies, common themes emerge that can be compared against broader observation. These themes, in turn, suggest data management techniques or requirements for data systems to improve access and use by modelers and generally improve understanding of the Arctic system. The study has an Arctic focus because of the rapid changes occurring in the Arctic, but the approach and results should apply generally to Earth system science. This case study based approach has proven to be a useful method for teasing out both general and specific data needs for different models. An overarching lesson is that greater short-term benefit to modelers and significant gains in efficiency can be achieved by improving the formats, convention, and consistency of the data rather than improved interfaces and analysis tools. A “data-first” philosophy can improve the data systems that support the overall interdisciplinary, integrative science necessary to understand the complex Earth system.     

Chasing a 'loose and baggy monster': almshouses and the geography of charity

This paper goes some way towards redressing the lack of geographical literature on charity through exploring the geography of the British domestic charitable sector. The size and geography of the third sector is outlined, followed by an analysis of how almshouses can be understood as inherently geographical and deeply embedded in local social networks of inclusion as well as exclusion.     

Structural improvement of age-hardened gypsum-treated bauxite red mud waste using readily decomposable phyto-organics

Several years (7) after incidental lithification by added gypsum, the effect of crushing and phyto-organics on some physical properties of red mud wastes (RMW) was examined under laboratory conditions. It was hypothesized that mechanically breaking up the hardened masses would increase their permeability. However, on the contrary, crushing these rock-hard masses into <2-mm particles did not increase permeability. A phyto-organic treatment on the other hand greatly improved permeability of the crushed entities such that infiltration rates increased from 0 to 10 cm/h. After 24 weeks of incubation, the structure of the newly formed phyto-organically treated gypsum-hardened RMW aggregates was also substantially improved over that of the controls, as a >3-fold increase in resistance to disintegration during wet sieving was exhibited by the phyto-organically treated samples. The increase in water infiltration by combining crushing with phyto-organic applications was ascribed to a decrease in particle dispersion, creating more stable clay clusters, which in turn stabilized channels and pores. Conversely, it was concluded that the singular crushing of gypsum-hardened bauxite waste in the waste dumps is potentially wasteful and ineffective.     

The composition and crystallinity of the near-surface regions of weathered alkali feldspars

Our ability to identify thin non-stoichiometric and amorphous layers beneath mineral surfaces has been tested by undertaking X-ray photoelectron spectroscopy (XPS) and transmission electron microscopy (TEM) work on alkali feldspars from pH 1 dissolution experiments. The outcomes of this work were used to help interpret XPS and TEM results from alkali feldspars weathered for <10,000 years in soils overlying the Shap Granite (north-west England). The chemistry of effluent solutions indicates that silica-rich layers a few nanometers in thickness formed during the pH 1 experiments. These layers can be successfully identified by XPS and have lower Al/Si, Na/Si, K/Si and Ca/Si values than the outermost ∼9 nm of unweathered controls. Development of Al-Si non-stoichiometry is coupled with loss of crystal structure to produce amorphous layers that are identifiable by TEM where >∼2.5 nm thick, whereas the crystallinity of albite is retained despite leaching of Na to depths of tens to hundreds on nanometers. Integration of XPS data over the outermost 6-9 nm of naturally weathered Shap feldspars shows that they have stoichiometric Al/Si and K/Si ratios, which is consistent with findings of previous TEM work on the same material that they lack amorphous layers. There is some XPS evidence for loss of K from the outermost couple of nanometers of Shap orthoclase, and the possibility of leaching of Na from albite to greater depths cannot be excluded using the XPS or TEM results. This study demonstrates that the leached layer model, as formulated from laboratory experiments, is inapplicable to the weathering of alkali feldspars within acidic soils, which is an essentially stoichiometric reaction.     

A Numerical Simulation of Residual Circulation in Tampa Bay. Part II: Lagrangian Residence Time

Lagrangian retention and flushing are examined by advecting neutrally buoyant point particles within a circulation field generated by a numerical ocean model of Tampa Bay. Large temporal variations in Lagrangian residence time are found under realistic changes in boundary conditions. Two 90-day time periods are examined. The first (P1) is characterized by low freshwater inflow and weak baroclinic circulation. The second (P2) has high freshwater inflow and strong baroclinic circulation. At the beginning of both time periods, 686,400 particles are released uniformly throughout the bay. Issues relating to particle distribution and flushing are examined at three different spatial scales: (1) at the scale of the entire bay, (2) the four major regions within the bay, and (3) at the scale of individual model grid cells. Two simple theoretical models for the particle number over time, N(t), are fit to the particle counts from the ocean model. The theoretical models are shown to represent N(t) reasonably well when considering the entire bay, allowing for straightforward calculation of baywide residence times: 156 days for P1 and 36 days for P2. However, the accuracy of these simple models decreases with decreasing spatial scale. This is likely due to the fact that particles may exit, reenter, or redistribute from one region to another in any sequence. The smaller the domain under consideration, the more this exchange process dominates. Therefore, definitions of residence time need to be modified for “non-local” situations. After choosing a reasonable definition, and removal of the tidal and synoptic signals, the residence times at each grid cell in P1 is found to vary spatially from a few days to 90 days, the limit of the calculation, with an average residence time of 53 days. For P2, the overall spatial pattern is more homogeneous, and the residence times have an average value of 26 days.     

ENSO impacts on salinity in Tampa Bay, Florida

Estuarine salinity distributions reflect a dynamic balance between the processes that control estuarine circulation. At seasonal and longer time scales, freshwater inputs into estuaries represent the primary control on salinity distribution and estuarine circulation. El Niño-Southern Oscillation (ENSO) conditions influence seasonal rainfall and stream discharge patterns in the Tampa Bay, Florida region. The resulting variability in freshwater input to Tampa Bay influences its seasonal salinity distribution. During El Niño events, ENSO sea surface temperature anomalies (SSTAs) are significantly and inversely correlated with salinity in the bay during winter and spring. These patterns reflect the elevated rainfall over the drainage basin and the resulting elevated stream discharge and runoff, which depress salinity levels. Spatially, the correlations are strongest at the head of the bay, especially in bay sections with long residence times. During La Niña conditions, significant inverse correlations between ENSO SSTAs and salinity occur during spring. Dry conditions and depressed stream discharge characterize La Niña winters and springs, and the higher salinity levels during La Niña springs reflect the lower freshwater input levels.