Comparison between cell-centered and nodal-based discretization schemes for linear elasticity

In this paper, we study newly developed methods for linear elasticity on polyhedral meshes. Our emphasis is on applications of the methods to geological models. Models of subsurface, and in particular sedimentary rocks, naturally lead to general polyhedral meshes. Numerical methods which can directly handle such representation are highly desirable. Many of the numerical challenges in simulation of subsurface applications come from the lack of robustness and accuracy of numerical methods in the case of highly distorted grids. In this paper, we investigate and compare the Multi-Point Stress Approximation (MPSA) and the Virtual Element Method (VEM) with regard to grid features that are frequently seen in geological models and likely to lead to a lack of accuracy of the methods. In particular, we look at how the methods perform near the incompressible limit. This work shows that both methods are promising for flexible modeling of subsurface mechanics.     

Computing derivative information of sequentially coupled subsurface models

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.     

Evidence for the occurrence of hydrogeochemical processes in the groundwater of Khoy plain,northwestern Iran,using ionic ratios and geochemical modeling

Rapid population growth, industrialization, and agricultural expansion in the Khoy area (northwestern Iran) have led to its dependence on groundwater and degradation of groundwater quality. This study attempts to decipher the major processes and factors that degrade the groundwater quality of the Khoy plain. For this purpose, 54 groundwater samples from unconfined and confined aquifers of the plain were collected in July 2017 and analyzed for major cations and anions (Na, K, Ca, Mg, HCO₃, SO_4, and Cl), minor ions (NO₃ and F), and Al. Magnesium and bicarbonate were identified as the dominant cation and anion, respectively. Several ionic ratios and geochemical modeling using PHREEQC indicated that the most important hydrogeochemical processes to affect groundwater quality in the plain were weathering and dissolution of evaporitic and silicate minerals, mixing, and ion exchange. There were smaller effects from evaporation and anthropogenic factors (e.g., industries). Results showed that the high salinity of the groundwater in the northeast area of the plain was due to the high solubility of the evaporitic minerals, e.g., halite and gypsum. Reverse ion exchange and the contribution of mineral dissolution were more significant than ion exchange in the northeastern part of the plain. Elevated salinity of the groundwater in the southeast was attributed mostly to reverse ion exchange and somewhat to evaporation.     

The stratigraphical signature of the Anthropocene in England and its wider context

The Anthropocene deposits of England, here regarded as those formed after ～1950 CE, are now extensive, take various forms, and may be characterized and recognized by a number of stratigraphic signals, such as artificial radionuclides, pesticide residues, microplastics, enhanced fly ash levels, concrete fragments and a novel variety of ‘technofossils’ and neobiotic species. They include the uppermost parts of both ‘natural’ deposits such as the sediment layers formed in lakes and estuaries, and more directly human-made or human-influenced ones such as landfill deposits and the ‘artificial ground’ beneath urban areas and around major constructions. ‘Negative deposits’ include the worked areas of quarries and regions such as the English Fenland, where thick peat deposits have ablated to leave a strongly modified underlying landscape, and extend beneath into the subterranean realm as mine workings, metro systems and boreholes. The production of these is still rapidly increasing and evolving in character, while the early signs of global change, such as warming, sea level rise, and modifications to biotic assemblages, are beginning to further modify the emerging geology of this new phase of Earth history.     

Surface Geophysical Methods used to Verify the Karst Geological Structure in the Built-up Area: A Case Study of Specific Engineering-Geological Conditions

This article presents a research study of complex limestone karst engineering-geological conditions in the municipality Valaská near Banská Bystrica in Slovakia. The aim of the study is to demonstrate the impossibility of spatial identification of cave spaces using surface geophysical methods due to the specific engineering-geological conditions of a thick surface layer of anthropogenic fill containing highly heterogeneous anthropogenic material. Its maximum thickness is 3 m. Another specificific condition of the study area is its location in the built-up area, due to which the applicability of geophysical methods was limited. The article contains methodological recommendations to be used in analogous geological conditions with karst structures topped with anthropogenic fill, which complicates the identification of cave spaces. The recommended solution herein is the identification of the cave system using underground mapping of the karst and its projection onto the surface for which surface geophysical methods have been combined.     

Small-strain behaviour of unsaturated silty clay: experiments and model interpretation

Acta Geotechnica - The experimental data dealing with the so-called small strain stiffness of soils are indispensable in developing and calibrating advanced numerical models. A literature review...     

Parametric rolling behaviour of azimuthing propulsion-driven ships

The study investigates the experimental and numerical analysis of the occurrence of auto-parametric rolling for large, high-speed pod-driven ships in waves. Considering unique design and performance targets, the aim here is to exploit susceptibility to auto-parametric rolling behaviour and to identify probable design and operational precautions. In order to achieve this aim, an existing non-linear time-domain software to simulate capsizing and other critical manoeuvring behaviours of slow- to medium-speed conventional and podded ships in waves is being enhanced for fast pod-driven vessels and then compared against the dedicated model test conducted in long-crested regular and random waves for a large, pod-driven containership model. This paper includes the presentation of current numerical modifications for pod-driven ships and the verification analysis.     

Seagrass as the main food source of Neaxius acanthus (Thalassinidea: Strahlaxiidae), its burrow associates, and of Corallianassa coutierei (Thalassinidea: Callianassidae)

