Interior boundary-aligned unstructured grid generation and cell-centered versus vertex-centered CVD-MPFA performance

Grid generation for reservoir simulation must honor classical key constraints and be boundary aligned such that control-volume boundaries are aligned with geological features such as layers, shale barriers, fractures, faults, pinch-outs, and multilateral wells. An unstructured grid generation procedure is proposed that automates control-volume and/or control point boundary alignment and yields a PEBI-mesh both with respect to primal and dual (essentially PEBI) cells. In order to honor geological features in the primal configuration, we introduce the idea of protection circles, and to generate a dual-cell feature based grid, we construct halos around key geological features. The grids generated are employed to study comparative performance of cell-centred versus cell-vertex control-volume distributed multi-point flux approximation (CVD-MPFA) finite-volume formulations using equivalent degrees of freedom. The formulation of CVD-MPFA schemes in cell-centred and cell-vertex modes is analogous and requires switching control volume from primal to dual or vice versa together with appropriate data structures and boundary conditions. The relative benefits of both types of approximation, i.e., cell-centred versus vertex-centred, are made clear in terms of flow resolution and degrees of freedom required.     

Contrasting sources of Late Paleozoic rhyolite magma in the Polish Lowlands: evidence from U–Pb ages and Hf and O isotope composition in zircon

International Journal of Earth Sciences - The Polish Lowlands, located southwest of the Teisseyre–Tornquist Zone, within Trans-European Suture Zone, were affected by bimodal, but dominantly...     

Do Macroalgal Mats Impact Microphytobenthos on Mudflats? Evidence from a Meta-Analysis,Comparative Survey,and Large-Scale Manipulation

Eutrophication and species introductions have resulted in increased macroalgal biomass in coastal ecosystems around the globe. Macroalgal mats may compete with microphytobenthos (MPB) for light and nutrients and, due to their position in the canopy, have a negative impact on MPB biomass. We tested this effect by conducting a meta-analysis of prior experiments, as well as a comparative survey, and a macroalgal-removal manipulation in the coastal lagoons of the Virginia Coastal Reserve (VCR) on the eastern shore of Virginia (USA). In all cases, MPB biomass was estimated using benthic chlorophyll as a proxy. While prior individual studies documented impacts of macroalgae, when effect sizes were averaged across studies, there was no consistent effect of macroalgal biomass on MPB biomass. In the VCR, a non-native red macroalga, Gracilaria vermiculophylla, dominates intertidal mats and attains high biomasses at some sites. Nevertheless, MPB biomass was unrelated to macroalgal mass based on a survey of mudflats. Further, when macroalgae were experimentally manipulated on a mudflat using a before and after impact design, there was no change in MPB. Based on the meta-analysis, survey, and manipulation we conducted, macroalgal mats do not have a generalizable effect on MPB, interactions seem context-dependent, and in the VCR, the effects on MPB appear neutral. This finding is important given the significance of MPB in supporting food webs and other estuarine ecosystem functions, as well as the increasing frequency and intensity of macroalgal blooms.     

Phytoplankton Spatial Variability in the River-Dominated Estuary,Apalachicola Bay,Florida

In shallow estuaries with strong river influence, the short residence time and pronounced gradients generate an environment for plankton that differs substantially in its dynamics from that of the open ocean, and the question arises “How is phytoplankton biomass affected?” This study assesses the small-scale spatial and temporal distribution of phytoplankton in Apalachicola Bay, a shallow bar-built estuary in the Florida Panhandle. Phytoplankton peaks were characterized to gain insights into the processes affecting spatial heterogeneity in biomass. Chlorophyll a (Chl a) distribution at 50-m spatial resolution was mapped using a flow-through sensor array, Dataflow©, operated from a boat that sampled four transects across the bay every 2 weeks for 16 months. Chl a peaks exceeding background concentrations had an average width of 1.3?±?0.7 km delineated by an average gradient of 3.0?±?6.0 μg Chl a L^?1 km^?1. Magnitude of E-W wind, velocity of N-S wind, tidal stage, and temperature affected peak characteristics. Phytoplankton contained in the peaks contributed 7.7?±?2.7% of the total integrated biomass observed along the transects during the study period. The river plume front was frequently a location of elevated Chl a, which shifted in response to river discharge. The results demonstrate that despite the shallow water column, river flushing, and strong wind and tidal mixing, distinct patchiness develops that should be taken into consideration in ecological studies and when assessing productivity of such ecosystems.     

Creep behaviour of intact and remoulded fibrous peat

This paper presents the creep behaviour of intact and remoulded specimens of fibrous peat obtained from a field site near Anzac, Alberta, Canada. The creep behaviour was investigated by means of long-term drained and undrained triaxial tests. The development of volumetric, axial, and undrained axial strain and strain rate during drained and undrained creep tests under variable stress conditions is presented. The stress _– strain _– strain rate (p′–ε _v–\(\dot{\varepsilon }_{\text{v}}\)) relationship is found to be unique for different stress and loading durations. The p′–ε _v–\(\dot{\varepsilon }_{\text{v}}\) relationship is analysed and represented by creep isotaches. The applicability of different creep models developed for normally consolidated clay is discussed and applied to define the development of creep strain in fibrous peat under varying isotropic and deviator stresses. The secondary consolidation coefficient for evaluating the volumetric strain rate of peat is found to be applicable with some limits. The drained creep behaviour of remoulded peat specimens differs from the behaviour shown by Shelby tube specimens, whereas the undrained creep behaviour in remoulded and Shelby tube specimens is similar.     

An orthotropic interface damage model for simulating drying processes in soils

The study of drying process in soils has received an increased attention in the last few years. This is very complex phenomenon that generally leads to the formation and propagation of desiccation cracks in the soil mass. In recent engineering applications, high aspect ratio elements have proved to be well suited to tackle this type of problem using finite elements. However, the modeling of interfaces between materials with orthotropic properties that generally exist in this type of problem using standard (isotropic) constitutive model is very complex and challenging in terms of the mesh generation, leading to very fine meshes that are intensive CPU demanding. A novel orthotropic interface mechanical model based on damage mechanics and capable of dealing with interfaces between materials in which the strength depends on the direction of analysis is proposed in this paper. The complete mathematical formulation is presented together with the algorithm suggested for its numerical implementation. Some simple yet challenging synthetic benchmarks are analyzed to explore the model capabilities. Laboratory tests using different textures at the contact surface between materials were conducted to evaluate the strengths of the interface in different directions. These experiments were then used to validate the proposed model. Finally, the approach is applied to simulate an actual desiccation test involving an orthotropic contact surface. In all the application cases the performance of the model was very satisfactory.     

A review on missing hydrological data processing

Like almost all fields of science, hydrology has benefited to a large extent from the tremendous improvements in scientific instruments that are able to collect long-time data series and an increase in available computational power and storage capabilities over the last decades. Many model applications and statistical analyses (e.g., extreme value analysis) are based on these time series. Consequently, the quality and the completeness of these time series are essential. Preprocessing of raw data sets by filling data gaps is thus a necessary procedure. Several interpolation techniques with different complexity are available ranging from rather simple to extremely challenging approaches. In this paper, various imputation methods available to the hydrological researchers are reviewed with regard to their suitability for filling gaps in the context of solving hydrological questions. The methodological approaches include arithmetic mean imputation, principal component analysis, regression-based methods and multiple imputation methods. In particular, autoregressive conditional heteroscedasticity (ARCH) models which originate from finance and econometrics will be discussed regarding their applicability to data series characterized by non-constant volatility and heteroscedasticity in hydrological contexts. The review shows that methodological advances driven by other fields of research bear relevance for a more intensive use of these methods in hydrology. Up to now, the hydrological community has paid little attention to the imputation ability of time series models in general and ARCH models in particular.