Burrows of the thalassinidean shrimps Neaxius acanthus and Corallianassa coutierei are striking aspects in tropical seagrass beds of the Spermonde Archipelago, Indonesia. Burrow construction, behaviour, burrow type and associated commensal community were investigated to clarify the ecological role and food requirements of these shrimps and their commensals. Gut content analysis and stable-isotope data were used to unravel the food sources and the trophic interactions among the commensal community.Individuals of Neaxius acanthus were caught on Bone Batang Island. In narrow aquaria filled with sediment they constructed burrows resembling those found in the field. During burrow construction and maintenance only little sediment was brought to the surface, most was sorted and compacted to create a distinct lining. Maintenance work by single shrimps typically took about 5 min, after which the shrimp walked up to the entrance and rested for a similar period of time. There were no differences in behaviour between day and night. Intrasexual encounters inside the burrow were characterised by a high level of aggression and all resulted in one participant being driven out of the burrow. Intersexual encounters led to coexistence with both animals taking turns in burrow maintenance and guarding the entrance. Offered seagrass leaves were pulled underground, cut into pieces and eventually integrated into the lining. Burrows of Corallianassa coutierei resembled a deep U-shape. Chambers branching off halfway down and at the deepest point contained seagrass fragments. All steep parts of the burrow were lined similar to burrows of N. acanthus.No commensals were found associated with Corallianassa coutierei. However, burrows of Neaxius acanthus in the field typically contained a pair of shrimps, up to 8 individuals of the commensal bivalve Barrimysia cumingii and large numbers of gammarid amphipods. Other animals found associated with the burrow were the goby Austrolethops wardi, a palaemonid shrimp species and two species of tube-building polychaetes, one of which was also found as an epibiont on N. acanthus.Stable-isotope and gut content analyses indicate that the diet of Neaxius acanthus, its commensal Austrolethops wardi, and Corallianassa coutierei is mainly derived from detrital seagrass leaves, with a potential contribution of sediment organic matter and seagrass epiphytes. In contrast the isotopic signature of Barrimysia cumingii suggests the presence of symbiotic sulphide metabolism bacteria. This study underlines that, besides their interactions with the surrounding ecosystem, thalassinid shrimp burrows play an important role as a sub-habitat with a unique associated fauna.     

Numerical study of baroclinic tides in Luzon Strait

The spatial and temporal variations of baroclinic tides in the Luzon Strait (LS) are investigated using a three-dimensional tide model driven by four principal constituents, O₁, K₁, M₂ and S₂, individually or together with seasonal mean summer or winter stratifications as the initial field. Barotropic tides propagate predominantly westward from the Pacific Ocean, impinge on two prominent north-south running submarine ridges in LS, and generate strong baroclinic tides propagating into both the South China Sea (SCS) and the Pacific Ocean. Strong baroclinic tides, ∼19 GW for diurnal tides and ∼11 GW for semidiurnal tides, are excited on both the east ridge (70%) and the west ridge (30%). The barotropic to baroclinic energy conversion rate reaches 30% for diurnal tides and ∼20% for semidiurnal tides. Diurnal (O₁ and K₁) and semidiurnal (M₂) baroclinic tides have a comparable depth-integrated energy flux 10–20 kW m⁻¹ emanating from the LS into the SCS and the Pacific basin. The spring-neap averaged, meridionally integrated baroclinic tidal energy flux is ∼7 GW into the SCS and ∼6 GW into the Pacific Ocean, representing one of the strongest baroclinic tidal energy flux regimes in the World Ocean. About 18 GW of baroclinic tidal energy, ∼50% of that generated in the LS, is lost locally, which is more than five times that estimated in the vicinity of the Hawaiian ridge. The strong westward-propagating semidiurnal baroclinic tidal energy flux is likely the energy source for the large-amplitude nonlinear internal waves found in the SCS. The baroclinic tidal energy generation, energy fluxes, and energy dissipation rates in the spring tide are about five times those in the neap tide; while there is no significant seasonal variation of energetics, but the propagation speed of baroclinic tide is about 10% faster in summer than in winter. Within the LS, the average turbulence kinetic energy dissipation rate is O(10⁻⁷) W kg^{− 1} and the turbulence diffusivity is O(10⁻³) m²s⁻¹, a factor of 100 greater than those in the typical open ocean. This strong turbulence mixing induced by the baroclinic tidal energy dissipation exists in the main path of the Kuroshio and is important in mixing the Pacific Ocean, Kuroshio, and the SCS waters.     

ALMA observations of the circumstellar envelope around EP Aqr

Atacama Large Millimetre/sub-millimetre Array(ALMA) observations of CO(1–0) and CO(2–1) emissions from the circumstellar envelope of the asymptotic giant branch(AGB) star EP Aqr have been made with four times better spatial resolution than previously available. They are analysed with emphasis on the de-projection in space of the effective emissivity and flux of matter using as input a prescribed configuration of the velocity field, assumed to be radial. The data are found to display an intrinsic axisymmetry with respect to an axis making a small angle with respect to the line of sight. A broad range of wind configurations, from prolate(bipolar) to oblate(equatorial) has been studied and found to be accompanied by significant equatorial emission. Qualitatively, the effective emissivity is enhanced near the equator to produce the central narrow component observed in the Doppler velocity spectra and its dependence on star latitude generally follows that of the wind velocity with the exception of an omni-present depression near the poles. In particular, large equatorial expansion velocities produce a flared disc or a ring of effective emissivity and mass loss. The effect on the determination of the orientation of the star axis of radial velocity gradients, and possibly competing rotation and expansion in the equatorial disc, is discussed. In general,the flux of matter is found to reach a broad maximum at distances of the order of 500 AU from the star.Arguments are given that may be used to favour one wind velocity distribution over another. As a result of the improved quality of the data, a deeper understanding of the constraints imposed on morphology and kinematics has been obtained.